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PREFACE 


The ZAP laser analysis program was written for the 
Naval Research Laboratory, Washington, D.C., under contract 
NO0014-70-C-0341. The Principal Investigator was Dr. M. Troost. 
The code was programmed and debugged by J.E. Welch. 
Dr. J.H. Alexander contributed strongly to the geometrical 
analysis, in particular the Helix. Dr. A.G. Petschek defined 
the statistical analysis. Dr. J. Trenholme o1 NRL made many 
worthwhile suggestions for improving the program and contri- 
buted the second sample problem. The excellent cooperation 
by NRL personnel-throughout the contract period contributed 
significantly to the success of this contract. 
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ABSTRACT 


This report presents a detailed description of the ZAP 
computer program. It is used to determine the distribution 
of energy deposition within a laser assembly. The Program is 
written in FORTRAN IV. The calculation uses Tay tracing and 
Monte Carlo techniques to obtain a solution. The program is 
applicable to a wide Tange of geometries and a spectrum of 


wave lengths. 
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1. METHODS USED 


1.1 INTRODUCTION 


Several methods might be employed to calculate the 
distribution of the energy deposition within a laser assembly. 
Analytical calculations are impossible due to the complexity 
of the geometry. Numerical solutions of the relevant wave 
equations should be possible, but the computer time cost for 
three-dimensional calculations is prokibitive. Ray tracing, 
which is frequently used in optical problems, seems to be the 
most practical method available. 


1.2 RAY TRACING: 


In. ray tracing, rays may be selected systematically or. 
at random. In simple geometries, systematic selection of 
Trays is advantageous, as only a few are needed for a good 
sample. In complicated 3-dimensional geometries, the number 
of rays required for an adequate systematic sampling becomes 
prohibitively large, and randomly selected rays from known 
distributions will give more reliable results for the same 
amount of computer time. As many laser designs are of con- 
siderable complexity, the ZAP program uses randomly selected 
rays from specified distributions. Each ray is defined by 
the three directional cosines, the coordinates-of the start- 
ing point, and an intensity distribution as a function of 


a 


wavelength. . : oe 


"A problem any method of ray tracing in a design 
with a lar amount of internal reflection and no black 


absorbers, is to decide when to stop following a particular 


“Tay. As each ray loses its energy very gradually via partial 


absorption, it is impractical to follow the ray until it is 
completely absorbed; without a black absorber this would take 


not waste too much effort on then. 


1.3 MONTE CARLO AND RUSSIAN ROULETTE / 
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a Rs 
a very darge. number of reflections. If a ray is ignored only 
after its energy falls below a negligible level (e.g., 1% of 
its original strength), the paradoxical situatior exists where 
most of the tracing effort is used for rays that contribute 
very little to the final result. If, on the other hand, each 
ray is cut off when its energy is still significant (e.g., 25% 
of its original strength), larger errors will be made. A 
method is needed that does not ignore weak rays, but that does 


/ By randomly selecting the rays to be traced, the prob- 
lew is transformed into one that can be solved by Monte Carlo — 
methods. A technique frequently used in Monte Carlo solutions 
is the Russian Roulette Technique. This technique allows 
each ray to be cut off while its energy is still significant 
without ignoring the weak rays. Basically this method 
involves statistically eliminating some weak rays, while 


vs. Strengthening other weak rays for further tracing. In prac- 


titey. the Russian Roulette method works as follows. 


e 


A ray is. followed until its intensity falls below a 
certain predetermined level (e. ao% 50% of original intensity). 
At that point, it is decided randomly whether this ray will 


“survive. If it dies, it is ignored and no record is kept. 


If it survives, its intensity is multiplied by. @ factor which 
is the inverse of the probability of survival | in the random 
test. The intensity of the surviving ray is therefore always 
relatively high (e.g., above 50% of the original intensity). 
As each ray may several times pass through this procedure 
before being killed, even very weak rays are statistically 
represented in the sample. Partial- reflections may be treated 
by creating a new ray of appropriate intensity. A pass 
through the Russian Roulette routine will then decide whether 
the new ray survives and with what intensity. . 

2 


Ne 


| 
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Energy deposition is determined by calculating the 
exponential attenuation along the path of a ray. The deposi- 
tion is assigned to the segment through which the ray passed. 
The intensity is checked at every apEET SCE P60 of the ray 
with a boundary, and a game of Russian Roulette is played 
when necessary. Eventually each vay A gia be killed 
in a game of Russian Roulette. 


1.4 ACCURACY 


The accuracy of the energy absorbed in the various 
regions of the system is dependent on the humber of rays 


- traced. Several rums may be necessary to determine the number 


of rays needed for the accuracy required. 


1.5 MULTIPLE WAVELENGTHS 


Multiple wavelengths are handled in the ZAP program by 
assigning intensities for each wavelength to a Single ray. 
For this technique to be applied it must be assumed that the 
angle of refraction at a material interface is independent 
of wavelength, which is an approximation. The absorption 
coefficients and the percent reflected may vary with wave- 
length. 


As the ray is traced through the system, the amount 
absorbed for different “wavelengths varies. Therefore some 
wavelengths willbe destroyed by Russian Roulette while other 
wavelengths Survive. The ray is traced until all wavelengths 
have lost their intensity either by absorption or by Russian 
Roulette. 


The energy deposited in a material by a given ray is 
the sum of the contributions of the various wavelengths. 
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1.6 GEOMETRY SPECIFICATIONS 


A flexible method for describing the geometry of the 
system has been developed which allows simple input and fast 
computation. At the same time, it allows for great gener- 
ality in geometrical configuration. 


A problem common to many Monte Carlo codes is their 
extreme slowness when applied to complicated geometries. 
This slowness is caused by the need to check for possible 
intersections with all possible boundaries. The running time 
thus becomes roughly proportional to the number of boundaries 
in the system. It becomes quickly evident that some technique 
is needed which eliminates the necessity for checking all 


boundaries for each intersection. 


One possible solution to this problem is to divide 
space into subregions and-to check only for possible inter- 
sections with those boundaries that exist in the subregion 
and. with the boundaries of the subregion itself. Some 
modern Monte Carlo codes use this method, but they are often 
inconvenient to use because of the cumbersome geometrical 
definitions required for each subregion. 


Systems, Science and Software (S*) has developed a 
method which uses the subdivision method recursively. The 
geometry is specified by a greatly simplified geometrical 
input, as the subregions used are the natural geometrical 
“ shapes of the system. 


1.7 SEGMENTS 


Each part of space in the problem to be studied is 
divided into one or more segments. A single piece of mate- 
rial may belong to several segments or a segment may contain 
several materials in the form of subsegments. Two kinds of 


segments are defined: Mother segments and Daughter segments. 
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The Mother segments act as subregions and the Daughter seg- 
ments lie within Mother segments. When a photon lies in a 
particular Mother segment, only the intersections with boun- 
daries of that Mother segment and its Daughter segments need 
to be checked. The boundaries of the Mother segment shield 
the photon therefore from the details in the space outside 
the Mother segment. 


It frequently is also desirable to shield the photon 
from details in certain areas within the Mother segment. 
This may be done by enclosing these details in a Daughter 
segment and declaring this Daughter segment to be a Mother 
segment at a lower level. We now have a recursive definition 
of Mother and Daughter segments, where any Mother segment may 
be a Daughter segment of a higher level Mother and any 
Daughter segment may be a lower level Mother segment. 


The geometry of the complete system fits into one 
Mother segment. This Mother segment usually will have several 
Daughter segments that are Mother segments in turn. 


This recursive subdivision may be continued to 
desired depth. Any Daughter segment, that is not a Mother 
Segment, must consist of a homogeneous material. 


A Mother segment is a segment that is partially or 
completely filled by one or more Daughter segments. That part 
of space within a Mother segment that is not occupied by 
Daughter segments has the material properties of the Mother 
Segment. The volume of the Mother segment is the volume of 
the space within the Mother segment that is not occupied by 
Daughter segments. 


Whenever possible, Mother segments should have Simple 
shapes, e.g., rectangular parallelepiped, cylinder, sphere. 
A simple shape allows an efficient check for boundary cross- 


‘ 


ings and reduces Computer time. However, in ea. 
\ 
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Mother segment may have any shape, even though that may be 
expensive in computer time. 


Although the possible number of Daughter Segments 
within a Mother segment is unlimited, in practice this number 
should-be quite low (e.g., below 5) to reduce the hurber of 
boundary Crossing checks. After all, the Mother segment was 
introduced to reduce the number of potential boundary cross- 
ings. 


A Daughter segment must lie completely within its 
Mother segment. Each Daughter segment is therefore associated 
with only one Mother Segment. Each Daughter Segment may in 
turn be a lower level Mother segment. 


For each Mother Segment, all of the external boundaries 
must be specified. Boundaries of Daughters need not be speci- 
fied for the Mother. For each Daughter Segment, all of the 
confining boundaries must be specified. 


In specifying the boundaries of a given segment, it 
is necessary to know on which side of the boundary the seg- 
ment lies. For a plane, one must specify if the segment is 
on the origin or non-origin side of the plane. This is 
accomplished by specifying a positive boundary number for 
the origin side and a negative boundary number to indicate 
that the segment is on the non-origin side. For curved sur- 
faces, a positive boundary number means that the segment is 
inside the surface and a Negative boundary number indicates 
that the segment lies outside the boundary. 


1.8 SURFACE DEFINITIONS 


The geometric surfaces of the System are defined con- 
veniently for vector Calculations. In theory, any 
mathematically definable geometric shape can be used with the 
ZAP program by adding the appropriate subroutines to the 
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computer code. The definitions for the geometry types of 
most general interest are summarized below: 


1.8.1 Plane 


A plane is defined by the directional cosines of the 
normal to the plane and the coordinate of a point which lies 
on the plane (Fig. 2.1). 


1.8.2 Cylinder 


A cylindrical surface is defined by the directional 
cosines of the cylindrical axis, the coordinates of a point 
on the axis, and. the radius (Fig. 2.2). 

, 
1.8.3 Conic Cylinder 

A conic cylinder whose major axis is perpendicular to the 
Z axis ig. 2.4) is defined by the directional cosines of 
the major axis of the conic, the X-Y coordinates of the focus, 
and the /X-Y coordinates of two points on the conic surface. 
The two points must not be, symmetric relative to the axis, 
although one of them may lie on the axis. . 


Whether the conic represents a parabola, hyperbola, or 
ellipse is defined by the value of the eccentricity, which 
can be calculated from the data provided. 


|parabola: e=1 
hyperbola: e, >1 
ellipse: ey. < 1 
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1.8.4 Sphere 


A sphere is defined by the coordinates of the center 
of the sphere and the radius (Fig. 2.3). 


1.8.5 Half-cone 


A half-cone whose axis of Symmetry is parallel to the 
Z axis (Fig. 2.5) is defined by the coordinates of the vertex, 
and coordinates of a point on the surface of the half cone. 


1.8.6 Paraboloid of Revolution 


A paraboloid of revolution whose axis of Symmetry is 
parallel to the Z axis (Fig. 2.7) may be defined by specifying 
the coordinates of the vertex and the Z coordinate of the 
focus. — 


1.8.7 Prolate Elipsoid of Revolution 


A prolate ellipsoid of revolution whose axis of Symmet 
is parallel to the Z axis (Fig. 2.6 ) is defined by the coor- 
dinates of one of the vertices, the Z coordinate of the focus 
nearest the -chosen vertex, and the coordinates of a point on 
the surface. 


1.8.8 Hyperboloid of Revolution 


One sheet of a hyperboloid of revolution whose axis of 
Symmetry is parallel to the Z axis (Fig. 2.8) is defined by 
the vertex of the sheet of interest, the Z coordinate of the 
focus of the sheet of interest and the coodinates of a point 
on the surface. This point may lie on either sheet. 


1.8.9 Helical Tube 


A helical tube must have its major axis parallel to 
the Z axis. A point on the axis is chosen such that the 


_——— 
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perpendicular from the point to the centerline of the helical 
tube lies along the X axis and points from the point to the 
centerline in the positive X direction. The perpendicular 
distance from the axis to the centerline is specified and 

the radius of the tube must be specified. The only other 
parameter required is proportional to the distance between 
successive coils. 


1.8.10 Boundary Intersections 


In order to determine which boundary is next inter- 
sected by a ray, it is necessary to calculate along the 
direction of the ray the distance to all pertinent boundaries 
and choose the boundary which gives the smallest positive 
distance. Because of the Segmentation, it is not necessary 
to calculate the distance to all of the boundaries in the 
System, but merely to a subset of these boundaries defined by 
the current Mother and its daughters. This is accomplished 
as follows: 


1. When a photon is inside a Daughter segment, 
the following boundaries should be checked: 


a. The boundaries of the current Daughter 
segment. 


b. The boundaries of the Mother segment. 


2. When a photon is inside a Mother segment, 
but not inside any of its Daughter segments, 
the following boundaries should be checked: 


a. The boundaries of all Daughter segments | 
in this Mother segment. 


b. The boundaries of the Mother segment. 
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1.9 NEW SEGMENT DETERMINATION 


After a photon has Crossed a boundary, the Segment con- 
taining the photon must be determined. This requires check- 
ing each boundary of the Segment to determine if the Photon 
lies on the desired Side of the boundary. If the photon lies 
on the wrong side of 2ny boundary of a segment, it is not in 
that segment. This search is Conducted as follows: 


1. Is the photon stil] in the Mother Segment? 
If the answer is yes, 


2. Check each Daughter segment until one is 
found that contains the ‘photon. 


If the photon does not lie in any of the Daughter Segments, 
it lies by definition in the Mother Segment of space. If the 


see if this Daughter segment is a Mother segment at a lower 
level. If such is the case, the Procedure must be repeated 
Starting at (2) for this new Mother segment. 


If the photon lies outside the current Mother segment, 
the procedure should be repeated Starting with (1) for the 


next higher level Mother Segment. F 


The search thus continues until the final segment con- 
taining the photon has been found. This might involve Climb- 
ing several Mother levels and then descending several Mother 
levels. The algorithm is written recursively, and is not 
Concerned with the absolute Mother level. When a Photon is 
found outside the highest level Mother, which defines the 
boundaries of the complete geometry, an error Stop will occur. 


Ser ne 
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1.10 NORMAL TO THE SURFACE 


In order to calculate the angle of incidence and the 
angle of refraction for a ray intersecting a surface, it is 
necessary to determine the directional cosines of the normal 
to that surface at the point of intersection. The equations 
used to calculate the directional cosines of the normal depend 
on the geometry type of the surface. , ‘ 


1.11 ANGLE OF REFLECTION 


When a ray intersects a surface, part or all of the 
ray reflect: off the surface at some angle which must be 
determined. In the case of specular reflection, this angle 
is the same as the angle of incidence, but in a different 
direction. With a diffuse or Lambertian surface, the ray 
reflects with an angle chosen at random from a known distri- 
bution. : 


1.12 ANGLE OF REFRACTION 


When part or all of a ray is transmitted through - a 
real surface, the transmitted ray may be refracted. The direc- 
tional cosines of the refracted ray must then be calculated. 


1.13 PERCENT OF REFLECTJON 


Each boundary is assigned a reflection type which 


defines the method used to determine the percent of the ray 


that is reflected. “Each reflection type is designed to. 
handle a particular type of boundary that might be necessary 
in setting up a geometrical system. All reflections are 
assumed to be specular with the exception of diffuse reflec- 
tion. The following types of reflection are currently in the 
ZAP program. 
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1.13.1 Fresnelts Reflection for Dielectrics 


The surface of laser rods and lamps or wny other type 
of glass component will generally be PEE Cente by Fresnel's 
reflection. 


1.13.2 No Reflection 


A non-reflective or dummy boundary is used primarily 
for segmentation purposes. The material should be the same 
on each side of the boundary. Hence there is no refraction 
either. “4 \ 


1.13.3 Total Reflection 


Total reflection boundaries are used primarily for 
boundaries of symmetry, but can-also be used for other pur- 
poses where one wishes the impending ray to be totally reflected. 


1.13.4 Constant Reflection 


Rays are reflected with constant percent of the inten- 
sity going with the reflected ray regardless of material 
properties or the angle of incidence. 


1.13.5 Table Reflection 


Rays are reflected by a percent which is dependent on 
the angle of incidence. A table is provided which gives the 
percent as a function of angle, and interpolation in the table 


gives the percent to be used. 
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1.13.6 Function Reflection 


Any number of functions might be used to give the per- 
cent of reflection. One of particular interest is 


percent reflection = 1-A cos” 6 oe 


where 
- @ = angle of incidence 


A=a déefficient between 0 and 1. 


i 
1.33.7 Diffusion Reflection 


In the case of a diffuse surface (see angle of reflec- 
— : a 
tion), the percent to be reflected is considered to be a con- 
stant which may be chosen for the particular surface under 


consideration. 


oO 


1.13.8 Wavelength Dependent Reflection 

The percent of the ray to be reflected'may be dependent 
only on the wavelength, or may be dependent on the wavelength 
and the angle of incidence. ae 


1.14 SOURCE DEFINITION 
1.14.1 Surface Source : 

One type of source that is available is a surface 
souree for a cylinder whose axis is parallel to the Z axis. 
A random point is chosen on the surface of the cylinder, and 
a. ray is started from that point. The direction of the ray 
is determined by a subroutine that is provided by the user. 
This subroutine provides an angle 9 with the normal to the 
surface at the point, and an angle ¢ with the cylinder 
axis in th. .’ane tangent to the point. The distribution 


of the rays is therefore under user control, and can be 
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chosen to fit the model desired. From these two angles, the 
directional cosines of the ray can be calculated. An option 
provides for the rays to be directed outwards from the cylin- 
der surface or inwards from the cylinder surface. 


1.14.2 Volume Source 


If desired, a volume source may be used. For a volume 
Source, a starting point is chosen at random within the scurce = 
‘segment. A direction is then chosen in such a way as to give 
isotropic emission. This is accomplished by choosing another 
random point within a unit sphere having the starting point 
as the center. These two points uniquely define the ray and 
the directional cosines are easily calculated. / 


More than.one source may be included within a given 
System. For each source the number of rays to be created is 
input and an initial ‘energy distribution for the desired 
wavelengths must be provided. This distribution may be read 
in from cards or may be input by providing a subroutine to 
produce the desired results. 


1.15 ENERGY DEPOSITION 


When a ray passes through an absorbing material, 
energy is deposited in the material for each of the wave- 
lengths being considered. For each wavelength the ray inten- 
Sity is decayed exponentially, and the amount lost by the Tay 
is deposited in the material through which the ray has passed. 
The amount lost by different wavelengths is dependent on the 


distance traveled and the absorption coeffic'ent of the mate- eee 


rial for each wavelength. It is therefore hecessary to pro- 
vide an absorption distribution _as a function of wavelength 
for each material. | 
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1.16 LOGIC FLOW 


The flow of logic using the ZAP program is as follows: 


1. 


2. 


Generate a ray or continue with a leg of 


the previous ray. | 


Determine the number of! the segment con- 


taining the ray. | 


| 


Find the distance to the intersection point 
or points for each of the appropriate boun- 
daries and choose the boundary with the 


smallest positive distance. 
at a 


Deposit energy in the segment» through 
which the ray. has just passed. _ 


Determine the segment ante on the 
opposite side of the intersected boun- 
dary. 

Find the normal to the surface at the 
point of intersection. 


Calculate the angle of reflection, the 
angle of refraction and the percent of 
the ray to be reflected. Calculate the 
directional cosines of the reflected and 
transmitted rays. 


Play Russian Roulette with each wavelength 
of the reflected and transmitted rays. If 
the wavelength survives, double its inten- 
sity. If it does not survive, zéro its 

intensity. The same-random number is used 


for each wavelength. ; 
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9. Determine if the transmitted ray still has 
live wavelengths. If it does, start again 
with step 3. If it does not, go back to 


step 1. 


The above procedure is repeated until the desired nua- 
ber of rays have been traced to extinction. After all rays 
have been completed, the final segment and material results 
are printed, and the problem ends. 


a 
1.17 OTHER USES 


The ZAP program is a technique for solving optics 
problems of many types. Although it was developed for the 
analysis of laser efficie cy it is also a general optics pro- 


~~ gram that may be~applied to a wide range of problems. 


An example of another use is in the area of designing 
reflectors for spotlights. The problem would be set up with 
a reflector and spherical light source within the reflector. 
At some desired distance in front of the reflector a black 
nonreflecting plane wouid be placed, and the material behind 


a the plane would be divided up into segments. The final 


_absorbed energy distribution in the target segments would 

give a very effective measure of the focusing ability of the 
reflector. Many optics problems involving light sources and 
eventual light absorption may be solved with the ZAP progran. 


/ 
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\ 2. THE PROGRAM 


mA 
2.1 PROGRAMMING PH'LOSOPHY 


The ZAP program'was developed on the UNIVAC 1108 and 
was written in basic FORTRAN IV and is as ‘machine independent 
as is possible. A version exists for the CDC-3800. 


Throughout the development of the program, every effort 
was made to make the program as general as possible. In those 
cases that require special geometric shapes not available in 
the present code, new geometry routines may be incorporated 
in the program by adding the name of the geometry to a table 
and writing four subroutines to handle the tasks required for 
each geometric shape. Likewise, new types of reflection or 
ray creation may be added in the same manner. | 


As far as is possible, each subroutine or function in 
the cede handles a given task without regard to what has 


happened before the routine was called. In that way*these™™’ a aialacnil 


routines may_be, called fxom,.varieus.sections of the program 
to perform the same task for different purposes. Obviously 
this is not always possible, but a definite effort was made 
in that direction whenever possible. One obvious drawback to 
this is the resulting large number of subroutines, some of 
which are so short as to be almost ridiculous. However, it 
was felt that the advantages far outweigh the disadvantages. 
Changing the way a particular task is performed usually means 
replacing a small subroutine rather than making changes to a 
Jarge subroutine. 


The remainder of this section presents the mathematical 
equations used for source generation, reflection and refrac- 
tion. The equations used for each of the allowable geometries 


are given next. 
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Z.2 SOURCES 


After a ray has been absorbed or eliminated by Russian 
Roulette, a new rily must be created. The intensity distribu- 
tion of the new rity is acquired from a table (ESOR(N, F)) that 
was set up at the beginning of the run, where N = the source 
number and F is the wave length index. 


‘The position and direction of the new ray depends on 
the type of source being used. 


Let 
x. 
i 
Y; = the coordinates of the Starting point 
of the ray. 
a 
\. 
ES a e \ 
$i 
Bs = the directional cosines of the ray. 
Vi 


2.2.1 Surface Source 


For a surface source the new ray is created at the sur- 
face of a cylinder whose axis is parallel to the Z axis. The 


Tay may be directed either inwards or outwards from the sur- 
face. 
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Let 

9) = an angle in the plane perpendicular 
to the cylinder axis 

8 = the angle of the ray with the normal 
to the surface 

o = an angle with the cylinder axis in 

: the plane perpendicular to the normal 

ZMIN = minimum z value for Starting point 
of ray 

4z = istance along the z axis for start- 


ing point of ray 


oY, = coordinate of a point on the cylin- 
der axis. 


RAND = a random number between 0 and 1. 


| aig Choose a random angle 8) O< ®, < 21 
| od ‘ x; =x + (4,\ + e) ® cos e, 
Pe 
7 | m= yy * (% ¢ e) * sin 6, 
J z, = ZMIN + Az |* RAND. 


In order to calculate the directional cosines of the 
ray it is first necessary tio find the directional cosines of 
+4 : 
the normal to the surface at the point (x;, Yi» z;). 


\ 


\ 


n 
6. = directional cosines of the normal 
to the surface at the starting 
\ point 
Yn 


* 
\ 
‘ 


_ The angles © and \ are determined by a special sub- 
routine. This subroutine must be provided by the user, and 
will be written to give the desired distribution of rays. 


\ 


‘| ’ > : I 
/ ; — 
0<o<5 


0 <@ < 21 “3 


7 1A 


a. =a cos @ + Sin 0 sin 
i n ny in ¢ 


Bs = (gees a sin 6 sin $¢$ 
at 
} 
Y; * COsct sin 8 
‘ 
If it is desired for the ray to point inwards from the 
surface, the signs of its directional cosines. (a;, B i? Y;) 
are reversed. 


| ‘A 
2.2.2 Volume Source ! a | 
' | < ‘} \ 
For a volume source the pew Tay is created at a random 
point within the source Segment ,and the direction of EBe Tay 


is chosen to give isotropic emisgion. 


.The random starting point ° Gs determined by creating 
a random point (x; Yi> Zs i? within | rectangular parallelopiped 
that encloses the source segment. en a check is made to 
determine if the point \lies within the segmeat. If not, new 


‘random points are chosen until one is found that does lie 


\ 


\ \ 


| "4 


' 


within the segment. 


‘ 
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To find the directional cosines, random components of | 
a vector are found. | 


A 
> 
ta] 
A 
_ 
—_ 


Ax = a random number (-1 < 


Ay = a random number (-1 « Ay < 1) 


A 
> 
N 
A 


1) 


Az = a random number (-1 < 


then . — Se 
R = Ax 
\ 


| If R <1, the random components are used. If not, other-sets 
or random components are determined until an appropriate set — 


aaa is found. 
ih Finally : ! 
oi : 
| = ‘ 
; : PY Se Ax/R 3 
B; = Ay/R 
: 5 = (Az/R , | Ne 
2.2.3 Ray Input ‘ K 
The third way to create new rays is to input the start- 4 
ing points and directional cosines\.of a Systematically selected \ 
set of rays. Then when a new ray is needed, it is’ available “ 
in a table. This penrune is primarily useful in | debugging. 


| 
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2.3 REFLECTION 


2.3.1 Specular Reflection 


When a ray intersects a specular surface, part or all 
of the ray is reflected at the same angle with the normal 
and in the same plane that contains the normal and the impend- 
ing ray. The directional cosines of the reflected ray are 
_'_ealculated as follows: 


Lai -2(a;0, * 656, * Yin ) 


ac =a. + ba 
i n 


r 
+ 
; BY. = B; bb, 
. + : 
Yy * Yi by, —e / 
where 
= 4 
B. = vhe directional cosines of the ray 
\ F 
Vi 
i : 
\ 
baal ; / ; { P} 
B. = the directional cosines of the normal 
t to the surface. 
i | 
Yn | 
. | \ | f 
\ : | 
\ | 
\ | . 
\ | 
22 ; | 
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That is not what thas been done in the ZAP code however. 
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2.3.2 Diffuse Reflection 


Diffuse (Lambertian) reflection is calculated, by giving 
the reflected ray & pseudo-random direction with a Lambertian 
distribution. Two angles (9 and ¢) represent the angle between 
the refiected ray and the normal to the surface, and the 
angle the projection of the ray onto a plane perpendicular to 
the normal makes with some line in that plane, re«pectively. 


By definition, a Lambertian distr*bution is such that 
the probability distribution P(6) varie. as cos 6, while 
P(¢) is constant. The problem then is to randomly choose 98 
and @ in a manner that gives the desired probability distribu- 
tions. One method of doing this is to choose a random 
¢( 0 < $< 2m) and a random cos” 6 (0.0 < cos? @ < 1.0). 


If the reflected ray is a unit vector, the probability 
distribution of the projection of that ray onto the plane 
perpendicular to the normal should be constant. This fact 
is the key to wie: technique used here. a 


With unifotn random numbers R, (0 < R. < 1} we get the 
diffusely reflected directional cosines. Of. the ray relative - 
to the normal aN) and two orthogonal axes (P, ® i in the perpen- 
dicular plane as follows: 


Let 


a 
" 


projection of ray onto p 


projection of ray onto ho] 


Yr = projection of ray onto XN 
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Then choose ry 


oR = 2R,-1 


Bp = Ria} 


If ap” + o. > 1.0, choose two new random numbers and calcu- 


late new values for a, and B,- 


Then 


ys Vio - (o,” = BR’) 


Normal WN has directional cosines a8 oY,: 


If , = 1.0: choose 


~ 


Yp = 0 


> 
Otherwise the directional cosines of P are 


l-a 2 


— eS 
SL SS 
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The directional cosines of 6 are 


Yq ~ “n’p ap8, 
Finally, the directional cosines of the ray are: 
a, = Gp%p \* Bet * Ypo, 
BL = apbp * Bebo * Yp8, 
1; = Sn¥n + BR¥Q * RY 


2.3.3 Percent of Reflection for Dielectrics. 


& When a ray intersects the surface of a laser rod, lamp, 
‘or any other glass component, the percent of the ray to be 
reflected must be calculated. The ZAP code uses Fresnel's 
equation for reflection to calculate this percentage. 


Le t 
cos a = cosine of the angle of incidence 
cos b = cosine of the angle of refraction : 
REF Ai = refractive index of the averted the ray 
is leaving. 
REF, = refractive index of the material the ray 


is entering. 
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then 


REF i 


REF > 


} c = a.a + iS + e 
/ os a | in BB. v7, | 


cos Pg la 5%, + BB, . Yr#'n! 


cos b - mcos a 2 


T = —_ - 
s cos 6 + m cos a 


cos a cox b - n(i - cos” al) 


refer ( 
P S\cos a cos b + n(1 - cos” a) 


‘ T + Tr oF 


Percent reflected = +5 P 


2.4 REFRACTION 


When part or all of a ray is transmitted through a 
surface separating different materials, the transmitted ray 


ds refracted because of the difference in material densities. 


The directional cosines of the refracted ray are calculated 
from the following equations. 


Let. 
orf 
Bor = the directional cosines of the 
refracted ray 
Yr 
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? 


the refractive index of the material 


the ray is leaving 


the ‘refractive index of the material 
the ray is entering 


the directional cosines of the normal 
to the surface at the point of inter- 


section . 


the directional cosines of the impend- 


ing ray © 


la REF Ai 
REF 12 


c= a.a_ + - ON 
in B56. Yaivn 


b= ac(-1 + Vi + (1-2)/a2e? } 
Oi = aa; + ba 
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2.5 DEBUGGING 
2.5.1 Printer Plots 

: A printer plot. of the geometry cut by a plane may be 
obtained for as many different planes as desired. In each 
case a rectangular section of the plane is defined as input. 
The geometry within this rectangle is the geometry that will 
be plotted. 


wey , 
THe technique used to generate the points that represent 
the geo tric boundaries is as follows: fake rays are created 
and traced across the plane. Each time an intersection with 
a boundary occurs, the coordinates of the intersection are 
noted, and a blank is plotted. The material number is plotted 
between intersections. Rays are not reflected, and are 
traced until they leave the area of interest. 


. The only major restriction in defining the rectangular 
Plotting area is that the rectangle must be completely enclosed 
by the overall mother region (region number 1). This is no 
real restriction since the overall mother can be made very 


large with no loss of computer time or accuracy. 


The input required to define the rectangular plotting 
area consists of three points Po» Py» and p>. Point Pp Will 
become the lower left hand corner of the plot and paint P} 
will be the lower right hand corner. Point P5 will be on the 
top line of the plot. . beat 


Fake rays are Created at even intervals along the lines 
that form the top and lefthand sides of the plotting area, 
and are traced perpendicular to the line on which they are 
<reated. ~The number of rays to be Created is determined in 
such a way that one Tay corresponds to each print line and 
one ray corresponds to each print column. 


\ 


NX = number of vertical rays 


NY = number of horizontal rays 


eas 


: , : i 
Once an intersection point (x,y,z) has been found, the 
coordinate system must be rotated and translated to correspond 
to the plotting area. The equations needed are: 


t, i V(1-*5)? ™ (v1-¥0)” Bs [=i-2o)" 


L. = length of line Po? 


directional cosines of PoP, 


Point p. must be found 


* = 8x (*2-%0) * Ay (¥2-%0) * 8.2 (22-29) 


é a 
“2° #3 : 4x2 
then 
ty = (x2-x9)” * (72-¥9)" + (22-29) ” 


£. = length of line PoP 


29 


Pd . 
3SR-405 


X5-xX 
: 2 "0 
a = 
L 
mY 
et 
Y2-Yo $ vs 2 
ay = ca directional cosines of PoP> 
| 
Zaz 
} 4y2 ™ Z o | 
nm ~ < 


Rays are created at equal intervals along the top and 
lefthand side of the plotting area. 


NX = 130 = 


number of print positions to be used 
As. = £./(N -1) = interval along horizontal direction. 
As, = 104s, = interval along the vertical direction 


y= 2 /As| + 1.5 
ne ea ety 
Horizontal ray number n has the starting coordinates: 


X= £2 - (mn - .5) AS ayy 


Ye = ¥, - Oe -5)As a 


Z. = Za - (n ig -5)As_a 


y yz 
and directional cosines 
B=a 
xy L 
: ies 4x2 
30 
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Vertical ray number n has the starting coordinates: 
X= xX, + (n-1)As, a. 
7, yy * (n-1)As,a,y 
27 2+ (n-1)As, a 


and directional cosines 


\ 


a= -a 
B = -a 
Y* 7%, 


For each intersection point x,y,z the plotting coordi- 
hates are calculated as 


od 
X= a, (x-x,) + a,y(¥-¥9) + a, (2-29) 


Y= ayy (x9) + yy (-¥0) + ay2 (2-29) 


2.5.2 Random Errors 


The random errors in the results of the ZAP program are 
due to statistical variations in the number of rays passing 
through a particular segment and due to the variation in the 
amount of energy deposited by each ray passing through this 
segment. The latter variation is due to variations in the 
path length through the segment and variations in the ray 


' energy upon entry into the Segment. It turns out that the 


variation in the number of rays usually contributes to the 
largest fraction of the random error. 


PX oe ene 
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The random error is displayed in the output as a per- 
centage error of one standard deviation. When the number of 
rays contributing to a segment is less than four, S228 is 
printed to indicate the lack of data. Large errors are by 
their very nature unreliable as their calculation is based 
on unreliable results. An error larger than 20% indicates 


unreliable results, where an error below 5% indicates usable 


results. 


The percentage error is calculated according to the 


equation presented below: 


ER; = 
Le 
where 
ER, Percentage error in segment i for 
one standard deviation 
N Total number of rays in the problem 
E Energy deposited in segment i by a 
ray ; 
> E Total energy deposited in segment i 


2.6 GEOMETRY 
2.6.1 Equations ‘for a Plane | 


A plane is defined by the directional cosines of the 
normal to the plane and the coordinates of a point that lies 
on the plane. An additional piece of information needed in 
many of the equations is the distance from the origin to the 
plane. The directional cosines of the normal must be such 
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Ste 


\ that the normal vector points from the origin to the plane. 
The calculation of the distance from the origin to the plane 
3 \ aye : 
- provides a check on the direction, since 


e distance will be 
mpeative if the vector points in the wrong 


direction. 
" 
Define: ; \ 
ue | : 
a 
| By = the directional cosines of the normal 
' to the plane 
r} 
' 
\ | 
‘ ‘ \\ xX, f 
al i 
me ® | Yy = the coordinates of a point on the 
ie | plane 
2 
\ 
a | 
\ 2.6.1.1 Distance from the Origin to the Plane 
ge A dQ yxy + BY, + 1% 
ri 
\ f, : 
‘yeatel2 Intersection of a ray with a plane: 
, 


| Let. 
| \ L = distance along the ray to the point 
| of intersection 
| 


| B- = directional cosines of the ray 
| 
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| x. 
i ; 2 
| ' | E: 
\ | t — i q 
: Yi = coordinates of the starting point | a 
of the ray a 
ig ‘ 
z i a! ! 
f 
: 


The point of intersection (x,y,z) can then be calculated: © | 


x Sx. + La. 
eet aa 0 

+ Ee R 
yY*¥ LB ‘ | 


z= 2; + Ly; / 


2.6.1.3 Normal to the Point .of Intersection ‘ 


In the case of a plane, the normal to the point of 
intersection does’ not require an actual calculation, since 
the normal to aa plane is the same everywhere and is Pare, 


of the definition of the plane. Hence: 


“ ” % 
8 By 
vn ~ Vy 
‘ 
34 
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Qos 
See 
A w 
. n 
Y=] a 
> i 
So 
wn 


Serres | 
- 
a 


oe 
a 


tao 2.6.1.4 Position of Point Relative to Plane 


: In determining if a point is inside a given segment, 
it is necessary tto determine on which, side of a plane the 
point lies. This requires finding the perpendicular distance 


from the plane to the point. 
: ‘ tee > 


Da ox + By t M7” ! 


i 
\ 


i, 'y. = the coordinates of the point 
i" | 


\ : a ‘4 1 


A. cylindfical surface is defined by the directional 


t 
| 
| Ly 2.6.2 Santioae for a Cylindrical Surface 
; 
Ld ' cosines of the cylinder axis, the coordinates of a point on 


/ 


the axis, and the radius of the cylinder. 


4 


ens _ Define i 


By “ directional cosines of cylinder axis. 


Y;"| = coordinates of a point on the axis. 
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Then 
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d, = the radius of the cylinder. 


: 
‘ 


Intersection of a Ray with a Cylinder 


Let 


“Ax 


Zz 


‘a 


= distance along the ray to the poin: 
of intersection. r 


= directional cosines.of the ray. 


= coordinates of the starting point 
of the ray 


0 * <p. + We aE gs a 
a has : Bi8y. ' YG¥y - ; 


a; Ax ” B, Ay a ¥;Az og 


f 
a, Ax + B Ay + e 
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(aay? = Ax? + Az? + Az? 


g 


(22) (xx) - yy t f Gyy-22*xx)? - (1-xx)?(aa?-22?-d, 2) 
L = 


a 


2.6.2.2 In or, Out of a Cylinder | 


_ To determine if'a point lies inside a cylinder, it is 
necessary to calculate the distance from the point to the 
cylinder axis and compare this distance with the cylinder 


radius. ; 
+ 
Let 
D = distance of point from cylinder axis 
x. oe 
i 
\ 
5 = coordinates of a point 
z. 
i 
Then 
Ax *s es con Vv 
- = 2572) - - 
salt * = ax? + by? + bz? 3 (a, 4x - By 4y + v,tz)? 
s ~ . : _ ae. “FE 
if p? < d’ point. lies inside cylinder. \ 
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Let | 
x 
y = @ point on the cylinder surface 
iz , 
a 
é Bo = the directional cosines of the normal 
é to the cylinder at the point (x,y,z) 4 
\ 
¥; 
Calculate 
Ax = X-Xy 
ay = ye 
Az = Z-Zy a 
C= a, Ax + By Ay + yz. 
Ax-C 
ig — bas 
n 
s Ay-CB, 7 
B= 7 = 
* Az-Cy, : \ 
38 
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2.6.3 Equations for a Sphere | 
: 2 } 


A sphere is defined by the coordinates of the center 
and the radius. 


Define: 
* 
yy = the coordinates of the center of 
the sphere 
1 
d, = the radius of the sphere 


\ 


& ' / 
2.6.3.1 Intersection of a Ray with a Sphere 


Let 


L = distance along the ray to the point 
of intersection 


a. 
* 2 ' 
Bs = directional cosines of the ray 
Vi 
x- 
y; = coordinates of the Starting point 
of the ray 
Zz. 
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then 


L = -B + ¥ B’-¢ 


- 2.6.3.2 In or Out of a Sphere 


To determine if a point is inside a sphere, one must 
calculate the distance from the point to the center of the 
sphere and compare this with the radius. 


Let y 


a DD = distance of point from center of sphere 
*i 
Yi ) = coordinates of the point. 
! 
zZ- 


‘ »? . (x,-x,)’ a (yi-%,)° + (2,-z,) 


If p? < ay, the point is inside the sphere. 
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Figure 2.1 A piane whose. normal is not parallel to any of the 
Major axis. _ 
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Figure 2.2 Cylinder whose axis is not parallel to any of 
the major axis. 
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2.6.3.3 Normal to a Sphere 


Let 
x 
y = a point on the surface of the sphere 
z 
A 
Bo = the directional cosines of the normal 
to the surface at the point (x,y,Z) 
& 
Yn 
Calculate, 


a7 (7% Vay 
nt (yy) /4y 


' 

‘ 

DW 
tt 


; - = los Aga 
x : Ta OKI 


ve 


2.6.4 Equations for. a Conic Cylinder 


The conic cylinder is limited to one whose conic axis 

‘is perpendicular to the Z axis and: whose. cylindrical | axis is 

parallel to the Z axis. In this case, only the x- y plane. must 

be considered in defining the surface. The conic is defined 

by the coordinates of thé focus, the directional cosines of 

the conic axis, and the coordinates of two ‘points which lie 

on. _the surface. In addition, the eccentricity and the distance 
from the focus to the directrix are calculated for use in the 
- equations. See Figure 2. 4. 
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Figure 2.3 Sphere 


‘ 


Figure 2.4 Conic cylinder with axis parallel to Z-axis 
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Define 
‘a 
\ k | 
= directional cosines of the conic 
By axis 
xy 
= coordinates of the focus 
as 
’ x} 
Yj 
; = coordinates of two points on the 
- surface 
x 
2 
- Y> 
-- ey s eccentricity 
dy = distance from focus to directrix 


2.6.4.1 Eccentricity and Distance from Focus to Directrix 


nt Vea)? + (yi? 


oa ¥(x,-x, + (y,-y,)? 


. 


ex 

a, =— 
1 Tr, 
3, = 71% 
“e | 
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ann" . 
[ vy o N 


¥2Y, 


2 tT, 


cos e a + 6,8, 

cos 8, = aa, + BL B : N 

ie, See Se { 
\ p 

k r,. cos o + rT, cos 8, 


TY) Tr, (cos 6, - Cos. 8,) 


i =a 
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2.6.4.2 Intersection of a Ray with a Conic Cylinder 


Let ‘ 
x. 
i 
Y; 4 =. coordinates of starting point of ray 
j 
. : x 4 
: a 
| i | 
‘ B. = directional cosines of tke ray 
~e Yj 
; , OL = distance along the ray to the point of 
intersection. 
46 


aa 


Ay = iM 


4 ee: 

eae 6 = 1 Vi 
a* : a; /cos >) 

Fd . ne i 

8° = B./cos 8 ! 
EE = a” Ax + B* hy 


‘FF = ay Ax + B) Ay 
\ 
GG =a, a” + By B 
2 


2 2 ne te 


BB” = Ax 


2 ect 
2 pat 
H = EE-ey GG(FF + dy) 


2 


2 
Q= BB“ -er(FF+d 


Pe sa 
oO 


‘cos 8 
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2.6.4.3 In or Out of a Conic Cylinder 


To determine if a point (x; , y;) lies inside of a conic 
cylinder it is necessary to calculate the distance from the 
focus to the surface in the direction of the point, and the 


distance from the focus to the point and then compare the 
two distances. 


Let 
a 
i 
= the coordinates of the point 
Yi 
D = the distance from focus to point 
R = distance from focus to surface 


then 


Ax = x5} 


1 
os ae ap 4 
D= ax? + by” 
a” = Ax/D 
B“ = Ay/p 
- dy 


if D < R, the point is inside the conic. 
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2.6.4.4 Normal to a Conic Cylinder 


Let 


33 = directional cosines of the norral 
to the surface at the point (x,y,z) 


then 


hey) 
‘ 

tt 
g 
SN 
o 
= 


cos § aya + B.8 


sin 9 = ws = an By 
— a, cos 6 - By sin € - a 8 
n 
1 + ef 7 ze, cos 6 
—_ ay Sin S$ + =y Eee 5 - eye, 
” 2 


1 + ey - Ze, cos § 
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2.6.5 Equations for a Cone _ 

A right circular half cone whose axis is Parallei to 
the Z axis is defined by the coordinates of the vertex and 
the coordinates of a point on the surface of the cone. If 
the point lies above the vertex, the half cone is a cup. If 


it lies below the vertex the half cone is a Cap (see 
Figure 2.5). 


Define: 


Y, = the coordinates of the vertex 
a 
* 
yy = the coordinates o a point on the 
surface 
aa | 
2 a 
_ (2, -2,) 
©” Gael - (pts ee 
(x, x, ) YY, (21-2, ) 
Qe. = 272 
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Figure 2.5 Cone with axis Parallel to 2 


-axis 
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2.6.5.1 Intersection of a Ray with a Half Cone 


Let 


— 


i 
Ye = coordinates of the starting point 
of a ray 
“i 
5 
B. = directional cosines of a ray 
aS 
L = distance along the ray to the point 


of intersection. 
Then 


Ax = X57 Xy 


ay = yi-y, 


Aor 
B = (1-4, Jy, az - d, (a ax + 8 ay} 


C = (1-d, Jaz? - d, (ax? + ay?) 


= Bt VB?-ac 
A 


+ 
> 


“ne 
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We must. now choose the smallest positive length that 
gives an intersection point in the proper half of the cone. 
Let 


@ 
t 


" i (2; ks y;L,) (e,) 


~~ (2, + y;L) (ey) 


i) 
I 


if e, 1S positive, L, gives an intersection in 
the proper half of the cone. 


if e is positive, L_ gives an intersection in 
the proper half of the cone. 


2.6.5.2 In or Out of a Half Cone 


To determine if a point (x; , Y;° 7;) lies inside of a 
half cone, the following procedure is weed, 


Ax = xix 
Ay = y.-y, 
Az = z57%y es 


EP = Az*e 


FP = Az? - dy (ax? + by? + tz") 


if (FP) (EP) > 0, the point lies inside the half cone. 
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2.6.5.3 Normal to a Cone 


Let 
a 
n 
a = the directional cosines of the 
normal to the surface at the point 
(x,y,z) 
Yn 
then 
Ax = X-Xy 
BY (= y-y, 
A. = z-Zy 
2 
R = a2 (ex + ay’) + (1-4, )? Az? 
a: d, Ax/R | 
By = dyay/R 


Y 


n (1-4) 4z/R 


2.6.6 Equations for an Ellipsoid of Revolution 


The ellipsoid of revolution must have its axis of 
revolution parallel to the Z axis. It is defined by the 
Coordinates of one of the vertices (x, , Yue 24); the coordinates 
of a point on the surface (x,, Y)> z,) and the Z coordinate 
of the focus (z,) nearest to the vertex chosen. See 
Figure 2.6. 
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Figure 2.6 Ellipsoid of Revolution with axis parzllel to 27 
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Figure 2.7 Parabcloid of revolution with axis 


Z-axis 
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Let 
f = 2 6° zy 
A restriction on the choice of the point (x, > Ye 24) is: 
2(2y-2,) £-(2y-2)? < (xy-x, 7 + (yi-%)° < 4£(2472,) 


The paraneters at and b are calculated as follows: 


w= (xy)? + (yy)? 


- 2 
Zz oy oo 


SQ = sign(f) * ¥or-£)7 + We 


n2 | 
b = (z*+£+SQ) /(s- “7 
a? = p* - (b-£)2 


2.6.6.1 Intersection of a Ray with the Surface 


Ax = X57X} 
by = ¥3-Y, 
&z = ca 9 
coy ee ew? 


a b 
a x+3 ty y; (4z-b) 
B = - —*— + a 
a 
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if B*-AC < 0, there is no intersection. 
Calculate »? 


L - BN B?-ac 


+ A 


The intersection of interest is the mininun positive 
value of L,. 


2.6.6.2 In or Out of an Ellipsoid of Revolution 


To determine if point (x5 0¥; 525) is inside of the sur- 
face calculate: 


Ll 


2 a ae 
ae (x5) + (v5-%,) _ (747%) 
PP) a = “wl oe ee a 
a 
if F(P) > 0, the point is inside of the ellipsoid. 


2.6.6.3 Normal to the Surface at a Point Ge wo) 


2 2 2 
(Nt Oo ey) 
a b 


The directional cosines of the normal are: 


Fi 


a = k 
n v4 
aR 
A “s Vv -~ Yy 
z a’R 


a7 
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Coordinate of the focus. if the focus is above the vertex 
the figure is a Cup. If the focus is below the vertex the 
figure is a cap. See Figure 2.7. 


Define: 


x 
% = the coordinates of the vertex 
2 
agli? 
Ze = the z Coordinate of the focus 
£ = Ze - zy ae 
call 


be) = directional cosines of the Tray 


58 . 


¥; = starting point of the ray 


Then 
C= 46(z5-24) - (xy-m)? = (vi-%E) 
jp Be 2f7; - a; («4-%,) - 8: (¥;-y,) 
; | A= -(a; + 87) 


if B? < Ac there is no intersection — 


L = 


i 
id 
Lg -B + VB?-Ac 
i + A 

| 


Choose the minimum positive value of L, 


2.6.7.2 In or Out of a Paraboloid of Revolution 


FOB) = 48(25-%) - (x-y,) - (75-%)? 


if F(P) > 0, the point (x; 5¥4223) lies inside. 


59 


3SR-405 


3SR-405 


2.6.7.3 Normal to a Point on the Surface 


Let 

x 

y = @ point on the surface 

z 
a 

n 
E. = the directional cosines of the normal 
to the surface at the point (x,y,z) 
Yn 
Calculate 


br ~] 


= ye-x,)? + (y-y,) + 4£? 


“ee 


2.6.8 Equations for a Hyperboloid of Revolution 


A hyperboloid of revolution whose axis is parallel to 
the 2 axis is defined by the coordinates of the vertex of the 
Sheet of interest (xp oy, 52), the Z coordinate of the focus 
of the sheet of interest (z¢) and the coordinates of a point 
on the surface of either sheet (x) 5¥)>2,)- See Figure 2.8. 
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Figure 2.8 Hyperboloid of Revolution 
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Let 


f = Zp7Zy 
A restriction on the choice of the point (x) 5¥) 2) ig: 


2 2 
4£ (21-2) <(%-x,)% + (71-Y,) 


The parameters a? and b are calculated as follows: 


Ww - (x,-x,) = (7,-%) 


SQ =, Sign (2#) * V(z0-£)2 + 2 


b= (z*+ sQ/( : 4) 
a? = (bee)? - p2 
2.6.8.1 Intersection of a Ray with the Surface 
Ax = x57 
Ay "Yay 


Az = Z5-Z, 
C = -1 - Ax*+py? s (az+b)? 
a ‘-b 
= YCAz+b) _ aax+3e 
oe ge ge 
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if B?-AC < 0 there is nc intersection. 


Calculate 


L. = B+ ¥B?-ac 
SN BAC 


I+ 


The intersection of interest is the minimum positive 
value of L, that will give an intersection on the desir-d 
sheet. 


2.6.8.2 In or Out of a Hyperboloid of Revolution 


To determine if point (Xj ¥; 275) is inside of the sur- 
face, calculate: 


Lies es 2.,.2 2 
| F(B) = -1 - 4X sy" , = 
LJ 2 


If F(P) > O and Az*f > 0, the point is inside. 


2.6.8.3 Normal to a Point on the Surface 
int on the Surface 


Let 


y = a point on the surface 


te 
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the directional cosines of the normal 
to the surface at the point (x,y,z) 


2.6.9 Equations for a Helix 


A helix whose main axis is parallel to the Z axis is 


defined by 


Yh 


coordinates of a point on the axis 
of the helix 


Pp 


Point on the helical wire in the 
positive x direction from the axis 
of the helix 


Parameter defining the distance 
between successive coils (D) 
D = |2m«| 


handed. 
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if «x < 0, coil 


is left- 


—. 


9p 


a 


The equations 
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radius from axis to helical wire 


radius of the helical tube. 


used are quite voluminous. They are derived 


in Appendix A.1. 
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3. INPUT - OUTPUT 


3.1 SETUP CONSIDERATIONS 


The use of any large generalized computer program re- 
quires a certain amount of knowledge of how the program works 
and requires a learning process the Same as with a piece of 
complicated machinery. The more general a program is, the 
more understanding and. experience is required to effectively 
make use of it. The ZAP code is no exception in this respect. 


In particular, a thorough understanding of the mother- 
daughter relationship and its consequences is very important, 
Since this is the basis of the ZAP technique. To try to set 
up complicated geometries without this understanding will 
invariably result in boundaries existing in areas where they 
are not supposed to exist. To cut off the extension of a 
Plane for example, requires that the area where the plane 
actually exists be enclosed in a Mother segment. Then the 
plane wiJi not extend outside that Mother segment. 


Unfortunately, with Monte Carlo techniques gross errors 
can occur in the setup of a problem without Causing any | 
indication that these errors exist. Boundaries may exist in 
areas where they are not intended without ever being obvious 
in the results. Therefore, extreme Care and thought should 
be put into the setup of a problem. To help spot a bad setup, 
test data points and geometry plots have been included in the 
program. 


The bad effects of boundaries extending into areas 


segues 


where they do not belong has perhaps been over emphasized in 
the above paragraphs. Regardless of the type o boundary 
under consideration, it is treated as a dummy boundary and 
has no ill effects if the same material exists on each side 
of the boundary. Conversely, it must be kept in mind that 
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if the same material or two different materials with the same 
index of refraction are separated by a boundary, that boundary 
will be treated as a dummy regardless of how it is defined. 


It is highly recommended that anyone Planning extensive 
use of the ZAP program should become as familiar as is 
possible with the inner workings of the code. 


3.1.1 Miscellaneous Rules 
~~ meeeous Rules 


1. Planes must not Pass through the origin 
(0,0,0). 


2. Each Mother Segment must also exist as a 
daughter segment (except for the Overall 
Mother segment). ; 


3. The overall Mother segment must be segment 
#1. 


4. Each Mother Segment must have at least one 
daughter. 


5. A Mother Segment must completely contain 
each of its daughters. Conversely a 
daughter must be completely contained by 
its Mother. 


6. Different materials that are to be Separated 
by a real boundary must not have the same 
refractive index. 


7. Planes defined for geometry plots must not 
coincide with. a real Plane in the geometric 
System. 


8. A segment must be “totally" confined by its 
confining boundaries . For example, an open 
ended cylinder does not qualify as a seg- 
Rent. 
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9. Any rectangular area to be Plotted must lie 
completely within the Overall Mother Segment. 


3.1.2 User Routines 
es _ Routines 


Provisions have been made to allow the user to provide 
subroutines to handle certain functions within the program. 
A description of these routines is given below as they apply 
to the functions desired. 


Srehic 2 sl Wavelength Dependent Absorption Coefficients 


In a multifrequency problem the user must provide a 
subroutine or subroutines that calculate and store the absorp- 
tion coefficients as a function of wavelength and material 
number. The name of the subroutine for each material is given 
"on the material data card and must be of the form: 


ABSX 
where X is a number from 1-to 9. 
The argument list for the subroutine is: 


SUBROUTINE ABSX (XLAM, FMAX , TSTG ,M) 


where 
XLAM = name of the array that contains 
frequencies. 
FMAX = the number of frequencies (integer) . 
TSTG = Storage area for absorption coeffi- 
Cients from F=] to F=FMAX. 
M = material number. 


This subroutine must determine and store the absorption 
coefficients for material M into the array TSTG(F), F-1,FMAX. 
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One subroutine may be provided for each material or 
the same subroutine may be used for different Materials. 


3.1.2.2 Wavelength Dependent Intensities 


In a multifrequency run, the ray intensity may Vary 
with wavelength. Therefore a subroutine to provide intensity 
as a function of wavelength must be provided. The name of 
the subroutine to be used for each lamp type is given on the 
lamp data card and will be of the form: 

: INTENX 
where X is an integer from 1 to 9. 
The argument list is: 


SUBROUTINE. INTENX (XLAM, FMAX , TSTG , LAMXO) 


where: 


: 


name of the frequency array 
FMAX = the number of frequencies (integer) 
‘am array to store intensity as a 


function of wavelength for lamp 
number L.AMNO . ' 


= 
en 
Sa | 
ir) 
" 


LAMNO = lamp number. 


This subroutine must determine and store ray intensi- 


ties for lamp number LAMNO into the array: ‘ 
TSTG(F),F=1,FMAX \ 


One subroutine may be-provided for each lamp number, 
or the same routine may be used for different lamp numbers. 


yy 
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3.1.2.3 Surface Source Distribution | 


When a surface source is used, a subroutine must: be 
provided that will give the angles used to calculate the 
directional cosines for a given ray. The name of this routine 
is given on the lamp data card and must be of the form: 


. ANGLEX 


where X is a number from 1 to 9. 
The argument list for the subroutine is: 


SUBROUTINE .ANGLEX (TH, PHI) 


where 


TH = an angle (6) with the normal to the 
surface. 


PHI = an angle (¢) in the plane perpendicular 
to the normal. 


This routine must provide vuiues for 6 .and @ where: 


0 


[A 


@ < 4/2 


0<6 


[A 


2a 


an 


The routine is called each time a new ray is created ( 
and should provide values of 9 and @ that correspond to some 
desired probability distribution. 


An example would be to provide a Lambertion distribu- 
tion. In this case the subroutine would generate a random ¢ 
and a random cos” 6. From cos” 6, 6 would then be found. 
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3.1.2.4 Wavelength Dependent Fraction of Reflection 


When the fraction of a ray's intensity to be reflected 
is dependent on the wavelength, a subroutine must be provided 
to give the fraction for each wavelength. The Name of the 
subroutine must be of the forn: 


REFLXX 
where XX is the reflection number of the boundary. 
There are two Separate cases that may be considered. 


1. Fraction dependent on wavelength only 


2. Fraction dependent on wavelength and 
angle of incidence. 


In case 1 the argument list is: 


SUBROUTINE .REFLXX (XLAM, FMAX , TSTG , DUM) 


where 
XLAM = name of wavelength array 
FMAX = number of wavelengths (integer) 
TSTG = array for storing fractions of 


reflection | - : 
DUM = dumhy argument. 3 


In this case the routine is called at the beginning of 


the problem and the information is stored for later use. The 


subroutine must determine and store the fractions of reflec- 
tion into the array: = 


—- 
_-TSTG(F) , F=1,FMAX. 
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| bow 
| In case 2 the argument list is: 


SUBROUTINE .REFLXX (XLAM, FMAX , TSTG, COSGAM) 
where 


name of wavelength array 


‘| FMAX = number of wavelengths (integer) 
TSTG = array for storing fractions 


COSGAM = the cos of the angle of incidence. 


~ 


In-this case the subroutine is called when an inter- 


section occurs and must calculate the fraction to be} reflected 


for each wav:: :length for the particular angle of incidence 
indicated. These fractions are then stored into the array 


“TSTG(F) , F=1, FMAX 


3.2 INPUT 


Input to the ZAP code is .in the form of data cards. 


/ Special ..keypunch forms have been designed for each of. the 
eae 


‘various types" 2. Anput cards, with the exception of the 


jeenerel Papi which” wsES a stmuare p BORTRAN data sheet. 


Oe eve 


alphanumeric formats. In the case of real cg Gu taeee input, 
either a blank or zero may be entered if the number is zero. 


All directional cosines entered are checked by _ the 
code and the sum of the Squares of the three directional 


j cosines must be within the range of 1.0 + 10 e If they are 
‘not, the program will terminate execution with a debug print. 


Care must be taken that no plane passes through the 

origin. If this happens, the program will terminate execu- 
/ 
tion. / 


‘ 
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General Input (three cards) 


3SR-405 


hie 


» 
Pear a 
Card #1 - 72 Character Comment or title Card 


Card #2 


Description 


FRAC = Base fraction. When a ray 


intensity drops below this fraction 


of its initial intensity, Russian 
Roulette is Played. Default number 
BD 


DZ = Distance along Z axis between 
Planes of continuity. If no Plzies 
of Continuity exist, this numbei is 
ignored. ¥ 


WAVMIN = minimum wavelength. (May 


_ be omitted for single wavelength.) 


WAVMAX = maximum Wavelength. If 
there is only one Wavelength, this 


umber is ignored. 
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Card 
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Columns 


= Number of wavelengths 


IOPRNT = Option for printing each 


Tay created. 


IDEBUG = Option for debug prints. 


0 = no 
1 = yes~ 


IR = starting random integer. 
Default = 1. 
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3.2.2 


For each special reflection a reflection data card must 
be provided. The reflection number must be greater than ten, 
since reflection numbers 1-10 are reserved for built-in 


Reflection Data 


reflection types. 


Entry 
Number 


1) 


2) 


Reflection number (greater than 
ten). 


Reflection type (left justified). 


CONST. = Constant reflection (frac- 
tion reflected will be input). 


FUNCT. = Function of the form: frac- 
tion reflected = 1.0-A(cos 6)? where 


A is input and cos 0 is the angle 
made by the impending ray with the 
normal to the surface. 


TABLE. = An input table of ten frac- 
tions which. indicate the unt 
reflected for different gles with 
the normal. Angles are ftom 0° to 
90° in 10° intervals. 


DIFFUS = Diffuse reflection. Frac- 
tion to be reflected is input. 
WAVCAL = Reflection dependent on 
wavelength. The fraction |of each 


wavelength to be reflected is calcu- 
as a function of the wavelength and 


the angle of incidence. A subroutine 


REFLXX must be provided where XX is 
the reflection number. Example, 
reflection number 12 means that the 
subroutine is REFL12. (See section 
on User Subroutines for details.) 


WAVINP = Reflection dependent on 
wavelength only. A subroutine 
named REFLXX must be provided where 
XX is the reflection number. This 
subroutine is called immediately at 


the time this card is Tread, and fills 


an array with wavelengt dependent 
fractions of reflection\to be used 
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y 


Integer 


Alpha 


later. (See section on User Subroutines.) 
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Description Format Columns 


2-4 


6-11 


} 


ON AE AA NED + ert SOI SIN Ne oa 


‘For CONST load fraction to be re- Real 
flected. 


For FUNCT load value for A. 


For DIFFUS load fraction to be 
reflected. 


For TABLE load fraction for 0°. 


For WAVCAL or WAVINP this entry is - 
ignored. 


3) 


For TABLE load fractions to be re- Real 
flected for angles 10°-99°. Blanks 
for other reflection types. 


4-12 


13-17 


19-71 


A blank card must follow the jast reflection data card. 
If no reflection data was loaded, the blank is still required. 


3.2.3 Material Data 


For each material a material data card must be pro- 


vided. 


Description 


Material number 


Name of routine to generate absorp- 
tion coefficients as a function of 
Wavelength. - (See section on User 
Routines). If this is a Single 
wavelength run, this entry is 
ignored. 


Alpha 


‘ 


Absorption coefficient (length +) 
If an input routine was specified 
above, this number is ignored. 


Refractive index (aust be non-zero) 
Lazing wavelength (for multiple 
wavelength runs only). This is an 
optional input value. 


39 character description to be 
printed Alpha 


A blank card must follow the last material data card. 
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3.2.4 Boundary Data 
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For each boundary a boundary data card must be provided. 


The form of the data on the card depend:: on the geometry type. 


feat estan 


Integer] 2-4 


| Entry 


Number Description 
Boundary number 


Reflection number 


1 = total transmission, no 
refraction 


2 = total reflection. May or may 
not represent boundary of 


Symmetry. 
3 = Fresnel's reflection for dielec 
trics. 
4 = Top boundary of continuity. 
S$ = Bottom-boundary of continuity. 


6-10 = Reserved for further built-in 
reflection types. 


1-S0 = Special reflection types defin- 
ed on reflection data sheet. 


Geometry type (left justified in 
Columns provided on data sheet). 


CYLIND = cylinder 

PLANE. = plane 

CONIC. = conic 

SPHERE = sphere 

CONE), . = half-cone with axis par- 
| allel to Z axis 


pamab: = paraboloid of revolution 
with axis parallel to Z 
| axis. ; 
HYPERB = hyperboloid of revolution 
With axis parallel to the 
Z axis. 
ELLIPS. = ellipsoid of revolution 


whose axis is parallel to 
—— the Z axis. 


HELIX 


helical tube with axis pa 
allel to the Z axis 


77 


serene ree 


3SR-405 


Boundary Data (continued) 


Parameters for the particular geometry 
type being described: 


PLANE: 


Directional cosines for the nor 
mal to the plane. Normal must 
point from origin toward plane. 


Coordinates of a point on the 


0-11)--Not used’ - 


YLIND: | 
Directional cosines for cylin- 
drical axis. 


Coordinates of a point on the 
cylindrical axis.. 


R =Radius of cylinder 
-- Not used 


) a aeecenona! cosines for axis. 


B 
Xe)Coordinates of focus. 
Yeo) 
X1\)Coordinates of two points on 
Y: (the conic surface. These two 
X2(points must not be symmetrical 
Y2/with respect to the axis. 


Coordinates of center of sphere 


Radius of sphere 
vot_used 


Boundary Data (continued) 


+ | Entry 
Number 


CONE: 
4) xX, 
5) Y, 
6) Z, 
7) x 
8) Y 
9) z 

0-11) -- 

PARAB: 
4) x, 
5) Y, 
6) Z 
7) Z, 


HYPERB 
4) x, 
5) Y 


8) Y 
9) Zz 


10) Zr 
ll) -- 


4) xX 
5) Yy 


7) X 
8) Y 
9) Zz 


6) zy 


Description 


| coordinates of the vertex 


surface of the half-cone to 
be considered. 
Not used. 


| surface of of a point on the 


Coordinates of the vertex 


= Z coordinate of the focus 


Coordinates of vertex for 
sheet to be considered. 


soordinates of'a point on the 
surface. This point may be 
Jon either sheet. 


= Z coordinate of focus for 
sheet to be considered. 


Not used. 


ELLiPS: 


Coordinates of vertex 


6) Zy 


Coordinates of a point on the 
{ surface of the ellipsoid. 
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reruns | extant 


L7=22 
24-29 
31-36 


38-43 
45-50 
52-57 


17-22 
24-29 
31-36 


17-22 
24-29 
31-36 


38-43 
45-50 
52-57 


59-64 


17-22 
24-29 
31-36 


38-43 
45-50 
52-57 


Boundary Data (continued) 


Entry "ye 
Number Description 


ELLIPS (continued) 


10) Ze = Z-coordinate of the focus 
; nearest the chosen vertex. 
11) -- Not used. 


Restriction: 
A<B<C 
where ; 
2 
A= 2(Z-2,)(Z,-Z,) - (Z-z,) 


B= (X-x,)? + (v-¥,)? 
| C = 4(2¢-2,)(2-2,) 


| .. # HELIX 
; 4) x, Coordinates of a point on the 


5) y, { helical axis. The nearest 
[ point on the helical wire is 
- 6) Zh assumed to be in the positive 
» X direction from this point. 


| 7) « = Parameter associated with dis- 
| tance between consecutive 
coils. D= 2m. If K<0, coil 
| Ss is lefthanded. 


8) p = radius from axis to helical 
wire. 


9) a = radius of tube. 
0=11) Not used. 
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3.2.5 Segment Data 


For each segment, a segment data card must be provided. 


Description 


Segment number 


Material number Integer] 6-8 


,Boundary numbers indicating the con- | Integer]10-63 
fining boundaries of the segment. A 
positive toundary number indicates 
the origin side of a plane or inside 
of other geometries. A negative 
boundary indicates the non-origin 
Side of a plane, or the outside of 
other geometries. 


Segment volume (if known). 


A blank card must follow the last segment data card. 


3.2.6 » Structure Data 
For each mother-daughter relationship, a structure data 


card must be provided. 


7 


- Description 
Sa See ah * On 


Mother segment i Integer 


Bs 


Daughter segments bas | Integer 


A blank card nust follow the last structure data card. 
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| 3.2.7 Lamp Data ( 


For each lamp type a lamp data card must be loaded. 
| This data describes primarily the method to be used for creat- 
ing rays. Many different sources may refer to the same lamp 
ic type. 


Entry ; 
Number | | Description ; Columns 


Lamp number 


AN 


Lamp type (left justified) . \ 


- VOLUME = voluwe source 
OUTSUR 


cylindrical surface source 
_+ with rays directed outward 
from the surface. 


_ INNSUR 


cylindrical surface source 
with rays directed inward © 
from the surface. 


directional cosines and 
Starting coordinates are in- 
put for each ray. 


RAYINP 


Name of routine to ba & d for gener- 13-18 
ating intensity as a function of | 

wavelength. This routine will have a 
name INTENX where X is a number from 
1 to 9 (see user routines) and must 
be provided by-the user. If this is 
a single wavelength run, this entry 
is ignored. 

Name of routine to be used for calcu- 20-25 
lating the angle with the normal for 
rays created if this is a surface 
source. If not a surface source this 
entry is ignored. See user routines 
for description of this routine. 


Intensity of rays if this is a single] Real 27-37 
wavelength run. If:multiple wave- 


length, this entry is ignored. 


A blank card must follow the last Lamp Data card. 
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3.2.8 Source Data 


A source data card must be provided for each Source in 


-~the problen. 


yy 


Source number | 


Segment number or boundary number of 
source. If the 1 type designated 
is for a vol source, load a segment 
number. If tthe lamp type designated 
is for a surface source, load a boun- 
dary number. 


Lamp number - 


refers to a lamp 
data card. 


Integer 


Number of rays for this source, Integer 


Parameters which depend on the particu] Real 
lar lamp type designated. 


VOLUME : 


5) XMIN 
6) XMAX Minimum and maximum 
7) YMIN Values of X,Y and Z 
8) YMAX which describe a 
9) ZMIN rectangle within which 
10) ZMAX the source segment is con- 
taincd. 


INNSUR or OUTSUR: 


5-8) Blank 
9) ZMIN Minimum and maximum values 
10) ZMAX of Z for Creating rays. 

RAYINP: . 

Entries 5-10) are ignored. 


Entry ! 
Nuaber Description | Columns 


10-12 


14-17 


For RAYINP, this card is immediately followed by Ray 


data cards for the number of rays specified above. 


A blank card must follow the last Source Data card 
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3.2.9 Ray Data 


If rays are to be 
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input individually, one data card 


must be input for each ray to be created. Use source data 


sheets. 


Description 


Directional cosines of ray 


Starting point of ray 


3.2.10 Test Data 


Test data Cards ar 
for errors. This consist 
and indicating that a spe 
outside of the segment. 
it does lie on the proper 
of test data as desired m 


e used to test the g€ometrical setup 
S of referencing a segment number 
Cified point lies either inside or 
The point is then tested to see if 
‘Side of the Segment. As many sets 
ay be loaded. 
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- imside segment. If JTEST<0O|, 


renee A TTT AEE HEI, CCAS TR A sce 


Description 


Segment number to be tested| 
if JTEST>0, point should lie 


point should lie outside of 
segment. oe 


/ 


Option for testing te see if 

point lies within the confin- 
ing boundaries of the segment, 
or to see if the point lies | 
within this segment rather | 
than some other segment. 


0, Test segment boundaries | 
only . * 
1, Determine which segment con 


tains the point and determine 
if it is the desired segmen 


oe 


=- 


3.2.11 Plot Data 


[ \ Each plot data card indicates a rectangular area within 
\ a plane that slices through the system. This plane may be at 
any desired orientation. The geometric boundaries sliced by 
the plane within the rectangular area are then plotted. The 
only restriction is that the rectangular area must be con- 

i pletely contained within the overall mother region (Region f1). 
This is no real restriction since the overall mother region 
may be made as large as desired without affecting the speed 
or accuracy of the program. 


Description 


) Coordinates of the point to be 
used for the lower left-hand 
corner of the plotting area. 


Coordinates of the point to be 
used for the lower right-hand 
corner of the plotting area. 


Coordinates of a point that lies 
somewhere along the top line of 
the plotting area. 


we a bo = = : - = 7 
A blank card must follow the last of the plot data cards. 


if 


— 
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3.3. OUTPUT | 
Output from the ZAP code consists of. printer plots of 


various cuts. through the geometry and printed summaries of 
the energy deposition for each segment and for each material. 


Under the heading SEGMENT RESULTS the following quan- 


Bi \ 7 
ye MN 
a ' 
© id Fearp rSeae Die oe 


tities are printed for each Soemente: 


Mother Segment number = mother of the segment 


Rs 
being considered. 
2. Segment’ number. 
3. Material number in segment. ad 
4. Total Energy = Sum of energy deposited in 
segment. 
5. Scaled Energy = Sum of energy scaled by lasing 
wavelength. 
6. Energy per unit volume. "i ina eee daieealliaie 


_—w 


-Percent—totat enerpy* The percent of the 


s total energy in the system that is in this 
segment. 
8. Percent error = probable statistical percent 
. error. (See Section 2.5.2) 
9. Number of legs = actual number of ray legs 
traced through segment. 
10. Number of rays = number of different rays | 
traced through Segment, even though each 


‘Sail : ‘ad 


Tay may have had several legs passing 
through the Segment. 
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Under the heading MATERTAL results the following quan- 


are 


Be. 


printed: 

Total absorbed encrgy = sum of segment ener-| 
gies. ,; 

Total source energy - sum of initial ray in- 
tensities. 7 


Total number of rays. 


Final random integer = Random integer that 
could be used to continue the, run without 
Changing the random number sequence. 


each material the following quantities are print - 


1 


Material number 


Material volume = sum of volumes of seg- 
Rents containing material. 


Total energy = sum of energies of segments 
Containing material. 


Scaled energy = sum of scaled segment ener- 
gies for material. 


Energy per unit volume. 


Percent of\total energy = The percent of 
the total energy in the system that is in 
this material. 


Percent error = Probable statistical percent 
error. (See Section 2.5.2) 


Number of rays = number of -different Trays 
that passed through the material. — 


FT ERR AORTA NES bh ARIANA RET fotttim mean 


= 


— 
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3.4 SAMPLE PROBLEM #} 


‘ .'.The first sample problem consists of a cylindrical lamp 
and a cylindrical laser rod encased in an elliptic laser 


cavity (Fig. 3.1). The area of interest is chosen to be in 


the body of the system so that end conditions are ignored. In 
this case each end is considered to be a boundary of Symmetry 


_ (a simple boundary of total ref” ection). 


A single wavelength is considered and parameters such 
as absorption coefficients and refractive indices are not 
intended to be physically real but were chosen arbitrarily 
for demonstration purposes. 


\ 


The laser system is encased in a 20 cm sphere that is 
considered to be the overall Mother segment (segment #1). The 
area inside the cavity was chosen to be segment #2 and is a 
mother containing two daughters.\ These daughters are respec- 
tively, segment £3 which is the cylindrical lamp and segment 
#4 which is the cylindrical laser rod. The laser rod is a 
mother and has daughters #5 and # which are simply the front 
and back half of the laser rod. Note that the area of 
Segment #4 is completely filled by\ segments #5 and #6 so that 
segment #4 has no actual volume of jits own. 


The material outside of the cavity is chosen to be 
material #1, the air inside the cavity is material #2, the 
lamp is composed of material #3 and the laser rod consists of 
material #4. These material numbers| appear in their proper 
places in the printer plots (3.2 through 3.6). 


Boundary f1 is the sphere encasing the system. Boun- 
dary #2 is the elliptic cylinder that | represents the cavity 
wall and is chosen to be a boundary of constant reflection 
with 90% reflection and 10% transmission. Boundaries £3 and 


\a4 ars cylinders defining the lamp and laser rod respectively 


\ ‘ 
\ \ 
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Figure 3.6 Top view of lower half.of laser system used in 
sample problem #1. 
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and allow Fresnel's reflection. The end Planes are planes of 


#7 and is a dumny boundary. 


The lamp (segment #3) is a volume source that emits 
rays with an intensity va 1.0. Only 100 rayS were emitted for 
this sample run. i 
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3.5 SAMPLE PROBLEM #2 


335: Introduction/ : / 


In order to ‘demonstrate the actual use of ZAP, the’ 
Simple laser system shown in Figure 3-7 will be analyzed. In 


‘this laser, a Spherical metal reflector contains a xenon 


flashlamp and a neodymium-doped glass laser rod; the lamp and 
rod are held at their ends by cylindrical holders. When the 
lamp is pulsed, power is produced throughout the volume of 
xenon; the amount. produced depends on wavelength. The power 
travels through the plasma to the wall of the lamp (some is 
absorbed on the way), refracts through the wall, and enters 


the reflector volume. There the power bounces off the walls 


and other Surfaces, travels through the glass laser fod and 
reflector surfaces. We are interested in the fractions of 
the initial power which end up in the rod, lamp, and walls, 
and also in the power distribution in the laser rod. 


3.5.2 Initial Analysis of the Geometry 


We can use the Symmetry of the laser assembly to improve 
the statistics of our results, by putting perfectly reflecting 
planes in the two planes of Symmetry (see Figure 3.7). This 
in effect gives four times as many rays per unit of run time 
(we could get the same result without the planes and their 
time-losing reflection by adding together the depositions in 
Symmetrically equivalent volumes, but the computer time lost 
to the reflections is cheap compared to the human time gained 
in not adding’numbers). The introduction of the Symmetry 
planes reduces the problem to a quadrant of a Sphere, as shown 
in Figure 3.8. |The decomposition of the quadrant into 
successive mother and daughter segments is also shown in 
Figure 3.8. (The rod is divided into eight segments of 
essentially equal volume so that the energy distribution in 
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Figure 3.8 The decomposition of the 
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3.5.3.1: Reflection Data 
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° 


6 


it may be found.) The decomposition proceeds in a’ straight- 
forward, way from the: all-enclosing segment 1 through smaller 
and smaller segments, always keeping the nuaber of daughters 
small at each level. In fact, the number of. daughters per. 
level may be unnecessarily small - a comparison run with more 
daughters per mother in the red gave essentially identical, 
run times. Once the mother-daughter decomposition of the 
geometry is decided, the preparation.of the input is easy 
(if se ea : , 


\ 


\ aa) 
3.5.3 Input Preparation \: 
\ . 
While reading the following explanation of the input 
information the reser should refer to the listing of the 
computer input and output reproduced.at the end of this 


section. | : yy 


! 


. Aside from the buiit- -in reflection types 1, 2 and 3 
(dummy boundary, total reflection, and dielectric boundary ' 
respectively) we will need | ‘three defined reflection types: 


1. reflection waryite smoothly from -85 for 
rays normal to the surface to 1.0 at 
grazing incidence to represent the reflec- 


tor and other internal metal surfaces; 


2. diffuse reflection returning half the 
incident li ght, 
eed et and electrode structure . 


to, roughly represent the 


in ek lamp 


3. reflecticn Vineaie smoothly from .95 for 
rays normal to the surface to 1.0 at 
grazing incidence to model the dielectric | i 
mirrors at the ends of the laser rod. 


ff 
ard 
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a These three types of reflection are entered as reflec- 
tion types 11, 12 and 13. aa 


3.5.3.2 Material Data S63 


The metal, air and quartz materials are entered in a 
straighiforward manner (the refractive indices for the metals 
and air are made slightly different so that ZAP? will know 
that the boundaries between them are not just dummy boundaries, 
as it would assume if their indices were the same - see sub+ 
routine BOUNC «:) : TRe xenon lamp plasma and Na>* - doped’ glass 


have absorption coefficients which vary with wavelength, and | 


the user-written subroutines ABS2 and ABS3 must be supplied. 
These routines (see listing at end of section) are self- 
explanatory. . , 


The abgorption coefficient of the metals is made very 
large so that any light entering them will be immediately 
cabsorbed. Note that a small ‘absorption has been added to the 
glass: at all wavelengths; this is to insure some absorption 
(even if small) for al] wavelengths, to avoid rays that 
dribble on forever. ‘ 


5 a es ‘ Boundary Data j Me oa! 


‘ Fie ge 8 on is rae 5) - 
The boundaries fin this case only sphérés, cylinders r 


and planes} are: just those nécessary ‘to Specify the, sepments - 


| previously decided*on. Because ZAP currently handles surfaces. 
through. the. origin, the cehter of the geometry is displaced 
“by .001 on’each axis. When two reflection properties are ime 


required for different: parts of a boundary, two coincidir.; oe 
boundaries are specified with’ different reflections, and the 
choice between’ them is made when the bouhdaries.of ‘segments 


are specified (note boundaries 5, il and 13) 7 


‘ were pee are! 
: , : : 6: 
» { a. ens 
ral . 
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3.5.3.4 Segment Data | : ' 


\ 


Once the boundaries are specified, the segments are 
defined by giving all boundaries of each (and the volume; if 
the energy density is desired for that segment). The process 
is quite ernEr es 7 


3.5.3.5 Structure Data 


The mother-daughter relationships are entered for each 
segment. The information is most easily prepared by aprerence 


to a decomposition diagram ‘such as Figure 3.8. ¢ 


3.5.3.6 Lamp Data 


For this problen there is only one type of source. 
This is the volume- -emitting xenon plasma in the lamp. Its 


properties are specified in the user-supplied. ‘Subyfputine 


INTEN1, which is reproduced at the end of this section. 


e ‘ \ 


a 


3.5.3.7 Source Data 


There is only one source, namely the xenon Tsim. It. 
is specified geometrically by one boundary and a rectangular 
paralielopiped. For this run 500 rays are called for from 
the lamp, in order to get acceptably low errors in: the segment 


_ deposition energies in the rod. 


= 


r / 7 : 
? gi. gs 3" 5 { 
/ 
Test points are specified for all segments: in this 


3.5.3.8 Test Data - 


example, although ae is not Necessary in general: Option’ 1 


; as used. except when the segment volupe is zero (because ~ 


G%, 


daughter segments am up ali its volume) ,-in «which case | 
Option 0 is used. egments 1 and 2 also hare points speci- 
fied outside ther (negative segment ‘Rhumber). -All points are. 


SQEEECES so no error messages appear. 
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3.5.3.9 Plot Data \ 


\ 
Five plots are called for. There is a cross section 


near the center of the sphere (but just above the symmetry 
boundary) perpendicular to the rod and lamp; detail of the 
rod and lamp is poor because of the limited printer resolution. 
Another plot parallel to and through the rod and lamp suffers 
from the same lack of detail, although the gencral configura- 
tion is clear from both. The third plot is a close-up of the 
lamp-rod-holders-sphere area oi the secord plot; it clearly 
shows the quartz lamp wall, rod segmentation, and holdér 
details. The fourth and fifth plots are close-up views of 
cuts across the rod and lamp respéctively. These plots are 
very useful in checking that the geometry has been specified 
just as desired, and that no silly errors have crept in. 


3.5.4 Discussion of Output 


The output due to test points and plot requests has 
been discussed above. Here the power deposition results will 
be discussed. 


3.5.4.1 Segment Results 


The most interesting information in the segment results 
is the pattern of power deposition in the laser rod. The 
half of the rod near the sphere wall (segments 14, 15, 16 and 


(17) and the half near the center (segments 20. 21, 22 and 23) 


get essentially the same amount of scaled power -. 47.13% and 
52.86% respectively (we get these figures by adding the’ nua- 
bers in the percent scaled power column). The lensing effect 
of the laser rod has concentrated more power in the core of _ 
the rod (segments 14, 16, 20 and22) than in the skin (seg- 


- ments 15, 17, 21 and 23); the core got 60.46% and the skin 
~ 39.53%. This effect is apparent in both the center and end 
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Ss 
sections of the rod. Because the 500 rays have been enough 

to give approximately 10% errors in the deposition powers in 
the rod segments (see percent error core) these conclusions 
are reasonably reliable. : 


3-5.4.2 ‘Material Results 


, The material results tell us the distribution of power 
among the various elements of the laser assembly. : From the 
percent total power colum, we see that the lamp absorbs by 
far the greatest fraction of the initial power. However, the 
lamp is a volume emitter and much of the power created in its 
volume may be absorbed before it ever gets out. To find out 
how much is internally reabsorbed in this way, another run 


| Was made with the absorption coefficient of the "air" in the 


cavity equal to 1000 cm rae The lamp was found to absorb 
76.59% (error 2.29%) of the radiation before it got out for 
the first time, leaving 23.21% (erfor 7.64%) of the initial 
power as the amount actually enitted from the lamp (the 
remaining .2% was absorbed in the lamp ends). Thus we find 
that the cavity and holders absorbed 9. 58% of the initial 
power (41.28% of the ‘emitted power), the rod absorbed 6. 70% 
(29. 30%), and the lamp reabsorbed 7.02% (30.25%). . The power 
returned to the lamp will act much like the lamp electrical 
input power (because the thermalization tize in the dense 
lamp Plasma is very short). If a fraction R of | the lamp input 
is radiated as light with the remaining 1-R lost as heat, and 


—if_R is. not changed. ‘greatly-by the-self-loading. due to returned 


power (probably a good assumption for the low self-loading 
here, but more questionable with higher loadings), we can make 


a rough self-consistent calculation of what happends to the 
returned power and what the deposition fractions are as 
; \ : 


—- 
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follows ‘assume all spectral regions act like the .4 - .92 » 


region we calculated): 


The input electrical power is P 


| The power returned to the lamp from the cavity 


, is Q . 
. 
| The total Oe thus P+Q 
| The lamp light out is R(P+Q) 


| Then Q = R(P+Q)F or Q = (RFP)/(1-RF) where F 
is the fraction of lamp output reabsorbed 
| | (F is: .3025 for the geometry we analyzed). 


Then if F! js the fraction of P re-absorbed, 
| | we have F'=Q/P = (FR)/(1-RF) 


And if L is the power absorbed by the laser rod. 


@ \ And if G is the fraction of lamp output which 
the rod gets (G is .2390 for our laser) 


| Then L = GR(P+Q) = (GRP)/(1-RF) and so 

G) = L/P = (GR)/(1-RF). Similarly, if C is 
| the fraction of lamp output the cavity gets 
| ; then C) = (CR)/{1-RF). / - | 


Thus we see that self-consistency requires that all deposition 
~~ fractions be multiplied by the factor R/(1-RF). - Taking the 
| reasonable value of .5 for R, ‘we find this factor to be about 
-59 for geometry we analyzed. ee ' 


° pereen, 


From the material power and scaled-power for the lasez 
glass (material 6) we find that 57% of the power the laser 
[ . glass absorbed was useful laser power (the remainder was con- 
verted to heat in the rod when the absorbed photon energy was, 
reduced to upper-laser level energy by radiationiess transi- 
tions after absorption). The self-consistent fraction of 
electrical input the rod absorbed is -18,, so .10 was useful 


| — oF i” 


ih: io ae 
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for laser operation and .08 went to rod heat (we have assumed 


\@ quantum efficiency of 1.0, which may be too high). 


| . . 
' | Our results, modified by the above self-consistency 
a ument , thus imply that: 


: The electrical input was 1.00 
_ The lamp heat loss was 58 
7 The cavity wall absorption was “+ 224 
‘The laser rod heat loss: was 08 
The uppér-laser-level result was -10 


3.5.4.3 Summary Results’ 


The total absorbed power and total source power are 
slightly different because of statistical fluctuation in the 
way rays are killed off. The execution time (on a CDC 3890, 
which has floating point add-substract, multiply, and divide 
times of about 2, 5 and 10 us) was, 1855:469 seconds. A run 


‘with one ray shows that the setup time is about 30 seconds, 
and thus the run time was roughly 


& T * 30 + 3.7 N seconds 


where N is the number of rays (each with 200 wavelengths) 
traced. The program and storage occupied about £0,000 words 
(48 bits each) of storage. i 


a 

3.5.5 Conclusidns 
We have found that about 103 of the input electrical - 

power is converted into a pumping rate into the upper’ laser 
level for the laser geometry and lamp conditions we specified. 
The actual laser output available will Be decreased by spon- | 
taneous emission and (for /pulsed lasers) by the fact that not 
all the upper-level energy can be extracted. However, the 
upper-leve1 pump rate i's robably the best measure of cavity 
efficiency and thus the best number to use when comparing 


cavities. im 
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SAMPLE PROBLEM #2 - USER SUPPLIED SUBROUT INES 
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SAMPLE PROBLEM #2 - USER SUPPLIED SUBROUTINES 
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SAMPLE PROBLEM #2 - TEST RUN INPUT 
EE RUN INPUT 


MULTIFREQUENCY TEST ON SPHENE WITH WOU aND Lane 
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APPENDIX A 


A.1l HELIX 
A.1.1 Introduction © 


This Appendix describes the procedure that is followed 
for finding the intersections of a ray with a helical pipe. 
Two related procedures for determining whether a given point 
is inside or outside of a given helical pipe and for determin- 
ing the normal to a helical pipe at a given point on its sur- 


face are also discussed. 


A.1.1.1 Notation and Conventions 


The region R of interest in which the helical pipe 
lies has a fixed coordinate system. An equation like 
¥ = [x,y,z] means that ¥ is a vector* with components x,y, 
and z. Vectors will be used to represent points, line segments, 
displacements of points and directions within the region R. 
The length of VW is denoted by |¥}. |[x,y,z]] = (x2+y2422)8, 
Curves will be represented in parametric form. For example, 
if P = [x,y,z] and D = [a,&,y] is a direction (in which case 
g24g? xy" = 1), then £(s) = Btsd = [x,y,z] + s[a,8,y] = 
[x+as, y+8s, zt+tys] is the line through P in the direction of 
D. More precisely, ,it is the point on the line that is a 
distance s from the point P. Similarly Hi(a) = [x+pcos6, 
ytesin’, z+«x8} is a helix with axis parallel to the z-axis 
and passing through the point [x,y,z]. (More properties of 
the helix will be given in the following sections, particularly, 


the next.) 


* 
The notation used is essentially that of Louis Brand, Vector 


Analysis, Wiley & Sons (1957) except that he uses hold face 
type instead of the arrow to denote vectors. 


Al 
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Surfaces will be represented in two Ways: by parametric 
formulas and by locus equations. For example, the plane pass- 
ing through the three points Pas Pas and Ps 1S given by the 


parametric formula 


Z(u,v) = By + u(¥,-5)) + v(8,-5,) (A.1) 


or the equation* 


o- (B-B,) - (F,-8,) = (8,-8,) (A.2) 


which is satisfied by just those points that lie on the plane 


and no others. 


A.1.2 The Helix 


A parametric fornula for a helix was given in A.1.1.1. 
With the addition of a subscript h on sore of the terrs, it 
is the one that will be used consistently in all that follows, 
viz. 


> 


H(§) = [sn + ecos?,y, + psins,=, + a | . (A.3) 


The axis of the helix is parallel to the =-axis and passes 
through the point Py = [xp >¥,o=pJ- The radius of the circular 
cylinder on which the helix is wrapped is <. The pitch of the 
helix is 27x, the distance traveled in the = direction per 
revolution. The helix is right-handed if «>0. The paraneter 
€ represents the angle of rotation required to move fron the 
point H(0) = [x + CsYpo=y] to the point H(z). 


——— 


* 
See Brand, §10 and 11 for definitions of scalar (dot} and 
vector (cross) products. 


3SR-405 


A.1.2.1 Tangent 


The tangent* t(3) to the helix at the point fic) is 
found by differention: 


t(9) = _ Hi(e) 


= [Se (*p + pcos) , a(n + osins), gy (2, + xe) (A.4) 
[-psin @, pcos @, Kk]. 


The length of t(3) is the distance the point H(8) moves along 
the helix ae 3) inex Eases by 1. The unit tangent 

T(¢) = t(8)/(o" + K 2's is the unit vector in the direction of 
the tanger:. It is the direction of the line tangent to the 
helix at HAG ), so the parametric form for the tangent line is 
Hc) + sT(é) where S, not @, is the parameter. 


A.1.2.2 Principal Normal 


The principal normal** to a curve at any point is the 
direction from that point to the center of curvature. It is 
perpendicular to the tangent. It may be found by differen- 
tiating the tangent once again. Since d?(8)/dé = [-pcos 8, 
-psin 8, 0], which has length ¢, the principal normal is 
N(2) = [-cos %,-sin 9, 0]. It is always perpendicular to the 


Z-axis as well as to t (3) and it points from H(9) toward the 
axis of the helix. 


---————————— 


a 
See Brand, §27. 


x 
Brand, §29. 
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The radius of curvature of the helix is given by 


antsy 


which can be worked out by straightforward algebra and differ- 
entiation. 


The binormal of a curve is the direction perpendicular 
to the tangent and the principal normal which may be deter- 
mined by taking the vector product of those directions. In 
our case, the binormal is given by 


M(3) = Fe) « Xs) 
7 [ oe ee 
= |Ksin 6, -«xcos #0 | Ue +K ) ¢ (A.6) 


A.1.3 The Helical Pipe 

The sur‘ace that we must deal with is, loosely 
speaking, a circular cylinder whose axis is a helix which is 
to say that a plane perpendicular to the helix at some point 
Bi cuts the Ssuerace, 20 the immediate vicinity of Py in a 
Circle with center Py and radius a. A parametric formula for 
the surface can be given as 


> > Eee is tes ps 
Q(é,9) = H(6) + a cos 4 N(é)+a sin ¢ HW¢(s) 
= [*n + p cos = - a cos ¢ cos § + «a sin ¢ sin =/ 
2 2,% om 
(e +K ) » ¥, + P Sin 3 - a cos ¢ sin 2? - ka sin ¢ cos $/ 


(02x? 2 2. 2 Pl 
p +x ) » 2 + K§ + pa sin e/(e +x ) : 


Ad 
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nOtice that when c=0 the formula reduces to one for a circular 


cylinder of radius 2 and if x=0 it reduces to one.for a torus. 


A.1.4 The Ray and the Helical Pipe 

The origin of the ray will always be denoted by 
e. cs [Xp >¥9>=9] and its direction by D = {[o,£,y] where, of 
course, a*+37++? = 1. We will always use P(s) = By+sb 
=[xp*s2,¥9tsi,zptsy] as the parametric forzsla for the ray 
where s denotes the distance from P, to B(s). Since By is 


the origin of the ray, s>0 for points of the ray. 


Let b denote the distance between the axis of the helix 
and line of the ray and let P(s,) (it is possible that S,,<9) 
be the point which is closest to the axis. Since B is per- 
pencicular to the shortest line joining P(s,) to the axis, we 
must have a (X},- Xq-S)2) + 304, -¥9 75,5) = 0. Therefore 


=, = [a (x,-x9) * 3(¥y-¥9) | o7+2?) (A.7) 
and . 
4° = (%,-%9-spa)” + (n-%9- Sp)? A (A. 8) 


As we soon will see, it is important to know whether 
the ray passes the axis- on the Tight or left. The ray passes 
the axis on the right if and only if the three vectors ee 
Dand < forn a basis for a right-handed coordinate systen 
which, in turn, is true if and only if the triple product 
a = (Por P)* @ x k) is positive. But 


Dx K = (2,2,y] « [0,0,1] = [2,-2,0] ; (A.9) 


¢ 
rq 
3 
f 
z 
: 
4 
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and hence 


g = [*o-*n-¥on’207Zh |” [é,-2,0] (A. 10) 


= 8 (x9-x, } - a(¥o-y},! 


The magnitude, as well as the sign, of ¢ is significant. 
Notice that % x Kk, being perpendicular to both D and k, 1s 
parallel to the shortest segment joining B. + sk (the axis) 
to B +s (the ray}. Hence, b must be the length of the pro- 


0 > > > 
jection of PoP on D x k, 
> _p c ca] 
P 1(Po P,) Bx k! , 
1B x E| (2.etls 


A.1.5 The Turns of the Helix 


When the ray and the axis of the helix are skew, there 
is a unique plane ry containing the ray which is parallel to 
the axis. Consider, then, the half-plane extending away fron 
the ray in a direction perpendicular to iL and having the axis 
as boundary. This half‘plane cits the helix into a series of 
individual loops, that we will call turns, whose endpoints are 


~~ 


as far as possible from the plane ae The turns are arcs of 


the helix Hi(s) corresponding to intervals +. « : < ae | where 


nt” = 
ny = Vot2n-, and \ 
sin vy ~ (y_7¥ 97 25,170 | (A.12) 
COS vp = (Xy7Xq735,,) 7% (Aco } 


Ao 


i 
/ 


i 


3SR-405 


The vector [cos Yo? sin ¥9 29] is parallel to the short- 
est segment joining the axis to the ray, which is perpendicular 
> 
to D. Therefore 


acos Yo + Bsin Yo = 0 (A.14) 


A.1.6 The Sequence of the Turns 


The basis cf the plan that is used for finding the 
first intersection is a procedure for calculating the inter- 
sections of a line with one turn of the helical pipe. For 
the sake of efficiency, it is desirable to know the order in 
which the turns are encountered by a point as it travels 
along the ray. The present section is devoted to the descrip- 
tion of a method for plotting rays that is well-designed to 
illustrate how they proceed from turn to turn. 


To see how the plot may be gencrated, consider a half- 
plane that is boundea on its edge by the axis of the helix 
and which cuts the ray at some point and let a curve be 
traced out on the half-plane by the cut-point as the half- 
plane rotates about its edge and the point travels along the 
ray. 


Let r(s) be the distance from the axis to P(s). Then, 
2 2\% 
Ee Soe Te een) ee 
4, 1 
2 2 2 
= (b on (s-s,) 


Let w(s) be the distance from P(s) to the half-plane that 
passes through P(s)). 


A7 
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Then 
w(s) = u(s-sb) (A.16) 
where 


2 52)" (A.17) 


i.e., wp is the length of the projection of D on the plane 
z=0 and it is positive if the ray passes the axis on the 
Tight as it does in Fig. A.2. It follows fron Eqs. (A.10) 
and (A.11) that 


vw =o/b . (A.18) 


Finally, let 
¢(s) = arctan w(s)/b (A.19) 
and . | ; a 
z*(s) = Ze ays = Me): (A. 20) 


The significance of ¢(s) is shown in Fig. A.2. It is the 

angle that the plane through P(s) makes vith the plane 

through P(s,). The sign of » in Eq. (A.17) was selected-so 
that ¢ and 9 are both measured in the right-handed or counter- 
Cloc!}.wise sense. The purpose of -(s) is to prbvide the terr 
 «O(s) in Eq. (A.20) which compensates for the tise (or fall) 

of the helix as the point on the Tay passes th axis. Thus 
z*(s) represents the height of P(s) above a that rises 
with the Helix. ~The desired plot, then, is generated by a 
point having r(s) as ordinate and z*(s) as 
plot, the turns 6f + the helix appear as fiked points with 


+: vr -& - =a -k ey ee ° la) Pa 
coordinat (rT, >= a? where Tht and = 2 = * <(~9* (2n+1) as 


AS 


Finire, A.1l <A view of the helix 


of PRE helix. 


A9 


cutting plane 


direction 
of X- axis. 


Y 
Z 
A 


and ray looking down the axis 


fh aa 
ny 


1 
rete be ee ne 


a 
' 
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Figure 


A.2 


Another top view 


of the helix and ray. 
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The equation of the trace of the ray is obtained by 
eliminating s from Eqs. (A.15) and (A.20). When «x=0 the 
result is 


2 | 
(a 2 (5) a (A.21) 


where’ 2), = Zot 5S,» SO that the trace is a hyperbola whose 
center has coordinates z* = z, and r= 0. The effect of k # 0 
depends on the relative values of xp/y and b. If xp/y<0 then 
the-ray.-goes counter to the helix and the hyperbola is . 
stretched by the amount mk, the distance the helix rises dur- 
ing a half-turn. In this case, the ray meets the turns of 

the helical pipe in their natural sequence. 


If xp/y>0 then the ray and the helix both rise or both 
fall and if the magnitude of « is large enough so that xp/y>b 
then helix rises or falls faster than the ray while the ray is 
close to the axis; the hyperbola is compressed by ™« and a , 
loop forms. , 


To verify this fact, observe that 


dz*(s)/ds = ¥ - 45 
I (s-s,) + b 
oF (A.22) 
Ky 
= ¥ 1 - yb. 


fu(s-s,)] 54 
(BS 


which changes sign if <u/y>b at s=s)*2 where 


| i fate : b) | (A.23) 


All 
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y, 


When O<xy/y<b, the hyperbola is still compressed by ~« but 
not enough to form a loop. , 


‘e 


Figure A.3 shows examples of the three types of 
traces that rays may have. The traces of the turns of the 
helical pipe are more or less elliptical. It is seen that 
when «p/y<b, the ray intersects the turns consecutively, each 
one at most twice. When b<<y/y, all turns except the two 
nearest z*=2) are still passed in natural sequence and each 
is intersected at most trice. On the othe> hand, the two 
exceptional turns (i.e., the two with z* nearest Zp)» may 
both be intersected four times by a ray whose trace resembles 


the one pictured in Figure A.3. 


A.1.7 Numbering the Turns | zi 


The strategy that is Ised to find the first intersec- 
tion of a ray with a helical ipe is to determine the sequence 
in which the ray approaches successive turns and then to 
examine the turns in that sequence to find the first that is 
actually intersected. In implementing this strategy, the 
turns are assigned ordinal numbers as follows. 


Let 


v= 2=sgn(xy) , %. = %9 * nv 


where 


n 
° 
4) 
« 
o 
I 


= (¥n-%07F5p)/P 
and | 


sin “9 = (Yn-So7oSp)/P 


Then, turn n is the turn for which 7_<i<v . 
n n+l] 
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Figure A.3 Traces of typica? rays: 1) Ku/y<0, 2) O<xi:/y<b,' 
3) b<cy/y, 4) ray with 8 intersections 
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This definition of the number of a turn differs from 
that of Al.5 in that here we require that the sign of vy 
agree with tha: of xy. This has the pleasant effect of insur- 
ing that n will increase as s increases (except for that 
segment of the ray corresponding to a loop in the trace). In 
Order to give precision to the above Statement, the function 
n(s) is defined to be the ‘number of the loop that is closest 
to the point B(s) = By + sb. The value of n which is assigned 
to points midway between turns is unimportant since such 
points are necessarily exterior to the helical pipe. 


A.1.8 Seven Special Points on a Ray 


There are certain points on a ray that are special 
enough to receive specific names. They are listed below 
together with a formula for the value s at which the point 
is attained. When the value of s for some special point on 
a particular ray is negative or complex, it means that that 
point does not occur. The points are 


Po» the Origin of the ray. Sg = 0 

Pi; ‘the point where the ray enters the cylinder 
‘that encloses the pipe. S$) = Sp - W(p+a)?-b?. 
| If b>p+a, the ray does not penetrate the 
cylinder. 


| 
P,> the poirt where the ray enters the cylinder 
| on which the pipe is wrapped. 


Sp = 5) - WY (p-aj7*b?. If b>p-a, the ray 
does not penetrate the inner cylinder. 


Ps, the point where the ray emerges from the 


aN — inner cylinder. S. = 6) V (p-a)?+b2. 
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Py, the point where the ray emerges fron the 
outer cylinder. Sq = S, + V(pta)?+b2. 


Ps» the point where the ray reaches the bottom 
of the loop in its trace. Sc = $,*4- There 
is no loop if \ is complex. 


6? the point where the ray reaches the top of 
the loop in its trace. S5 = S,7A- 


The number of the turn nearest Py is my = n(s,). Note 
that 


"2 = 52 5 53.5 54 4 S5.> 5, 
mE = "2 503 Sm» ny <n, 


whenever the indicated quantities are real. 


The subroutine INTHEL which we consider next is 
not exactly simple but it is straightforward. The complica- 
tions that arise have to do with the Singularities that arise 
when one or more of the quantities y, yp and b become too 
small. It will be observed that when Ng?n, these two turns 
are treated Simultaneously because they can generate inter- 
locking intersections. | 


A.1.9 The Distance fron a Point to a Line 
==" ._fO a Line 


The distance fron the point P = [x,y,z] to the line 
{(s) = Pytsb = [x+tas, y+3s, z+ys] is attained where D-(£(s)-3) 
= 0, i.e., the segrent that joing P to the Closest point on 
{ is perpendicular to the line. If S=s,, at that point, then 


B- (6, -B) 
mB 


= a(x9-x) + E(yo-y) + ¥ (29-2) 


Ss 
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and the distance from the point to the line is 
d = (x + as -x)? + (y + Bs -y )? + (z + ys -z)? : (A.25) 
0 m Yo m 0 m site 


A.1.10 The Distance from a Point to a Turn 
Sete trem a Point to a Turn 


It is often necessary to determine the distance fron 
@ given point to a given turn of the helix. The distance 
> 1 
from P = [x,y,z] to Hi(e@) is easily verified to be (2G)? where 


2 2 2 
G = 5 (x-x,, - cos @) + +(y-y,- sin 6) + 4: (z-z, -xe) ; (A. 26) 
A change of @ will reduce the distance if it will reduce G. 


The distance will be a minimum when dG/d9 = 0 and d2¢/d6259. 
The first two derivatives of B turn out to be 


a 


do ~ esin 8 (x-x, -pcos )) - pcos 0(y-y,-psin 9) 2 x (z-z, -«6) 


= psin @ (x-s,) - pcos &(y-y,)- K (2-2),-«6) (A.27) 


a. pcos 6(x-x,) + psin 8(y-y,,) +x? 


‘The problem reduces to that ef finding the nearest root of 
dG/de. A given estimate of the 6 of the closest point is 
improved using Newton's method «hich is to set 
dG 
7 
) = @--— . (A. 28) 
new dg 


de” 


It is well-known that Newton's method can fail and the present 
application is no exception. The trouble arises when 


a’G/de2<9 which usually causes a maximum rather than a ninimun 
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to be approached. The Problem is overcome by setting 
eee = 0-sgn(dG/d6) whenever d7G/do2<0, i.e., by moVing one 
radian toward the minimum. 

The question of when to stop this iterative approach 
to the minimum cannot be answered out of context and there- 
fore we leave it to be discussed when it occurs in a specific 


application. 


A.1.11 The Inside of the Pipe 
0 the Pipe 


It is sometimes necessary to know if a given point 
P = [x,y,z] is inside or outside the pipe. The point B is 
inside or outside the Pipe. The point - is inside if its 
distance fror the helix is less than a. This assertion 
follows from the definition of the surface of the pipe which, 
according to the parametric formula given in A.1.3, consists 
of those points lying at a distance a from the helix. Thus 
the problem reduces to finding the distance from P to the 
helix. 


| 8, = ATAN2(y-y,, x-x,) + 2mm (A. 29) 


The value of n which minimizes [P-H(0 | is the one that 
minimizes [z-z,-«a| which is, in turn 


zZ-Z, -Kd 
ae [3 + | (A.30) 


where [x] denotes the greatest integer that is not greater 
than x. Caution: in FORTRAN, INT(x) = [x]+1 when x<0. 
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Using H(6 ) aS a starting point, a procedure like that 
of Al.10 is used tec search for the closest point with the 
refinement of a quadratic term. Thus 


a°c = -psin a(x-x ) + pcos @(y-y (A.31) 
and we write 
(°nen"®) 
2 6. 76 . 
dG d°G new d“G 
O= + (6 -6 + —— -2 (A.32) 
6 ( new ee 2 de> 


Writing G', G", G"' for the Successive derivatives of G and 


for O ew 08 we get 


é 7 -G" + (coy? - 2G'g": )% 


: (A.33) 
2G! t 


=o" , 
Magan? |. Saab othe 
((6") 7 26g") 
The: second form, which is obtained from the first by ration- 


alizing the numerator, is preferred, especially when G'" is. 
small. . 


The iteration may well be Stopped when [S|<.01 because 
the error will be less than a(I-cos&) < 5x10 74a. 


The method tails for points on the axis but these 
points are outside the Pipe because ax<p. 


‘ 


/ 
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A.1.12 The Distance fron the Ray to a Turn 


A ray approaches every turn at least once and never 
more. than twice. If it approaches turn n at two peints, e. 
and eb then Os and ey will be separated by Yn + &v, the mid- 


point of the turn 


The procedure for finding where and how often a par- 


, ticular turn approaches a Yay is a combination of. techniques 
described dbove. First, the point 6, = vite is selected. 
The equation (A.24) is used to find the point on the rav 
Closest to H(s). Let its parameter be called s(0,). Next the 
method discussed in Al.10 is applied to find a 8, such that 
H(e,) is closer to P{s(9,)) than H(8,) and then s (65) is 
calculated. The’ process is repeated until a convergence 
criterion is satisfied. When the procedure was first designed, 
the plan was to apply just one step of. the algorithm described 
in Al.10 per iteration but certain test cases showed Signs 
of instability. The dynamics of the instability are not 
understood but it was found that repeating the algorithm 
until [6 ew Ol<-01 was sufficient to damp the oscillations. 
The whole procedure is repeated with 8 =u + ay. 
When the ray approaches the turn on both Sides of vty. the 
points of closest appfoach are found with an accuracy of 
better than .01 in S. I¢ the ray approaches the turn on only 
one side, the procedure on both sides will usually converge 
to that point but occasionally » when $, is set to the "wrong" 
Side, it will drift around the back of the helix to an adjacent 
turn. 
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A.1.13 Finding the Intersections 

If a ray is found to be within a distance a of a turn at 
a point where it is closest, it can be inferred that it inter- 
sects that turn. Tro procedures are used for localizing 
intersections, one is used if 1D-T(2) {<0.3, else the other is 


used. 


In case [D-F(s) [<0.3, a first approxination to the 
location of the intersections is made by calculating the inter- 
section of the/ ray With a circular cylinder of radius a and 


axis through i¥(8) in the direction ¥(s). If c is the distance 


from ii(e) to the ray then the intersection of the ray with 
the cylinder will be at £(8)+4s where 


22 <, 
As = aaa (A.34) 


The first approximations: to the intersections of the 
ray with the turn of Pipe are further refined. | Let Ba 
= [x,5¥,°Z,] denote one of the first approximants. The pro- 
cedure begins with. the application of A1.10 to determine the 
6 corresponding to the point on the helix that is Closest to 
the point Po: Then the closer intersection of; the ray to the 
tangent cylinder =t the new 9 is determined: The equations 
used are discussed below. First, the equation of the cylinder 
is 


[P-fice)|* = a%+ B-ii(sy) + Fee) 2 (A.35) 
Next, we substitute Btu, for P in Eq. (A-55) and solve for 
u, which yields i 
2 =a 2 > 4 2 = 
u°(1-(B-T) ) + 2ub-(B,-i) + jP_-Hj* - a? = 9 oo (A.36) 
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a~ & 


*o b(1 + (1+a(a?-<2)/s?)*) se 


HFo-Al’ , d-1- @-f)?, b= 5-(8,-A) 


where 


and 7- (Bo -fi) = 0 (because the segment joining By to the helix 
at A is a ee to the tangent T). 


\ Then, ‘Stu is substituted for s and 6 is further refined. 
If \the change in 6 is less than .01, the current value of s 
is accepted, else another step. is taken. 


When |D-T|>.7, the approximation in Eq. (A.35) will 
often yield such poor accuracy that the refinement technique ae 
fails. For that reason, an entirely different method is used 
when [5-T] is large. The criterion fOr ‘the alternate method 


is |B- Fl>.3. ‘ ; 


’ 


In the beginning of this section, s (8) was defined to 
be the vaiue of the parameter of the point on the ray closest 
to fi(e). We now define a new point on the ray to be associated 
with @; and that is the point Where the ray intersects the 
plane (8) that is normal to (6) and is’ cut by the helix at 
H(9). The equation for (9) 


P-F(8) = Hi(e) -F(e) (A.38) 
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substituting BytsD for BP in Eq. (A.38) yields an equation in 


s, the parameter of the desired point: | 


(Ears) (-osin 9) + (ygrts)ocos 0 + (2q¢v5)« 


xX, +pcos 6) (-psin 3) + (y,+psin : oe 9+ (z+ xE)}x 
h ) 0 h 


| 
Denoting the solution by Ss, (8), we get 


s, (8) = 


-apsin 6 + Bpcos 8 + YK 


ps that the denoninator of es “RHS of Eq. (A.39) is 


(Xy-¥o) Cosin 8) + (y,-yp)esin © + (2,-z9%«9 )x 


(A.39) 


(p? +K a a5. “Ff, making sp (8) large when b- T is small. The impcr- 


tant thing about the point B(s (8) i il that (2) is locally 
the closest point on the helix. . Ther fore, P(s (8)) is a 
point of intersection if the distance [PCs (9)) -Ai(e) | = 


Let B(®@) be defined as 


B(8) 4 (xptas, (9) - x, ~ pcos e)? 


+ 


's(yg*8s,, (2) - ¥, > sin e)? 


* h(zotys(6) - 4, - <6)? 


(A.40) 


or, in vector notation, B(S) = &1B(s (¢))- H(6){? or, in words, 


hhalf the square of the distance fron B(s (3)) to fic 2). 


the problem of finding intersections reduces to that of soly- 


ing the equation 


(A.41) * 
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When [D-F (98) | <<a, Eq. (A.41) is hard to solve because B(8) 
becomes very large in the vicinity of the root and consequently 
the cylinder approximation is employed in such cases. 


A.1.14 Findiag the Roots of B(3) = a2/2 


The process of finding the roots of Eq. (A.41) involves 
dB/d@ which we denote by B', as well as B itself. The formula 
B' follows from Eq. (A. ); 


—~. B'(0) = (xptas_ (8) - X- pets a) (as. (8) + psin e) 


+ (yo*8s,, (8) ~ Y,7osin a)(3s * (8) - pcos ) \. (A.42) 


+ (Zo*vs, (@) - z,/- x9) (ys*(8)-x) 


where 


s(0) = (px'-xp') sp? | (A. 43) 


N = (x,-x9) (-osin 6) + (y,-Y9) pcos ts) . (A.44) | : 


¢ 


ce it i + (2,729 **8) x 


| 


[Nt = (%,-%g)(-ecos 9) + (¥n-¥9) (-esin 6) + «@ (A.45) : 
D= “apsin 5 + Spcos 8 + yx (A. 46) 
D' = -agcos 9 - Boesin 8 | (A. 47) 


The procedure is to examine B(5) and B'(2) at Steps of 
about 1/4 looking for intervals that might contain one or more 
roots. When a change in the Sign of B(8)-a“/2 is observed 
between the endpoints of an interval, a rather Ordinary root aa \ 


} 
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IN 
} 


finder is evoked to find the location of the locatioh of the 


intersection to the desired accuracy. But a sign Change in 
B(6)-a7/2 will occur only when a Single root lies within the 
interval. Two indicators are used to detect the occurrence 
of double roots. The first is a sign change in B'. If B! 
changes sign as an interval is traversed, then an extrenun in 
B must occur within the interval. The extremum is ignored if 
it is in the wrong direction, e-g-, if it is a maximum and B 
is greater than a*/2 at the end-points. But if the extrenum 
is in the right direction, a linear interpolation of B' is 
applied to provide an estimate of.@ at the extrene point. If 
B is on the other side of a*/2 at the estimated extreme point, 
then two roots are bracketed and the root-finder is called 
forth for each. Otherwise, the reduction of [B-a7/2] is con- 
puted and if it is reduced by more than 1/2, the interpolation 
is repeated, if not, the interval is abandoned. | 


The second indicator of a double root is used to detect 


Pree ane ATR ty) 


intervals in which a double Sign change of B' has occtirred. 


_ B(9,) - B(e,) 
(91-82) (B°(9,) : B6,)) 
where 8, and 8, are the values of 6 at the end points of the 


interval. Thus, c is the ratio of the difference quotient 
of B to twice the average of the derivatives at the endpoints. 


Cc 


For a quadratic, c = 1/2. If c<1/4, the interval is considered 


to be too large and it is subdivided. The subintervals are 
then processed in the Same way as full-sized intervals. 


In the above, it has been assumed that no interval can 
contain more than two roots of B-a*/2. To be sure that such 
is the case, each turn is divided in the middle into two half- 
turns which are treated separately. Each half-turn contzins 
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at most two roots and at most one Minimum of B'‘. 


fact is also used to reduce the number of branches 
points in the progran. 


The latter 
at a few 
| 


The techniques described above would occaSionally fail 


in intervals containing a zero in the denominator of the RHS 
of Eq. (A.39) for 5 (8). The point to be avoided is found by 
solving the equation D(®) = O where D is defined by e 
(A.46) The solution is si 
8 where 89 = vith. Then 


quation 
mplified by substituting Oyte for 


D= “apsin(09*¢) + Bocos (6, +z) + YK 


= -apsin 8 Cos € - apcos 85 Sin € 


+ 


Bpcos 8p cos '¢ - Bpsin 85 Sin 0 + yx 
ave Now kv = +7 So cus 89 = -cos Vn = \cos ¥o and sin 6, = . 
“sin 3) * -sin vy: It then follows from Eq. (A.14) that 

e@ cos °0 + 8 sin 89 = 0. Therefore, when D = 0 


se YK = . YK 
ites c= p(asin 85 - Bcos 99 pu (A. 49) 


If |y«/pu{>1, then f(°) is never 


perpendicular to b, 
but Otherwise, B-F(9) = 0 for 8=6y+c. 


At those two points, 
Plane (8) defined by 
To determine which 

Side of T(E 9+z) contains the ray, we calculate d,=T(99+r). 
(By -Hi(35+2)) which is the distance from the plane to the ray 
except that d 1S positive or negative according as the ray 
is or is not on the side of U(Sp+c) to which fF points. 


s, (5) and B(8) become infinite and the 
Eq. (A.38) becones parallel to the Yay. 


| | 
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The next question is how close to 69+ can the inter- 
section (if any) come. The following approximation will be 
adequate for our purposes. Recall that the radius of curya- 
ture of the helix is ptK /p. The intersection Might be as far 
as pta from the axis or ptatx?/p from the center of curvature 
which would minimize the distance from Oyte to the intersection. 
In that Case, if the plane moves a distance d, to go from 
85+ to the intersection, ® must Change by more than 

& 
0. = _prk*/p ae dy (7+x? ) 


a A.50 = a 
ptarx*/p (p+?) p +k +pa ( ) 


where (p2+x2)% = Jt (6) | = [dH(0)/d0], the distance a po#nt on 
the helix moves per radian\ increase in 9. Thus, of 8nin?? we 
look for intersections between 6*c+O,; , and Oot" and if Sin? 
we look between 8 and 8 and Te a In either case, B(é) 


is well-behaved over the range. 
. \ 


\ 


Similar image ne apply to 69-5. It should be 
emphasized that if O Lin is véry small, then b-F is small and 
the problem of Searching for roots of B(@) are avoided by 


using the cylinder approximat on. 
7 
4i, 


A.1.15 Assembling the Pieces : 
: | 
Having considered all the components of the procedure 
for finding the first intersection of a ray with a helical 
Pipe, we proceed to describe how, they are put together. The 
flow of the Calculation g0es as follows. 
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A.1.15.1 Step 1 \ 
‘Fundanental constants (b,uy,A...) are computed. Six 

special values s., 1<i<6 of the parameter s are evaluated 

together with n;=n(s,) and Ny = n(0). The number n of the first 

turn to be investigated is selected using the following 

decision rules. / 


l. if $170 set n = n, ¥ 


\ "ae ‘if $,<0<s, set n = No 


oe BE $2<0<s, set n= ny 


i 


4. if $3<0<s, set n = Np 


Ri: 

se 
5S. a£ S4<0 then there is no intersection 
6. if N=Ng as set by\one of the above, change 


n to Ne. 


f \ | 
A.1.15.2 Step 2 ! ‘ | 


Having selected a turn fot investigation, the method | 
of Al.12 is applied to determine ‘whether the line of the Tay 
comes close enough to: either half tof the turn to intersect 
the pipe. If so, go to step A.1.15.4 otherwise go to step 
A.1.15.3. \ \ 

a \ ; \ 
A.1.15.3 Step 3 __ \ : 

It has been established that there is no intersection 
with turn n (s<0 does not fount). ‘ . 


| 1. if n=ne<ne, Set n=n¢ and go to A.1.15.2 
otherwise, set n=n+] and continue. 


\ : ‘ 


| 
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2. if n>n, there is no intersection. 


3. If nj<n<n 


2 3? set n=n.. 


4. Go to step A.1.15.2 


A.1.15.4 Step 4 


An intersection has been identified. The value of 
B-F is determined at the closest known pair of pdints between 
the ray and the half-turn and on that basis, the B-nethod or 
the cylinder approximation is used to determine the location 
s* of the intersection to within an acceptable tolerance. If 
0<s*<s** where s** is the best previously determined inter- 
section (initially; s** = 1071) then set s**=s%, Both half- 
turns are conside.. d. If n=ny<n¢ then set n=n¢ and repeat 
this step. Otherwise, one or more intersections with s>0 
have been found and the procedure for ordering the turns for 
consideration assures us that the minimum of these positive 
S's must correspond to the first intersection. 
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APPENDIX B 


B.1 AUTOMATIC SEGMENTATION 


A separate computer | program has been developed to auto- 


matically segment a generalized circular cylinder. The code 


is called LAZY and produces data input cards for the ZAP code. 


It is assumed that the cylinder is horogeneous and is 
confined between two end planes. The outer surface of the 
cylinder and the two parallel end planes are considered to be 
Fresnel's reflection. All internal cylinders and planes are 
treated as dumny boundaries. 


Either a full or half cylinder may be generated. If a 
half cylinder is used, the plane through the center of the 
cylinder is treated as a plane of symmetry. 


An astion exists for treating a cylindrital shell. In 
this case the segments are those confined between an outer 
cylinder and an inner cylinder. 


The segments are formed by generating a specified nun- 
ber of equally spaced planes between and parallel to the two 
end planes, and a specified number of concentric cylinders 
within the outer and inner cylinders. In addition, planes 
are generated at equal angles through and parallel to the 
cylinder axis. In this way each segment is confined by four 


planes and one or two cylinders. iz 
In addition to generating boundary and Segment, cards, 


the program will generate structure data cards. 


The cards generated by the program will be only a por- 
tion of the cards needed to set up a ZAP problen and are 
intended primarily as a means to Segment the laser rod. 


| 
} 


| 
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B.2 INPUT FOR LAZY 
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‘ \ . 
Input to the LAZY code is in the fourm of dzta cards. 


Standard Fortran Data keypuiich forms are used :nd formats are 
1216 and 6E12.4. 


Output from the code consists of punched cards and a 


listing of the card images. 


Care must be taken that no plane passe+ through the 


' origin. 


1) 


2) 


3) 
4) 


6) 


(Card #1 (1216) 


IOPSYM 


NBAR | 


a 


? 


f 


Option for generating only half of 
a cylinder with a plane of symmetry. 


Plane of symmetry with segments on 
non-origin side of plane. 


Full cylinder, no plane of sy eury. 


Plane of symmetry with segments on 


"' origin side of plane. 


P | 
Option for generating nother segments 


and mother-daughter structure data 
cards. 


No 
Yes 


First boundary number to use 


Number of slices to generate (MBAR+1 
parallel planes are generated). 


a 


Number of rings to generate: (LBAR or 


LBAR+1 concentric cylinders are gen- 


erated). 


Number of wedges to generate (for a 

full cylinder, NBAR must be an even 

Number. The number of planes gener- 
ated is NBAR for a half cylinder and 
NBAR/2 for a full cylinder. 


B2 . 


’ re # 
7) JPAR = First -segment number to use. 


8) MAT- 


Material number for segments gen- 
erated. . 


9-12) Not used 


Card #2 (6512.4) 


1) a Directional cosines of the normal to 


2) a’ _ end planes. If the aistan.e to either 
p end plane is negative, the negative of 

3) YD these cosines is used. 

4) a. 

5) B Directional cosines of the normal to 

Fs s the plane of symmetry. 

Oy %. 


Card #3 (6E12.4) 


1) xXx 1 . #3 
2) ye Coordinates of 2 point on one of the 
pl end planes. 


3) 251 

4) X,2 . 

5) P Coordinates of a point on the other 
p2 end plane. j 

6) 


Card #4 (6E12.4) | 


1) a 
2) a. Directional Cosines of the cylindrical 
c axis. 
oY, 
4) X.. 
5) Y. Coordinates of a point on the cylindrical 
6) 2z axis. $ 
c 4 . Sw 


BS 


Card #5 (6E12.4) 


1) RC = Radius of outer cylinder. 


2) RI Radius of 


3-6) Not used 


inner cylinder (usually 0.0). 
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ZAP FLOW DIAGRAM AND PROGRAM LISTING 
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Paocnan ZAP 
Cocccepnecocececcocceseccvseucecoccesces eccencocencesenccesesesooneaeece 
‘Ceeectevee Zap 1S Tre NAIM PROGRAM AWD [5S PRIMARILY A CONTHOL REGION, 
CO ee cececescecesseeseeeesOed Se eeeSeSeleesneeeSeESOSaSOSSSeesoceesasre cee 


Ceeceetese 

Ceecescces ThE MEAMING-AKD- MAKINUN SIZE OF THE SUBSCRIPTS IS AS FOLLOWS 
Ceoesetese ‘ 
Ceeeeesese F WAVELENGTH rai) MOTE THAT F IS INTEGER 

Cesesesese 7 STORED RAY « 

Ceecesesoe 1m [PUT RAY see 

Cescecetese 3 SEGMENT 200 

Ceesscossce «x BOUNDARY 2ee 


Cee-cessee CP QUURDRY PHTIR 900 
Cececceeeee 4 REFLECT TYPE 286 


Ceeceovece Lamp TVPE 2e 

Coecceecsss ¢ * MATERIAL 2e 

Cecesesese / ; SOURCE se 

Ceeccesece *' 

Ceecccccce Tie MEANING OF THE VARIABLES Iu COMMON IS AS FOLLOWS . 
Ceesestocoe 


Ceecsstee A8CO (oF) saBSORPTION CUEFFICIEmT OF MATERIAL © AMD WAVELENGTH 
Ceecccrces ats ALPHA OF RAY BEING FOLLOWED. 

Ceercceces alK(K)s KX DIRECTIONAL COSINE ASSOCIATED WITH BOUNDARY Ke 
Ceecceccce ALNZALPHA OF RORMALS 

Ceeresecce AL PuACTIS KX OIRECTIONSL COSINE OF STOREO MAY J. 
Coeneerroe ane ALPHA OF REFLECTED CAYo 

Ceoseece -¢ ALRAY(IN) SARRAY OF ALPrtA FOR [MPUT RAYS 

Cescenecce AL RFs ALPHA OF REFRACTED RAY, 

Ceccosscce ges RETA OF BAY BEING FOLLOWED. 

‘Ceeceeeoee BEKIK} © Y DIRECTIONAL COSIME ASSOCIAILD WITH SOUNDARY Ke 
Ceecceteope SExsSETA OF MORMAL. 

Coeccetece R= BETA OF REFLECTED RaVe 

Ceeresccee SERAY(IN)SARRAY OF BETA FOX IMPUT RAYS. 

Cescseccce RERF=s BETA OF REFRACTED RAY. 

Coeccececoes BETALI) = Y OIPECTIONSAL COSINE OF STORMED May I. 
Ceseeseooe CCORDS(IMDeM)s SIX COOMDINATE LIMITS FON SOURCE Ne 


—~—€sveseecoe PEVIJIS DEVIATION OR ERROR FOR SEGFENT J, 


Ceseceetoce NISs DISTANCE FROM FO[aT TO INTERSECTION, 

Ceeccosose Ox(K)= DISTANCE ASSUCIATED WITH ACUNDARY Ke 

. CBeveetecces Qwyns Mimsy JNTEMSECTIOCs DISTANCE FOUND. 

Ceeeseseoe O2s DISTANCE BETNEEN TG? AND BOTTON ROUMDARY. 

Ceoeenerene E(t oF}= EMERGY OF STORED RAY 3 AT WAVELENGTH F. 

Ceescenesce EMASE(F)= BASE EXEQLY FOR RAY AS A FUNCTION OF WAVELENGTH, 
Coecceceece Exqk)s ECCENTRICITY IF BOUWaRY K IS A ContCe 5 

Cesceersce ELaniJ)= SCALED ENERGY DEPOSITED IM SEGMENT Je 

Coeesescce ELAMATIM) = SCALED EMERGY FOR MATERIAL Me 

Ceeeeesece EMATIN) = ERERGY DEPOSITED IN MATERIAL Me 

Ceceeseece FAERGYI J) = SEGMENT ENENGY. 

Ceescerece EPSs SHALL DISTANCE FON MOVING RaY POSITION OFF OF BOLNDARY. 
Ceeseoeres ap agen lle RAY ENERGY FOR LA@P TYPE LN AnD WAVELENGTH Fo \ 
Ceesceeoes Fraxs BER OF WAVELENGTHS (INTEGER). : 
Ceoccesese Gis OF RAY SGEING FOLLOsED. ‘ - 
Ceaceeeece Gax(K}s Z\ DIRECTIONAL COSINE ASSOCIATED ¥ITH BOUNDARY Ke 
Ceeccooves GAMMACIDS ¥ OINECTIUNAL COSINE OF STORED RAY I. 

Coccccerse GancGaml NORMAL» 

Ceocecese-ce Gitis Gants REFLECTED Ray. 

Ceeseceeces GARAY( 1K) SARRAY OF Garna FOR lePuT RAYSe 
Ceoccesccs Ganfs GAnnA OF REFRACTED RAY. 

Ceecceessce IFLAG(I)S FLAG IMDICATIAG IF STOREO RaY. 1 IS LIVE GR DEAD. 
CHeeeseooe IPAXs NUMBER OF. RAYS SIORED. 


C3 


 fcecmetaa | 
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pees 
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[= 
L 


1nOws PRESENT RAT INDEZe 

1OPazTs OPTICS: FOR PRINTING Each AAV CRESTED. 

IPOIKTIN}sA POINTER POIATI“G TO FIRST STORED RAT FGR SOURCE ts 
IRs LAST S@AADGM INTEGEH USED. 

Tals RAV MBER BEING TRaCECe 

ISMAR(N)s RUYBER OF RAYS TO CAEATE FOR SOURCE me 

ISORILN) = INDEX OF SGUACE TYPE FOR LAF? TYPE Lite 

JS (J)s FIRST DAUGNTER OF POTHER SEGHERT Je - ; 
WLtJ)s FLAG IMDICATIAG IF Gav DEPCSITED EMEAGY Im SEGMENT Je 
AL (Jis LAST DAUGHTER OF MOTHER SEG"ENT J, : 
WUJ)s MOTHER OF DAUGHTER SEGMENT Jo 

SPARS WAKIMUM SEGHERT BUMGER, 

St Jie KET SISTER Ix mE IRANCHYS 

JKOUS SEGHERT RAV IS 1m NOde 

JCLD= SEGMENT RAY WAS Ihe 

WPtJie PMEVICUS SISTER Im nETRAACKY. 

Bav(s)s MUPAER OF RavS wnICe have PaSSEQ THROUGH SEGMENT Je 
JSOR(M)= SEGMENT RUMBEH OF SOURCE MIPDER fe 

KER (KP) = BOUADARY AUMSERS FOR SEGMENTSe ; 

KFLAG(K)= FLAG IMDICATIAG IF GouacaRy K MAS ALREADY BEEN 
KINT= GOUNDARY BEING CONSIDERED FCR INTERSECTION. : 

KYAEs PAL INUN GOUADARY DUPBER. | 

4ACes BOUNDARY WHERE RAY IATERSECTED. 

KPIi= 2MDEX POIRTING TC GOUNDARY HUNBERS FOR SEGPEXT Je 
KYV>‘Kie GEOMETRY TYPE FOR SOUMDARY Ke , 
Le-timde LarP TYPE suvER FOR SOURCE AUMEER No 

LecLE(LN)s IRDER GF ANGLE ROUTINE FOR Lap? TYPE LM. 
UREFLIK) © REFLECTION AUPBER FOR BCUMDARY Ke | 

LIVPECL) © REFLECTION IVPE FOR MEFLECTION MUMBEN Lo 

Salt Je MATERIAL WER FOR SEGMEAT Jo 

ELEM =FLAG INDICATING IF Ray DEPCSITED EAERGT Im maTERIAL Me, 
Prats PALIN MATEAIAL RUPSERS 

AvaT(P)= RUWSER CF AayS TrAT PASSED ThROuUGH MATERIEL M, 
AAOes SQUSCE KOW SEING LSEDJe 


RSORs “ARI@LY SOUSCE AUPSEXe 


ALPQAT( J) 2-2 °dEM OF GAYS IraT PASSED THROUGH SEGFENT Jo 
PER= PERCENT OF RAY TRANSPITTED. 
PTDtLeF)= REFLECTIC’ OATA FOR REFLECTION L AAD VAVELENGTN Fo 


@ Scr(mpe IMDER OF BEFRACTIOW FOR MATERIAL F, 


St2waTimps TERS MEEDEC FOR ERROR BXALVSIS. 

Siva = SUTTATION OF GAY EHERGYe 

Sivaavig)s TOTAL EXERGY DEPOSITED Im SEGuERT J BY ORE RAV, 
SuMSORL) = SU" OF RAY EXERGY SQUAKED IM SEGPERT Jo 
TIVLECI2)© VARIABLE ALPrARUPERIC TITLE SntCe DESCRIBES RUA. 
weaTtM)s VOLUVE OF MATERIAL Me. 

VCL{ Js WOLUYE OF SEG>EST Je 

ZUI)s % COOADIMATE OF STORED Rav Io 

Riz. POIST USED FOR FInOING POSITICN AbD IATERSECTION. 
AK(K)® X COCADISATE OF POINT ASSOCIATED wITH BOURCARY Ke 
RLAS(F} = ARREY OF BAVELENGTrS. 

EAQss &% OF ISTERSECTICS POIATS 

ECLOS % OF STAITIAG PCIAT OF Sav, 

EQAVEIA)s ZX COCROISATES OF INPUT BATS. 


EGEFs Z OF REFLECTEO GAVe 

RTgaMe X OF TRANSMITTED Rave 

VEI)s Y COOADTMATE OF STCHEO Rav To 

Vis POLAT USEC FOR FIACING PCSITIC’ AAD IATERSECT ION. 
Veends ¥ COCSCIAATE OF PORT ASSOCIATED WITH BOUASARY Ke 
waOes v OF IXTERSECTICA PCINTo 

YCLC® VY CF STANTIAG PCIAT OF MarT. 

WRAVEINDS VY COORDINATES OF INPUT \ATSe 

REF s VY OF REFLECTED Rave 
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sae 
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~ Ceoeeesoce yimaxs ¥ OF 1 ITYTEO may. 
Coecceeoee Z63)=8 2 COOZD] TE GF STOSED May 1. 


Coecscscoe 723s POINT USEC) FOR FInOIA6 POSITICN ano INTERSECT Sox, 
Ceeeesecee Za(x)s 7 Cocacraate oF Point associatEep wits Eolacary 
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Ke 
Ceececesoe 7: saminda LASER WAVELEAGT" OF maTeaiag nm, 
Coocecesee ZaCe= Z OF IRTERSECTICN POIRT. 
Ceeeee-see Zope 72 OF STARTING POINT OF may, 
Ceeeccoces Zaay(Ixjpe 7 INATES OF Input Rays, 
Ceeccesose ZGEF=s 2 OF REFLECTED maYve 
Cesecesoune 772400 Z OF TRAYS=ITTED May, 
Comemgu 7697 (200290) eALRAY (500) oGCRAY (S08) 
I eCOCASS(6.50) oF 6400200) oEBaS2 (200) eESOR (200290) 
2 GANA (565) RED (200200) oTZTLeE (12) ovat (29) 
3 evor(zee e (200) oZRAY(S00) oYRAY (See) - 
@ e2MAY(500) ; 
| COmme~ saRRsvS7 pas oAlLPraqéep BEX (200) ! 
{ E oF T0149) eDEV (230) eOX (200) ‘ e€K (209) 
| 2 e€Ua- z20@) oEL Aras (29) oE™AT (20) eE®EB6y (206) i 
+ 3 eGax«zees oGar~2 (40) oREF (20) eSOetat (26) | 
© oSuPmavi200) «= oSumscl1zee) oK(4e) oKK (286) ' 
5S evtee) o¥K(200) o2(ae) o2K (200) | 
© e%laxt29) 
taRmays IFLAG (60) oIPOINT (SO) e 
@lSmax «Sep o1Son(2e) oF (200) oJFL (200) ® 
@0L (200) oJ" (200) odm (200) eS? 6200) ° 
@JNeY (260) oJ$02 (56) ,0RBAD (990) oXFLAG (290) e 
@nP (260) eKTYP (209) eLarwiSe) oL ANGLE (26) e 
*UREFL(2003 «= eLTYPE(29) omaT (Zee) oMFL (20) ° 
Lee LATS o&uraay (200) 
ARES 7s a olin oR oALaF eet ° 
Ste SER o@kar 0IS - —onIm 4 EPS ln 
6a eGam e6aa eGanr ° ePER ° 
eSueaa = gE eRACe eKCLD oXRGF okTaan evi ® 
evans e¥CiD oVREF oVT@Qam oZI ° 
eZwne e2CLD oZatF oZTRan 
Com wy lPasaes Fmax olmax eixow © 
*lopaay .icmmar gta olaay oJnak oJupt eJOLD 
@ayaT okPax okay oMmMAk 


ote oNSOR 
IuteGEa FeFrax 
Coecccccce Te 12-56 a mBERED Com—On LOCK IS FORCED 1p 0640 OvER Ine 
Cesecesece LeAreP 1% pant Be Tus KEEPIVG Int PuCGlars Ih wang De mugen 
idindiaas 3S FastEa, 

Baa (0200 1697 
Ceeeceecee necarse Tx MUPBEKED CiPPON MLOtee WLdag O'-48 2¥eea.e. s3 
Coweeseese acy SET 19 ZERO OW LOA0INGs IT must ag Sracateue Wu-es, 

OFeEnSICK Cw~v 123032) 

EQUIwALERCE (Lm) o£aCOr3)?D 

Tee @ TIMELEST OND 

89 106 [Ae] 523032 

Ouvenv(34) 2 @. 

Ld ive 
Gtbccctsce mech Tapul Cara. AnD SET Up PRON EH. 
Tue 

giasianes ae Te€ SExT May T2 RE PaCCESSES, 

Caut Saycleice 
Geateccake IF TeE=€ a-~€ 20 “CRE BavS 6 TO 382. 
ctsocseces Ge TpiiaE anton SEGMENT CenTanss STATTING POIST OF Ray, 
cosencsans IF ste~tat NUPOER IS ZERO Tr€ way IS CUISIDE IME SYSTE~. 

Ie d106,0 00 de 7 

cabaseects ie POIST v-€2E “av JATEM@S. “TS Te€ MEAREST BOWERY. 


aes 


: hee 10% Cart sano - 
‘ Ceeeeccces DEPOSIT ENERGY Im SEGPEAT RAY WAS WUST ClS55E0. 
POS 


; - Car O€ 
| Cowccccsce CETEMNIRE WHICH SEGHERT CONTAINS TRANSMITTED Ray, 
: Gomst ; 


CML SF 
Coceeeesece GET OIRECTIOCNAL CCSINES ano EME“GY CF REFLECTED &a9 
| , Ceesessoce TEanSMITIEO RAYS, PLaY &USSIay ROULET iF £liver R-%a5S 
' Ceoccesece fagacy LEVEL IS SELOw Cmosen VALUE e 
905 Cart scunce 
_ Coceeesscee OF TERNIRE IF nEw May 3S *EEDEO, 
{ COLL OECIDEIIFLG) 
IF CIFLG) LOT lOrei08 
102 Convince 
Ceseseseoe Paint £ TAL (OM SEMZoF sna) ANSvERS,. 
| ~~ Cau €orr 
| Tie = TINE « TIMELEFT C3) 
PRINT 106 TIE 
te corsa iin oLIMERECUTION TIPE 1 oF 8.30GH SECONDS) 
oF 
Cecccetene ie COmDITzOm. RAY HAS ESCaPEeD SYSTEM, 
‘to. Tee Caul Expmurion — Zap) 
‘oP 
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PS 


a se ; \ 
| : / 
f leds 
SvazouTy ANGLES (TH gPngy  / 
Coeceoecceosens secsacseaseccccsecescccncecoccccesecccecccecencs t@ePeocee 
Ceascestece TIRE ANGLES CAULS Int DESiaEep ANGLE ROUTINE, 
Cooeoe nner cceehecanssccscosesesevecasecescesocones cent 
[ Coseo:, 7697 ABCOtZ2es200) = sa mav(5@9) sSEAAT (509) 
/ 1 eCOC3LS 16050) 0 E14600200) oEBaSE (208) sESOR(2¢ e200) 
: 2 eGaRay1S00) oRED (200200) eTIYLECI2) oUMal (29) 
3 evo (Zee) oXL 4m 1290) oXRay(5e0) eYRal(Se@) 
4 eZ8av«See) ei 
i Corson saRRayss ALK(200) 042 P8148) oBEK (209) 
1 oPETALée) oDEV (200) 00K (200) oEK (209) 
| 2 e€t 01208) > EL Anat ¢2e@) oEMaT¢ 20) oRNERGY (269) 
3 eGaxi2ee) \ @GArrta (48) oREF (20) oSQWMaT (29) 
4 eSu™=Ray i200) oSUX56R 1299) oX(ée) oXK (200) 
j S evtses eo¥K (20%) 02 (ee) o2K (206) 
| 6 027i 4%(2¢6) 
Conewgns | anmays FLAG (408) 0 IPOINT (58) e 
©1S~4x(5@) 03S0a(2¢e) oSF (200) oJFL (200) e 
©. 1200) 0200) own (Zee) oJP (200) e 
{ @JraV (200) eJ50R (50) oROeD (908) oRFLAG (290) ° 
en? (290) oS TYP (Zee) ola? (Se) oL ANGLE (20) ° 
| @LREFL (206) ol TYPE (20) oMal (200) o™FL (20) ° 
otal i2@) oMUPMRAY (289) 
BRANSS a OAL oALR oALF ome ° 
| REX o8ER oBEQF eDIS oDnln 00Z ——sEPS 5 
i 6a oGan e6ar aCe oPER ° 
CSu--e4 KI oX™Ow oZCLD oXREF oXTaan evi ry 
o1ma eYOLD oVQGF oVitan ry 4 | ° 


e2unw eZ0up oZREF oZTRAg 
Com Ouy 1P Anas Fman olpax o Inow ° 
| @loraay e10°ant oiR olaay oSMaAx oe eJ0LD » 
eKjarT oKRAx oKwOe ommax . obmenN oMSOR 
Imtretm F oFoax . 
Ceeceeccese Cex Lavo MUNBER FOR SOURCE tuOw. 
LY at 4p (aa) 
coe GET MUPBER OF ROUTINE TO cau. 
L4sL ANGLE (LT) 
°° GO TO DES1-=n Cau. 
60 To {100 o2C* + 30206000500 0600.700 08000000) ola 
200 CALL ANGLED (Ie, 207) ; 
60 TO lle 
200 Cart ABGLE2 (TetoPHz) 
Go ‘to rile a 
300 CALL ARGLE3( TIN PHT) 
‘60 TO lle 
400 CALL ANGLES (TH.PHA) 
60 To lle 
See CALL AAGLES(THPH}) 
60 70 lle 
600 CALL AAGLEG(THt.PH}) 
6 To lle 
700 CALL ARCLETITH.PH}) 
GO To lle 
Gee CALL AAGLESC THOM) 
GO To ile 
000 CALL AAGLES(THoPH}) 
lle R& Turn 
Emp 
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SURaouUT INE ABSORP (748.0%) 
Cocvcesesccsencceccccocccocoscoccccccccescceccs. CG4 Se2 250806 POGeCOEL CE 
Cocooecees SUBROUTINE ABSCi? CALLS THE DESIRED ROUTINE FGA dtazazrtng 
Caeescecece ARSOAPTION COEFFICIENTS as A FUNCTION OF wAVELENSTH, 
Coeeseescocosceocecsece CPecnsseereseeosescoesercoeecsete @ocaseecon: 


Comings 7697 ABCO(200280) = satmayisee) 
I sCOORTS(6055) o£ (400200) sEBASE (200) 
2 eGARAV (586) oME0(200208) eTITLE 12) 
3 ever (200) oXLan{290) oxRay (500) oTRAY (50z; 
' @ e2MaT (500) ‘ 
| Cowpuou saRRAYSs ALK(200) oAL Pra (40) oBEK (200) \ 
' I oMETACA@) oDEVI200) o0n 4200) oEK(2ge0) = \ 
l 2 cELAN«200) sEL AMAT (29) oEMKT (20) vEMERGT {3- :\ 
3 268x(200) oGArw-a (48) oRES (20) oSQQMAT (2-5 
© eSimRayi200) sSuptson(2¢6) oR(a0) eX (200) 
S evises o¥K (200) oZ(4e) 02K (200) 
6 eZLani2e) 
TaRmays IFLAG (40) o SPO1KT«(S@) e 
@lSuax (Se) o3S0R (20) oF (200) oJFL (200) ® 
\ OSL (200) | 011200) oJh (200) oJ? (200) ® 
| @IMaYi200) 0 JS0R (50) oKBAD (990) oXSLAG (200) ° 
\ om (200) oR TYP (200) olan? (59) oL AKGLE (20) e 
| CLMEFLIZ00) oh. TYPE (20) omaT (209) oMFL (20) e 
| ofa T ezes sRUPRAY (200) 
| *CommegeP amans/ au oALn oAlR oALRF oME ® 
' e0En ; eSER oBERF 901s eOnin oDZ eEPS ng 
{ @e6a Gan oGan oGanF oPER e \ 
{ e@Suwen our oXnOw eXOLO oAREF oXTRAN oVI e 
7 @Ynow eVOLD eTREF oT Totang eZ] ® 
\\ @Znou oZ0L0 o2hEF oZTRam 
i rPanans Frak o Tmax 9 Btw e 
) @lopmay .f0PRnNT oR : oImay oJMax oJngw eJOLO » 
eXiat oKax eERlel onwax oMthOe eoMSOR 
INTEGER FoF max 
OIFENSION TST61209) 
; SIENSrOw 1eptasi9). 
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Oatac 148T49=60143S) eGMABS2 0674853 


e 
é lence’ rayne or neste anna eCHAESO 66n28S9 » 
IND MMUPRE = Rout TO Cau 
00 161 wel, : 


| IF (7a8—-]asTagins) 10101020301 
COnT Seg 


e¢ ROUTINE MAME KOT Im TABLE, 
\\parer 160348 
10, F ossaT (M1 e21MaaSORPT Iom SUBAOUT INE o 1x AG i 3H nOT In TABLE) 
cer 


Cc we GC TO DESIRED CALL : 
102 |Go 10(100020003000400050006000700.80C,900) om \ 
106 t 48S) (XLAMoF MAX eTSTG ond ‘ 


10 
Cc 
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To 1le : 
260 COLL ABS2 (cL AM, FmaKe YSIS on) 
£0 To 11e 
309 FALL ABS3CKL AMF raxeTSTC om) 
0 TO lle 
£00 CALL ABSA CILAN.FmAxeTSIC om) 
To 130 
300 CALL A8SS (aL aM FmaxetSTS er) 
| To 110 
\ 608 CALC ssSetnLanFwaxetste om) 
To ile 
Vee CALL A8S71xLaMeFmaxeTSTG on) 
$ re le | 
| 000 CALL COSStZLAM.FmaxetSTC ond 
To lle 
(900 CALL AASOLILAM, Frans ISTSG on) 


lle Coatime i 
c Oe. sg ORE AGSOSPTION COEFFICIENTS FCR MATERIAL NO, My 
| 00" 103 Fal oF max 
(1F OTSTC (FS 
103 Proll 
eRe 
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PUBROUTINE OCALC 
Conesepeceseececesceccoegses ccueeseoueesosese TT Sa 
Ceecessece BCaLC CALCULATES FUNCTION BITH)» CERIVATIVE BP ETH)» am 


Ceoseeesos OISTANCE SinsTet eo \ “0 - 
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& eFRAat (500) 

Cor=0m saRMAYSs ~ALK(200) oAL Pua (46) o8€K(200) 
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2 ef0an¢200) cEL AMAT (20) Ematc20) oEMERGVEZOR) . 

3 eGax(290) eGARAA (46) RES 420) oSORnat 12¢) 

4 eSuU"MavIZ200) © 0 SU"SORIZ00) oX(40) oRK(200) 

S evcse) eYK(200) eoZlee) o2K (206) 

@ eFLanczer 

Cones 1 ARRAY IFLAG (40) olPOIMT(5e) “i 
©XS0ax(50) eTSOR (20) oSF (208) ouFl (200) . 

©A. (700) 0086200) owmt2ee) oSP 1206) . 
@JMavi2ee) 0 JSON 150) Kae 1900) oXFL 061200) ® 

nr 1200) oxTYP (200) oL ahr iSe) eL ANGIE (29) ° 
PLREFLI200) «= o LTYPE (28) omati20e) oMFLIZ0) . 

@fona T1290) oMUMRAY (260) | 
COmnOusPaRanss aL oaye oALn Ale ont . 
OE ms oBER eBERF 07S hd 002 PS ls 
e6a 96g eGar oGang PER x) : 
eSuwan KI oINOu eXOLD omnes =| ekTAAm gE 
eTnQw eYOLO eYREF oY TRan eZ] i) 

ezZuou eZ0LO o2ktF oZTRAn 

Cow~gus 1 PaRans Fnak olpax oluge ° — 

®10°2ay fIOPRNT = g IR olaay oJeax oN oo OLD oe 
@XIaT okmax eke olmak o8De oNSOR 


IsJEGER FoF MAX 
HEL COW 74 o ABAR SH eBI 00 IP eBPoCAPoDEL THe DEL RX oDELYeDELZ eDulT, 
OPS1Z oREV RUNDE STHe THe THY oSORKS oARS CS Tuo SMTH ORCS THREAT cOSTie 
@COSTAMe THZ eDINZ 
Ceeseescee CALCULATE PARAMETERS, 
SToweSia rT) 
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OO 163 Keolycnzx 
HFL EG (4106 
103 CoOaTince 
102 CAatIwe 
RYweXK To] 
Ceoecccroee C2 ALL wIAIZCATAL FAYSs 
OO 106 Lelsay~ 
w1s0 
Ceeeeeesoe SIT STLATING POINT CF HAY, 
Rh exk Tale (JoeS) OSS TORYK 
Visvielte(Jo.5) *SSYesavyY 
ZleTT ele (Jo.5) OSSVOAyZ 
AU eSEK 
GEesiv - 
Geotx7 
OLes.6 
OD | 
Coeceeetcce Finn SESPENT COATAISING POCIAT. 
133 CLL SES~NT 
IF ¢ ACodLTOsSTIGSACS 
Coeeooeoes FIAT SCUNDARY INTESSECTICNS 
10S Cace 87U45 
DO 325 Kat emuan 
BFL S6(<) 58 
106 CCaTca_se 
Ceeeeeeoos SFY PRIAT CeeBSR2ACTEA Flat WATEQI EL ALPSERe 


MseTC (het CNCole] Jase) 
MPATs 1 Cres tee) 
Ors %Lerer% 
Coeccecrcoe JF OFSTASCE TRAVELED GHEATER THay LENGTM OF RECTANGLE, 
Ceeeecesesee GO TC 106. 
IF coe e¥L) 19701670108 
Ceeeecrooe CrCiLATE PEST COOADIABTES. 
WOT MAeVsAsy CAT -AZERCI CALVO CY POTZERC) oarZe(ZJo7ZLEAO) 
VARI=kvs LH ATEQOICAVVO(YJoVZERO) oavZe(ZoZLERD) 
Ceeeestece Sii ri uty RUMZER OF AUOUNTANY IATERSECTICA, 
O22 (NE 28/55701.5) 
Ceeeceseee SET PRIAT CuARACTER TC BLANK, 
Cine an2ypsko 
MUNEPza Je] 
WT." vtX\2—e] 
LF ce “SE PONT AON) 10901090112 
Ceescoeooe STORE MATEQIAL NUMAER IN ALL PRINT POSITIONS BETWEEN 
Coessecsce CCAS. CUTIVE INTERSECTIONS. 
109 09 1113 n=.CrEPnTuCc™ 
LIAE (Rp) erat 
Mlle CCATIAVE 
112 Nle%2 
Ceaeescecoe GC LOI“ FOR NERT INTERSECTION. 
69 TO 113 
Coeeessooe SICH. MATERIAL AUMBER FOR REw.rNDER OF PRINT POSITIONS, 
108 AGAEP cA} 6] 
LF CNOSEP OKIVIGCL1460116 
114 OO 115 wawCrFP oun 
LIAE CX) owPRT 
115 Coarr.ve 
116 ALceava = 
BEa-LvyY 
Ganeav7Z 
Ceeneseeee (CCP TkQQuSH ALL VERTICAL RAYS, 
O09 117 Nolonx 
Ce. seeeoee arv—E RLY STARTING POINT BY Crt LIAE WIDTe. 
XCOL (ON) SKCOL CN) eALODSY 
YCOL CN) CYCCL ON) CBE COSY 
ZCCL C4) EZCOL EN) 0GAOUSY 
OCS CN) =CCOL (8) OST 
Ceseeceecce IF INTEASECTICN WAS ACT CROSSED GO TO 1176 
AF CRCOL (NID LIB e1LITe! .7 
116 XjencOy(%) 
FievVCOL (%) 
ZleZCOr tw) 
INO az SCOL fh) 
COeeceeece F1X) KEW SEGMENT KUMAER. 
Cart SEGYNT 
IF C800 913901100119 
Ceeseeeece FIKD AEXT BOURDARY INTENSECTICHe 
119 Cart s8cuwd 
00 120 Kel oKmAx 
KF. AG (K) 08 
120 CoaTIweE 
Ceseceectece SET BLAWK FOR PRINT PCSITION, 
CINE (sy sla 
OCrr (md EDMIN 
117 CoaTis € 
Ceeecetoce PalnaT LISE 
PRIST Le Lint 
10 FOe~AT(im o13¢43) 
Ceeees*ees GO 00 NEXT HORIZONTAL RAY, 
@ 


104 CCatiace 
Coeeeeeece GO LCCK AT AEXT DALA CARD, 

GO 10 123 

121 Jwt5e} 
RETUAaDK 

11@ Care E®PANT (6rNGEOPLT) 
Stro 
Exo 
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SUS2CUTINE HEL NDA: 
Ceeeeetooce SHSSHOSSHSSCHSHSHESHSSSESSEOHSESE SDE LSEOS SOSHOCHSEHS. ceseseasesete 
Ceeeseroce Hii .CR CALCULATES CIAcCTICN&L COSINES FCR Tré »Camat 19 
Ceeeeeeevee ¢ POIAT Om Tei SUMP ATE CF A mELI ay 
eT rir ii iriirrririrrrrriry) 


Correos 7697 ABCOI2507-5) o2l247 (500) eBERZY6556) 
I 0€9047S866050) «=e El 460205: eEFaSE (205) oESCAI200¢2¢0) 
2 eGeSav(500) eREU(25 0250) eoTETLE C12) ovr ali29) 

3 ev%Lt200) oXLan(250) oXRaV¥ (500) eVSAV(S595) 
4 o7=aveSo0) : 

COse0N sERRAYSS ALK(260) oAL Pra (66) o5Ex (200) 

1 eee Tae) eDEV«200) efx¢200) oEK (256) 

2 ef A™(200) oELA~aT(20) oEwat (Zo) oE%E~57 (206) 
3 06AK1209) oGAMnA (46) eoREF C20) oSe** 47420) 
& eSuUMA2YI200) oSLKS5R 1205) oXlad) oXK(226) 

S ev40) 0 ¥K(200) oZ(4G) oZK 1206) 

& e2t a™i2ce) 

Coe~nns TaARRAY/ IFLAG(40) eIPOINT(SO0) ° 
@®I1S~an (S90) 01S0a120) oJF (299) oSFL (200) ° 
OSL £760) oJM1205) oJh(2660) oJ? (263) ° 
®JRaV(200) 0 JSOR (S60) oK@nd (900) eKFLAG(206) ° 
@KP 4709) eKTYP 1265) el am? (So) eLAAGLE (20) ° 
@LREFL (200) oLTYPEt<¢:? oMaT(200) oMFL(26) e 
@r“MaT 120) eMUPRay i250) 

COwMOK/PARAMS/ AL oALK ALR oAL AF oAE ° 
CFE DN BER eBERE o0IS oOnIn 062 oEPS 6 
GA oGan oGar oGakF eoPER ° 
@Syvoa oXI okhOw oX¢CLo oXREF oXTaan eV] ° 
Coxe] eVYOUD oYREF oVTatan ry 3) ° 
@Za76< eZCLO oZREF oZTRAan 7 

Coe-nws]PaRews Prax olwax oluOu ° 
ClO>aay gICZanT =f oJRay 0 MAX ovrcd eJXC « 
@KIGT eae ax oKwOe omeak oMADu eKSun 


InTEEEA FoF een 
Ceceesrcce FIND POINT ON MELICAL IRE CLOSEST TO POINT OA SuAFaACE, 
Catt PCINT (Xe ALPT eBEPT oGAPT oXNCd ov YRDOg ZAC) 
RaStQAT (AL PTee2eREPTeo2eGarteoe?2) 
Ceeeeeeoes CALCULATE OIRECTIONAL CCSIAES OF RORMAL » 
ALK=ALPT/A 
BEasAEPT sa 
GAn=sGaPT/R 
AN 42a Oa e (ROK (K)) OALN 
YNOWZYAOue (ROKK) ) OGER 
ZNOdsZnOU- (R-DKIK)) GAN 
CALL CRECK (ORMELNOR 2p ALA BEN OGAR) 
CALL OERUG (OHNELNOR) 
RETURN 
Emp 
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SUHASCTINE wyPasa(a} 
Cowocecccccessecccascsccccosecscsccseccevocescsececececesevousecceszeeto 
Ceccetces mtanfa FIe_§ Tre STSECTICASL CASTAES CF tre ACSMAL TO A 


Ceoeeeeteee S-10T ON & wYeEMaCLE Ch WEWCLUTIsa. 
Coescescccsccccesccccccccccccccanesasccccccncceeseeonecescecececcesoeete 
COwwt~ 7697 mee 250268} e4Laay(See) eBE2ZAY (590) 
T eC OC5°S$160S6) =e EL ete200) e€5aSF 6266) e£SIR(25 0299) 
2 eG22aviSctad oRED (22-200) oTITLE C12) ow"al (20) 
3 e¥OL 6255) eRLanc26c) e£Aay(56e) oV4a¥ (566) 
@ e72av(506) 
COwrcn saraenvSy ala(26e) oAL Pra (46) oBEK (236) 
1 e&f T2667) eOE (255) 00«(263) oEK(205) 
2 c€tauc25e) eELAMAT (26) o&MaT (20) eENERSV (206) 
3 efen(2eg) oGAMwA (49) oREF (20) oSGRMAT (20) 
& oS av(26c) oSUes<-{2¢8) 0X60) oXK(25¢) 
S ev(4e9) e¥K(209) oZ(4@) o2K(206) 
6 e2L4"629) ' 

CGmwtusTaRaays IFLAS (66) eo TPCINT (58) ® 
@1Swax(¢9) e1SS= 623) 0 JF (290) oJFL (206) ° 
OF 1709) oM( 200) e123) 03° (260) ® 
@yaevi2co) eJSCQtS9) oKAND (930) oKFLAS (208) ° 
@€2(209) exTve E258) ola? (50) eLAAGLE (26) ® 
®LOzFL 6200) eLTYTE 620) oMaT( 260) oPFL (20) e 
oneeT 625) etMEAY (200) 

CCH=OnsPariusys aL eALn e4LQ eo SLAF BE ° 
eHEA ee£Q edERF eCrs eOnIn 902 EPS r) 
“6a oGa% o6aR oGsiF ePEa ° 
eS waa oXT ek\ce Xe. oXREF eXTRAR eVI ry 
eYnta eVCio everr oVTiin e2t e 
eZura e209 oe Zhi eZTRAn is 

Clestys[pacevs Fear olpax olndu ° 
OClCHIav eIClANT eI olmay oJMaX 0 JNCe eKLO 2 
@KIaT eCheK oXwOu oMmak ohCw ehSOR 


INTEGER Fel oRR 
Coeeeeeoos CrrCLeTE COMPS~ENTS CF NOAMAL VEC.72 

DEL case reeere cme so< ia) 

OF. Vaz (V'.~- vegan rslOC ns 

DE, ZA 62-04-26) OEREKID/ LER (KR) 902) 
Coemeee.-o CrrCULSTE LENGTs OF acawas veEctca, 

Wee laT :DELRACO2eCELYACe2eC EL Zaee2) 
Coeeeecoos Ci CULATE DIRECTICKAL COSINES CF AORMAL VECTOR, 

ALA=lEL RAs—R ; 

BEasCELvaser 

Ganz-Cf. Zasa 
Ceeseeeeoe CHurcK COSIAES 

CALL CHECK (6-1 °AcReALNeDEne GA) 

Cat CEaLGIE-~rY?R0R) 

RETURN 

Eup 
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SUACTATIAE LR An« 
cone ooercceeececerececcosccesccoctcocccocsccccoocetoecsocesoccccacconces 
Ceeeeetces aur GEIS Yer «fay hie Rate 
Ceeesetoccersscecsesececes «OOOO COO OC CODD CoCo ee Oe BOCee COLL E Ces cos eeccels: ta 


COe~sn 7657 Mal Otese2Ch) eALa2y1559) emiahv(Secy 
1 0€632°516056) =o 60259) e€2AaSE (206) e€SS4 (250258) 
2 eG6274v(559) oFE0 (250252) eTIVLE C12) ovat (23) 
3 ev9L (289) oXL8"(295) eRALY (SSG) ovr £1 (506) 
@ e7ALV1569) ! 
Cowmce saRRAYS/ aLa(zee) oALPnma (66) eGEK(205)3 
1 e®€Ta¢66) eDEv(285) 0D (200) o€K (250) 
2 EL anr2c9) efL awat(29) °£*2T¢20) eENE ROY (2Ee) 
3 eGaxn(z99) eGl~wtiec) oREF (20) oSSa47 (29) 
@ oS'm2a7(290) oSu"Sc2 i295) oX(ee) eXK (250) 
5 svtet) eYK(256) 0240) 92" (208) 
© o7ta~-(29) 
S/T aRRAYs 1FL46 (40) ol PCINT (58) ° 
@lSwaxccg; 01802026) oF (200) eSFL (295) ° 
OSL (206) e255) owet260) oF (260) ° 
@Jeav( 279) 0 JSCR (58) eKPxD (966) oKFLAG S206) ry 
CP 1787) eaTYP (260) else (S59) oLAnGLle (26) ° 
@LOFFL (200) oLTVPE(235) o"al (200) oMFLI20) ° 
Onmat(29) OMMRAY (220) 
COMPSON P86 2657 aL eal eALea eo ALAF oRE Py 
eBER 28€2 e652F els oOmln SZ ef PS ° 
6a eGan 26aa Gar oPER ° 
*Sueaa eX] ekas5e eXCL9 eXerr oXTaAm ° 
evans eVOLo eYarr eVT2an eZ ® 
©Z%% Zin of-iF oZTRay 
Commons Pacts s Fuas olwsn - ol kOe ° 
*loeeay ,JO>aaT ela elaay oJMAX oJNce eJOLD 
K1AT oKrPax oanOe oPmax otnOn eXSQR 
INTEGEQ F eFuex 
Cowsqns82S€/ Fuar eIDEELS 


Data (COATee) : 
Coeseceoes COUKT NUMER OF RAYS FROM SOURCE ausgEeA RAGE, 
Ico w~T=1COUNTe] 
Coeccesoce IF alt RAYS FOR SOURCE Pave BEES CREATED, 6O TO AEXT SCLACE. 
IF (ICO NT-1590x (NNO) 919191619199 ; 
Ceeceesece Gr 10 SEAT SOUNCE. 
1660 NO e=he bat | 
3Co'At=] ' 
IF (WH -NSOA) 10101810116 
106] IF C1S#2x(NR5e)) 16601500185 
Coeerecooe |.CREASE RAY INTEX, 
165 lwans]raxe] 
Coeesesoce GFT Lem NUMBER, 


LT st ave (ante) / 
Coeesesoes EXT SOURCE TYPE, . 4 
ISsa=:SeaiLt) : ‘ 
Ceceeseccce E=ARCe OW SOumCE TYPE, { 
GO 101103020002800430' o1SCuK ' 


Coscsseoce E21 RAY FROM VOLUME SCUACE. 
103 CaLL vsoece 
GO TO 166 ; 
Cteseseeee Git RaY FRoa SURF ACE SURCE. = 
200 Cari SsSoRcE 
61 10 166 ) na 
Coeseeeooe GET RAY FROMIAPUT Taste 
£00 CaLL azYSOR(1COWT) 
GO TO 166 


. 
7 
af 


Ceecceeeee STORE Ray EXERGY AND UPDATE TOTAt ENERGY 4h) Ray Count, 
164 OO 166 FaleFuax ; a 
ECIMAn Fy ZESTRILT oF) 
EGaSE(FYZESCIOILT oF) OF RAC 
IF (ZEARSE (F9116701079108 
167 ERaSE(F e166 
106 Suszassuwaeee (YMAKSF) 
106 Coatisce 
IFLaG(Ivax)z} 
TRavza]sevey 
Ceercesooe YE OFAKT GREATER THan Os PAINT nay PARAvCTERS, 
IF CICPRNTSIIGeL3Q014] 
141 PANT 1OsTRkY 92: Fa (IMSK) oSETALIMAX) oGAMMALIPAX), 
OXC TAX) SY CT% 525 62 (1 MAX) 
10 FoOawAT (ime 1S CREATED REY ACoolG/4rt Akay £1206,0H BE®, £€12,4, 
OOH Gaz, E120403H RSe E12.403K Ya, El2e403n Za, El2.4) 
Wle RETuSs 
Eno 
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Fuacticn IACK(«) 
Ceatet*tesccescesseases POP SOS SORES SO SSO SOHSHSSOSOOSSOOEHSOESES CER SOVEVER 
Ceeeveteoce FUXCTIC.. lwo CETEW2A=S Jr PCLT IS INSSLE & CONE, 
Ceecestcccocccsenes FOSHS SSS SHSSSESSEHSSPSSHSSHSSHSSSSSHSSESCSSFESL SOE SSSORSCE 


Covw~on 7697 a$CO(260206) e4L=47(509) ebixs (S560) 
I eC20RLS(6050) = E 1400250) eESASE (200) eESGR(200205) 
2 e€4Rarv (500) eRES (20 ec St) oTIYLEC12) evVrat(2c) 
3 sV0OL(700) eXL2“%(266) exRav (S50) e¥#a¥(590) 
& e724v(5C0) 
COewox sZaRRAYS/ ALK (200) oALPria (40) »BE« (200) 
1 e@F7A(45) eDEV(2690) \ sOK (200) 0£% (200) 
2 c€LA4(259) oELamal {2¢) o£MaT (20) eENENGY (266) 
3 0Gan(253) oGaw~vatee) oREF (20) oe SCRMAY (20) 
@ oS:eRay 1700) eSi"SS4(260) oX(as) e 2K (200) 
5 eV¥t(40) eVK(200) - 0Z(40) eZK (206) 
6 e7L4¥(20) 
COe~gn 71 4RRAV/ IFL46 (460) eT POINT (50) r) 
354x150) e1S02(20) oJF (209) eJFL (200) e 
. #5 (200) oJ46200) ewh(209) oJP (205) e 
@JRay (266) eJSO2 (53) eK6%0(990) eo KFLAC (200) e 
@KP (700) eKkTYP (2259) elas? (Se) eLaAGLe (26) ° 
@LEEFL (250) el TYPz 12°) e“al(200) oMFL (20) e 
onus T(25) eMUMRLY (296) 
Cowon sP2=4uSs aL o4LN e4LR sALRF eBE ° 
| «8E£a eaeRe e0IS eOMt hs edZ eEPS e 
oo : o62 «GaN eGaR oGanF eFEa ® 
eSu-la |e KI eXncd eXcid oXREF oXTRAN ev! ° 
¥ @YNC« eYDLO eo Vat eoYTaan eZ ® 
ry elnta eZ0LD eZacr oZTxay 
Coe~ox.s1PsRius Fuax olwax ol0u ° 
2,3 eTozazay «1 OFFAT ek ofeasy oJMAK oJhow eJCLE 9 
@KIAT OMe th oan00 ebhwaXk eNAOW ehSOR 
INTEGER Fees AX 
a * . DEL7F=Z1-2K 64) 


Ceeeeetece IF POINT ON BAONG SICE OF YEATEK GO TO 101. 
IF (MEL ZOEN (KD LOTS 1010103 
203 CELK=KT<xa(K) 
CEL vsvleva (a) 
Ceeceetoece Puy POINT INTG EGUATICN GF CONE(FP). 
FO () Ooh C(K)) DEL ZO@2-DK (K) O (CEL KOO 2CLEL YOO?) 
Ceeee:eece IF FP LESS TnaN Os POINT YS CUTSICE CChE, 
IF 1EPYTOLeL22e102 
Ceeecece-ce PCINY IS INSISEs 
102 Incas! sy 
RETYAN 
Ceseseetece PCINT IS CUTSIDE. 
201 IAcxse} 
RET IAN . 
Eno a 
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Fuacries INcouiay 


Ceeeeetooe OHOOS OG SOS OSOSSHHOSSSSSSHSSHHHSHEEHSSSOSO SESH SSOHOSOCE Peeseesonese 
Ceeeccesoe FL QCTICW INTON LETEGMINES IF a POINT(X2eV30Z1) IS INSINE 


Cheecesooe (2° 15S. DRNILAY Ay rogkn Ry 
@ YaCo.ze}] FOR YESe 
\C@eeeoeone BAC (i20-] FOZ AD, 


CPO80 SF 00000 0000000005006 OF 0888S OBOES OOS DODO ODOOOEHODDOOOS L000 C8OOEECE 


Cowes. /69/ a8C0(230276) -@ALAA7(599) 
I 0C9547S(6059) =n 1050256) 9£3aS£ (200) 
2 eGa2av (596) eRES (290206) oTITLECI2) 
3 ev ae oXL4- (260) oRRay¢S0e@) 
&@ ePtayv(se2 

Cos~wsy parnaySy aL<{200) : oALPns (40) 
1_ sSETaceg) oCfvi2ce) 00<(290) 

2 e€L4°(250) o£L4* 51125) sEsat(20) 
3 eGax(zoc) oSa"~. (46) y @REF (20) 
@ oSPAMV(2ZES) =o Su 54206) oo A140) 

S 04s) o¥(259) eZ (66) 

6 e7La~(20 

CO--SN/ 1 aRRays 1FLAG6(668) 
@3S~2" (56) e1Scac2s) : o JF (293) 
ef. 1230) os (2ct) oSk(250) 
@.-av6209) oJSTI573 oKBX5 (900) 
on> (705) eaTv2(2e5) olaw? ($9) 
@LAFFL (200) oLTYFE (25) oMal (200) 
onwal (20) oWSMARAY (2360) 

COwsSy sP 2226s aL oALr oALR 
@BEn o8Es oScQF eDIS oD™Th 
6a oGL% oGa2a oGarF oPER 
@Sc~2BR og BI eXNCa oXCL0 oXREF 
ornca e¥Sid evar oY tian ez! 
eZmee 6ff.0 eZacer oZTRaNy 
Cowmnys[Peases Fran olmax oINOn 
elorazavy elCra%t i olaay oak 
, ORIAY na ohNGe ormAax ona 


set SYTESES F oFrax 


o6ER4Y(506) 
sEScR (200200) 
eWMeT (29) 
eVRaY(530) 


o6EK (20) 
ofA (230) 
e£EXERSY (200) 
oSudwat (20) 
oZK(2C9) 
o2K(200) 


eIPCIAT(SS) r 
oJFL (206) ° 
od? (2055) ® 
oKFLAG (260) a 
oLaKGLE (20) ® 
oMFL (20) ° 


oALRF ok. 
oDZ oEPS 
® 
sATRAN ovt ® 
° 


9 JON oJLG » 
eMSOR 


Coecoesoee- Cap. Te GISTANCE Fale FCCUS 16 POIaT, - 


Oxstfora(Ky, 
Ovsvi-vaca) >.» 
RQeS2AT (Ties2eryeugp,, 


7 


IF (QS-1 £6) 19201520199 


Coaeceereoe FIKD C5Si~ES CF VECTOR Face FOCUS TO POINT, a < 


10@ ALL 20475 
OEUstysa2 7 
CALL CHECT(4™ INCCUOALL EEL 0608) if 
DEL] OME KIC) O CELL OAL KIC) OFELOBER(K)) | 
e “ FF (AES (DENS 1 Eek DIGZ0LS2e8TI 
Ceeseetoce Fen TISTENCE Face FOCUS 10 SURF ACE» 
BO) 4072 (94 (4) ©6404) 9 /0ES 
1F (Qa-5aTV T2001 20103 | 
Ceeccceeoce POINT 5S INSICE CoNICe. ; 3 
162 Tecanz1 A 
- 6O TO 104 
Ceaeeoeooe PEIST IS OvTSiCe CONICe 
163 Yacoxse: 


106 COSTINCE : e 
Revuas 
Eno 
C47 


L—- 


ede 


say 


aie 


Fuscticnw Jacr. (a) 


Ceesante seee Pereyvyrirt SOS eees eeoteceesoeseoeese 


Ceoveccersre FUXCTIO® INCYL DETEnes: S IF a COURTEAToVE ZI) IS IKSINE 


Ceeeecs-c@ Cri INORICAL BOP.DART Au ct Ke 
COeesce see [2 CYL=e] FOR VESe 
Ceeenescoe J[*CYLs-] FOR AD. 


Coecceccccccccnsencccoosocsonncocoosevee 2000 een eSeeeeees elses conHORaOCe 
ef tev (S53) 


eatesi ta 


ovrsnti203 
eVRay (500) 


eBEK (200) 
oEK (296) 


eEREHXSY (200) 
eSGRtAT (28) 


oXK (200) 
o2K (200) 


elPorwT(Sé) 
oJFL (206) 
oJ? t20) 
eoKFLAG (200) 


eLaARGLE (20) 


COv=aGn /69/ aBCi (260286) eALRAY (506) 
1 eCOORDS (6050) 0& (403200) oF ARSE (260) 
2 eGankayv(5ee) eREDI260200) eTITLECI2) 
3 ev (209) oXL£nt206) okRay (508) 
& e7@ayv(See) 

as CO~mON sARRAYSS ALKI269) o ALPHA (0) 

1 emETA(4O) eDEWC2to) oda (206) 

2 of LA%1 200) e€L A447 (20) o&™aTc¢20) 

3 e6ax(2700) eGAreth (46) eoREF (20) 

@ eS may (2680) eSUNSOR (200) ox(4@) 
S ev(46) o¥K (260) oZ(46) 
6 o7Lan2e) i 

COned« 7] ARRAYS r IFL4G (<9) 

@3Seak(59) e3SoRn«20) oJF (200) 
OR 1(769) ePt(262) oJu(280) 
e@umay (229) 05508 (50) oKBND (909) 
KP (268) exTYP (200) eLanP (50) 
@LAFFL(IZC3) eLTYFE (26) ematc20@) 

enmaTi2c) enmaeRaAY (290) . 

Cown0n /P ARAYS/ al oALN eALa 
even. e€E2 eBERF eDIS eDa] ts 
Ga Gar e6anr Geer ePER 
eSuwvna =e KI eXx00d eATLO oKREF 
eYxoe eYOLO eYREF eYTRAK ry 2 
eZnos eZOLO oZREF eZTRan 

Cos-2ous [PARAMS Fax olmax oIsou 
CL0E22Y =e TOPRAT 3a eoyanay ovmak 

oteod 


exyat ekmax eKkNOd ommwaAx 
INTEGER FeF eax 4 

COwereecce FIND DISTAKCE FROM 
Dxxal—xK(K) 
Dysvl<-va(K) 
ZOsZ1<Z<(K) 


oMFL(20) 


eALRF 
7 4 


° 
etTasn 
° 

° 
ood 
eoNSCOR 


DOENKOS2ePy2020/D#82= (ALK (K) MOXOBEK (K) ODT eGAK (K) ©7D) Se2 


Coezeesece 15 DISTANCE LESS Than RADI US<< 
IF (DO—0K (K} O07) TOL el: 302 
Cesecetses POINT IS INSITE CILINCER. 
16] Iscy s1 
RETURG 
Cesecetece POINT 1S OUTSIDE CYLIACER, 
102 INCYL*-1 
RETURN 


Exp 


C48 


CYLINORICAL AXIS TO POINT, 


~ 


SS 


Fuacticw ISELL 4) 
Ceeeeeteooce 8S OO 8 OOS OO HS OEOSSESS OSSESES OOO SE SOSO SESSESSESSOSSOSESSESSESOREECSE 


Ceesceeeose F.n-TICN TNELL FINTS IF POINT 1S iNSICE ELLIPSOIO.W 


(0000000. 00000 00000000800 OES OOS OHSS EOESEOS CESS SOOSOESESSSSSESESEROLOSOD 


Cowon /69/ €6CO (200269) eALRayv(S09) eBERAY(S30) 

1 eCOO=rs(6eS50) 0£ (4002068) eE€GaSE i200) 0£S02 (20.290) 

2 eGenayv(see) oR2- (200200) oTITLE C12) ovrat (20) 

3 ev (200) eX 4~ (206) omRay(506) oe YRAY(500) 

& o7401(590) 

Cow“cn saanaySs AL%(200) oAlL Pra (460) oGEK(200) 

‘d omFTacee) e0E (200) e0« (200) eEK(20C)— 

2 e€(2°(209) vECL amAT (26) oE£at (20) sENMEASY (200) 

3 0&4<(2700) oGAowh (40) oREF (20) oSCum al (20) 

& eSurdayv(200) oSLMSGn (2090) 0K (60) okK (200) 

5 evi4e) evx (200) 0Z(468) oZ&(200) 

—& Fl 2120) 

Cowan /TARAAY/ IFLA6(60) eo IPOINT (S56) P) 
@15+%12(56) e1S0R(29) oF (2068) oJFL (200) ry 
PK (275) ost 200) oum(200) oJP (266) e 
@Jaav( 206) 0 JSOR (53) 988% (909) oXFLAG1200) Py 
KP (789) eaTV¥P (200) eta (Se) ol BAGLE (20) ° 
@LREFL (209) eL TYPE (20) oMaT(200) oMFL (20) ° 
enwaT (20) ommnayi2te) 

Cow sow sP2536S/ aL eal ALR - 9 ALRF oRE Py 
@AEn e6Ea edERF eDIS so ENIN oDZ sEPS os 
oGa e6t~ e6an eGakF ida Py 
eSu+sa oe XE oAnCw eXCuo oXREF oXTRAN oVI r) 
evaced eVOLD | o VREF oVTRan ry 4j r) 
elnuis e2CiD o2aEF oZTRAS 

Cowrons 1 PaRees Foak olmak oihkOw ° 
elopaay sl0rrsT ofa oTRay oJMAX o INO eJOLS 


e@nrat oKPAx enna oMmAk ote e&SCR 
IMTEGEQ FeF “an : 
Ceeccessoe CALCULATE FUNCTION FP, 
FOm) O-( (X14 (4) ) O20 LVI o YK (KR) ) OO2) JDK IK) 
On ( (ZI $7 (4) ERI) SEK(K) DOC? 
Cesceceesece [F FP LESS Tham Oe POIMT IS OUTSIOE 
IF CFP) 1020101 o LOT 
Ceeeeetsece POIAT INSIDE 
103 IMELL#I 
RETURY 
Ceeseeseee POIAT OuTSIDE 
102 IMELL =<} 
RET AN 
Ex 


c49 


Hier i palsewante dic panied en 


OR tte gerne mtn 


Ory 


eee eae 


sages 


Z 


SuUeaoQuTIVE INFLCIFL) 


Ceeceece SA SSSSOSRSSSSSSHSSSSSOHSSSSSSSSSSONS | SOOSHOSHSHASOSSSSOSSCE ICCC OES 


C@eceetoce JRFL LOG™S FCH AN INFLECTJON FCiAT. 
Ceoesecooe Pagr CF Fit IATERSECTICN Catcir.sTline 


Ceeseetoccssaceece » 2eeee POCCEOOCOREHEHSCECOCEHOCOOOSCSE COED EEDCRESE eseeete 


Coueon 769/ ABCO(Z0e20G) = e AL AAV(S66) phEaAT (585) 
“E eCCORTS(6050) 06 (400200) eEbaSe (200) eESOR (252° 75) 
2 064Rav (508) 08E0(20e200) eoTITLEC12) o¥MalT (29) 
3 ovOL (200) oXLAM(206) oXRaY (500) eVRAY (S00) 
@ e78a7v(See0) : 
COMMON JARRAYS/ ALK(200) o AL Pria (40) eBEK (200) 
a e®FTA(4Q) eDEV(260) eDa (200) oEK (200) 
2 cFLA4(268) eELAMAT (20) oEMaT(20) eERERGY (206) 
3 068K(200) oGArra (46) oREF (20) oSsarat(2c) 
@ oSuUMRay (200) oSUMSCH (209) oX(49) oAKk (209) 
5S e¥(40) oV¥K(206) 0240) oZK (206) 
© e7L4%(20) 

COv#Ons LAQRAY/ 1FL4G (40) olPO1MT (58) ° 
@1Sax (50) 01S0R(20) o JF (206) oJFL (200) ° 
OL (700) oe ( 206) o Je (299) oJ? (206) ‘ 
@JRav (200) 0 JSOR (56) ©%8N0 (900) oXFLL6(295) ° 
ox (200) eK TYP (200) eLamP (So) eLAAGLE (20) e 
@LREFL (2060) eoL TYFE (20) omaT(290) oM™FL (20) ° 
enmat(20) oNUMRAY (200) 

COrMOnsPAQLeS/ al o4LN eALQ eALRE ene 
BEY o6ER eGEnr e01S oDeIn 002 e£PS 
@6A eGAn eGaR oGanr eFEQ ry 
eSuewa ox] okade eACL0 oRREF eXTRAN evi 
eYnOd eYOLD eoYuHEF eVTRan ry 4 | ry 
©Z4O« e200 oZREF oZTetAm a 

CowvensiFsrsrs FuaAXK elwax eo IkGe ° 
e10eaaY 6 10?RNT eIR el&aY o UMA tt +JOLC 
@KIAT oKMAX oX.0e oMrax oMNOE ehSOR 


INTEGER FoF =ax 
COMMON JHE LCOM/ Ae KBAG2B 0B? eBIP CEP oCAPeDELTMeDELNeDELYOCELZ o7-~17 0 
OPS1T o RE Ve Ree ST rg Tite TH] 9p SCANS eo hAKSy SToty SNTH RES THe HSNTH ES” 46 
PCOSTANS FHZ oCTHZ p 
Ifuize 
Ceaccesece JF FUACTIOW LCOXS LIKE A P£k Ae 6O TO 162. 
1F ((G=R}) 7 (CHP OB1P) © Tre TH 1) 2-625) 10101020102 
Ceeccetece 1F BOOT IS SQACKETEC. F1NG AOOT. 
102 1 ¢(8-a8£9)6(81-A8A8) ) 10301030110 
103 CaLL Fixca 
6O TO Je | 
Coaeecsecee SAVE STARTING VALUES. 
163 TH221H 
6228 
S2rpsEP 
TH3sTo} | 
63-41 
83°=<R21P 
Ceweeecece CALCULATE KEW ANGLE 
1046 Teen ( TrieleH1) /200 
Ceeeseccee CALCULATE MEW B AND BP. 
Cart eBcarc 
CHeeeeseoce JF SIGN REVERSAL» CALL REVS TwICe. 
1F (91 PeBP) 10501850106 
105 CALL REVS 
IfL#1 
TH) Te 
6128 Ps 


B61 e=ae 
THe Ta? 
b=q? 
82s22P 
Catt REVS 
GO 70 Ile 
Coececetece IF FUNCTION LCOKS LIKE A PARAACLA 
106 IF ( (R219 /( (820410) ®{THoTH]) 1-025) lose lOre tes ia 
Ceeeescese TF ROOT 1S BAACKETED:e Fl) ROOT 
107 IF ( (Ro £92R) © (B81 —482A) ) 10601680109 
gee Cait Fyxca 
Coaecceece SET LIMITS FOR OTHER PART OF FUXCTION. 
109 T3104 tien 
8)=9 
BlD=t>a 
Tea” 2 
B==- 
82-:2> 
Ceesccecce IF FUNCTION LOOKS LIKE PaRA@cLa 
IF ((8-021/ (127 0A29) 0 (Tea Tm) 9025) loved Deedee 
Céeecescoe JF SCTT 15 ARELKETED. FIND ROOT. 
VIL UF (8-225) 0 (51 -48229) 11201120920 
W2 Cary Ty: 7a 
110 BETUaS Po 
Eun 


C51 


Wdvenbss ash nss 5. 


Fuacticn IWs€L («) 
(occcccoccccessccccscccooccsooooooasosesoooss® epecgessococcooccccorsese® 
Ceneccccee Je _ DETEMWIAES IF PCIAT KleVIeZ! JS IASISE rziis. 


Ceseeossecece: 3 20000800008 60000 SOS SOOOCOSSOS J SECCESESSCOSORO OOS SOOSEOOF ee 


Cowres. «(34 asco 220206) eoaraayviSte) eGtaltviScd) 
1 eCOCa: S' 1056) oF (400255) o&saSE (263) of SC" (250200) 
2 eGAaRev (36) e2E9(26 e265) oTITLCECI2) evralt(2¢) 

3 eva (200) oXL£e(205) oxaav(563) eVaav(So0) 
@ oZhsv (S00) 

Cowon saRRAYS/ ALa(200) oAL Pn (40) eBEx(2¢¢) 

1 oMETA(4E) eCEvi26c) 00 (200) oEX (268) 

2 c€L ANC 206) eELa~27(20), o€"aT¢29) oEMERSY (226) 

3 eGax(760) oGarra(ec) oREF (26) oSGamaT( 26) 
@ oS PR: V1200) oSu"SOr 1250) oX(40) 0X8 (208) 

5. evt4e) ova (200) oZ(46} oZA(26u3 
@ o7iani26) 

COw~Ons LARRLY/ IFLAS (46) oIPOINT (S92 ° 
@1Ssax (£6) o1SCRi20) oJF (263) oJFL(206) ° 
#1206) oP259) . 05m (200) oP (200) ® 
omav (260) 03502156) 0X50 15060) oBFL 61255) ° 
oP 1206) exTY? (200) ela“ (59) el a®GLe£ (25) ° 
PLREFLI26€) eoLTvPe ize) oPat(230) o™FL (26) e 
oomraTict) onurmiaY (200) 

Cows. /72R64S/ aL oALh o4L8 oALAF ont ° 
OBE eaten evtRF eDIS o0r In eDZ oEPS e 
ea 06a e622 oGaRF 6, PER ° 
eSueaa =o KI okhow oXCL0 olREF ekTaar ev! ; 
evxcy ov eo YRee oYTRan ry 4 r) 
eZucw ) olaté oZTALE 

Cownnns [PLGaws FraAL olwan olkte ® 
elcraat eIC?=~.T ok oIkAaY oak ssvcu eJOC 
exiaT ° oKnOw ommak one e&SOR 


IntEGER FoF wax 
OELErKI-" ««) 
OELY=VI-vA(C) ‘ 
Ceoecesseoce CALCULATE DISTANCE FROM AXIS TO PCINT. 


Rar EL KeeFesz_veoe2 
AIns(REK (KICKER) DCC? 
aasTotnretts coated o02 
Cocccesecss IF PCINT SOT SETCEES IANER ANG OUTER RADILSs GO TC 182, 


IF (FRR) LOZ 1 O30 TEI 
103 IF taere TH1S402O40TI92 


yes CCATIar 
Ceeesevoce LCCATE KEASEST POINT CH PELiIca, wIRE., 
Cait PCIAT (SALT SEPT oSAPT OXI oVI ZI) 
Ceceepecce CALCULATE DISTANCE StUARED GETEEEA TaO PCIATS~ 
DSaALPloe2ecEsTee2ee PTS“? 
Cececeecce IF OISTANCE $2~2ED LESS Tsk TUBE RADIUS SOUA5ED» SO TC 10I. 
1F (05-24 (03 002) TET oTO3 9162 
Ceceesooee PCINT 1S IASICE 
101 IsM-€Ls 
RE TUAX 


Ceescesooe PCINT IS GUTSIDE 
302 ImELs~1 
RETURN 


Exp 7 


FUSCTICN Innve cay 
Ceeecetor > 000000080 000008000000000000000 0050008 008000b8C0e Feeegoscesete 


COoccecteoe Fine TICN Inedy>? FIAOS IF PCINT 1S INSICE HYPERDDLL ID. 
Coscccecorcccconeccesccccccccccooccsetoccccceccccesoeecooeccccocoecccete 

Comess 7497 48C0(260200) eALQ4V(S59) oBERAT (S00) 

1 e€CIATS (6050) §=00 E1400200) oESaSE (209) eESOR (22,290) 
2 eGa5ayv(S20) oRED (260209) oTITLE C12) o¥ral (20) 

3 evfL(200) oXLa%(200) oXRAY(SC60) oVRAY(S50) 
@ eZRAV(5901 

COw“0N saRRaYSs ALK(206) oAL Pca (40) oBEK (204) 

1 om©Tacéep eDEVI(209) 0041200) oEK(200) 

2 oft 4"(250) oEL AMAT (20) o&*aT(20) ENERGY (200) 
3 064<1202) oGatHa (4c)  OREF (20) oSCanat(20) 
@ 0S~287(260) oSUMSGR (290) oX(460) o XK (206) 

5S ev(40) eo ¥K (200) 0Z(40) 02K (200) 

6 e2.2"(20) ; 

Cee-c ../Jaaaays IFL46 (460) 0 IPOTNT (S50) ° 
@lSwexn¢to) e3S0R(20) oJF (200) oSFL (200) e 
AN (200) oI4(250) 0JN (290) oJP (200) e 
Rav (200) o JS02 (50) oKBN0 (990) oKFLAGI200) ° 
@xP(200) = «= eo KT YP(200) La (S90) eL ANGLE (23) ° 
@CLOzFLizecr eL TYPE (20) eal (200) oMFL (20) ° 
oneaT( 2c) omUurRAY (200) 

COwvOa sr eemss AL ALN e4La eo ALAF oRE ° 
BEA eAtR eSfaF oDIS oOMlns 00Z eEPS r) 
Ga oGA% oGak oGaar | oPER ° 
eSucza eX! oARld oACLO oXREF oXTRAN eV! ° 
ernce eYCLO 0 YAEF oYTRan oz ° 
ezNne ry sere oZREr oZTAAN 

ComvonsyPrsavs Fouax olmax o know ° 
@10P2ay ef SCPENT ota oTaay osax oJNOM eJOLD 
exIat Shan oKA0w oMMAX onnOw oNSOR 


INTEGES © Fess 
Cecoeeeroce Ce _CULATE FURCTION Fo, 
OELZ22 7-72 (n) 
IF (CEL Z@GAK (K)9 10201020153 
103 CrLaxeoxq(K) 
OEL vavJovn(K) 
Fae) Om (DELXOOZ DEL VO*2) /0K (K) 0 (CEL ZOEK (K) 1002) / (EX UK) 00) 
Coeepoe-oe IF FP LESS THaN Ov POINT IS OUTSICE. 
IF (©P910201010101 
Ceececcetece PCINT InSIDE 
101 IsnyvPa) 
RETURN 
Ceeeereeu0 PEIAT CUISIDE 
102 IkevPse}] 
RETURN * 
Eao 


PeislavPaiuee 


FUACTICN Ia?2Q(K) 


( 


coe ee een neces ecereccccncccsccecscsceetcensccccssccccescceccccicccsceecs 
Fu.CTiis 7.P2R FIADS IF PCIAT TS INSICE PePascrcsc, 


Ceasectcooccoeeen: «ane OOOO OOO On eececeeseseacescoceeconccecelseseccecccece 


eBERLY (SC) 
vESSA (250290) 
oveal (20) 
oYnAYV(Sos5) 


o3Ea(2¢09) 
oER (256) 
eENE~Sy¥ (209) 
oSCR*4T (20) 
04K (200) 

oZK (200) 


oe IPOlmT (50) 


‘oJSFL (200) 


ComeCn 7/657 4£8C0(260290) oSL2£7(505) 
‘ I eC9O2ATS (6.50) oF (4002CC) 0E64SE (200) 
if 2 e€akAy (500) oGED( 260290) eTITLE (32) 
3 evVOL (765) , oXLan(200) oXRay (562) 
/ . @ o7RAV(5.9) 
! peek: COMMON /2RAAYS/ §=ALK(200) eALP~A (49) 
1 oRETA(40) oDEV(200) 00K (299) 
2 s€L4™ (260) oEL AMAT (20) oEMaT(29) 
3 064K (2¢0) oGAmta (00) oREF (20) 
7 & oStmRay (200) oSUMSSR (260) oX(45) 
a 5S ev(4o) e¥K(2CO) oZ (40) 
By 6 07 40(20) 
) COwsONs 1 aRaays IFLAG (469) 
ms @1Smax(sg), 01S0R(20) oJF (200) 
@A1(200) oJ(200) 0 Je (252) 
@JR AY (260) | 0 #$0a(S50) oKGax0 (969) 
F @KP (790) oKTYP (2-5) ol artery) 
*LREFL (200) oL TYPE (29) oMaT (256) 
@nR4AT (20) ORUMRAY (200) 
Cowson /PaRnauss at o4LN eALQ 
eSEn o8Ea eSEar 0318S oOMIn 
M @GA oGAN 0GAa oGanF oPER 
@SuUrnA og XI oXR0e X20 oXREF 
; @YNne eYOLD oYREF oVTaay ry 4 | 
7: ‘ e2now eZOLN oZREF oZTaiay 
CO¥ 40K 719 4c48K7 Fmax ofrax oIhCu 
sw @C1O0taay ,ICPRAT oR .  oFRay o JMAX 
oNOw 


eKInT oKM ax oXkOe OMMAX . 
INTEGER F oF Max : F 
Ceeceseree CALCULATE FUACTIGA FP. 


FP28 96K (K)©(7162K(K) Do(XTORK(K) )@02e¢ 
Cooceoecce IF FP LESS Tha Oo POINT 15 OuTsic 


2 

Coeeccecce PCInT 1S INSIDE, 

102 InPaks] 

RETURN 

: “Ceeawecece Point 1S OUTSIDE. 
i Ie] Inpaks.] 
RETURN 7 
Eng 


C54 


5 ited 


° 
° 
oJP(200) r) 
oXFLAG (206) ° 

oLANSLE (20) ° 

oMFL (20) ° 

oALRF oBE ) 
002 oEPS og 
° € 4 
oXTaan ov! ° 
r 


5 , 
pees eJOLD 6 
° 


YVICYK (x) poe? 
E 


FUACTICN IAPLEN AK) 
Coccccccccccccesecccccsccccesccccccenccccccccccccecccccccccoccccoccoccces 
Ceoceecsesse FLACTICK ISPLAN DETENMIAES IF & PCINTIAI oVEeZI) 15 On 
Ceeeecescos Tuf ORIGIN SIDE OF PLEnt BOUADANY RUNGER x, 
eee [AP_Bhze) FOR YESs 


ce ee [8.7L skze] FOR NO, 
CPeeeetcoocense (OOOO RCO OCC C ROC CEE OEO CDEC CCONOe COED EEEOCCEECCCCCCOESeS 
Co~~o% 7697 adC0(2002¢c0) #ALRayv(S500) eBERAY (500) 
J 0€90275516050) =n E 1400200) o€9aSE (230) sESOR (200290) é 
2 eGhRay (504) eRED (26.2003 eoTITLECR2) o¥™al (20) 
3 ev (2G0) oXLan(230} oXRayv(500) eYRAY( 500) 
@ e7827(506) 
COPHcn saRmarss ALK(209) oAL Pg (40) oBEx (2900) 
I e@F T2149) eDEV«I2sCc) eDK (200) eEK (208) 
2 c€L4%1200) eEL AMAT (205 o€£aT (29) sERERSY (206) 
3 0642200) 0642160) oREF (20) oSGInal (29) 
@ eSu"27(200) oSuM"S0e (200) oX(49) oXh (200) 
S eves) eV (200) oZ(46@) 02% (200) 
6 e7t £120) 

Cowman s 1 aQRays TFLAG( ae) oe TPOINI (S50) ° 
elSvanicgy ersoai2c) oJF (200) oe JFL (200) e 
OA 70) ot(200) oJu(200) oJP (200) ° 
OS-20 1250) 0JS02(50) eKBND (900) oXFL AG (200) e 
OK? 57559 eKTYP(2C9) olen (50) eLAAGLE (26) ° 
W557. 1209) oC TYPE (20) oMAT (200) oMFL. (203 ® 
onus (20) erUmRav (260) 

CCwrONsPARE4S/ AL oALN e4La eALAF one e 
BER 96E2 oUERF eDIS oOnza oDZ eFPS e 
e6a oGaN 0G4aR oGanF oPER ® 
eS yung oXE oXnOw eXOULD - oRREF eXTRAN evi e 
evnce eYOLD eo VREF vVTRAN ry 4} e 
oZ%ow eZCLo oZREF oZTRAN 

COw0Ns PAR aes Fax odmax oInOw e- 

@ICsaay yg fCPANT eI o3Ray oJMak oJnQe 2eJOLD 5 
@KIAT oMhax ohKOe oPmAK ohmeOu e™SOR 


INTEGER € oFeax 
Ceesceeesee FIRM DISTANCE FROM POINT TO PLane. 
_ ODBALK(K) OX] OBEK(K) OVS 0GAK (K) OZ eDK (Ky 
IF (MOD Letele2cree 
Ceeececege POAT NOT OW CRIGIN SIDE. 
100 INPL anne} 
RE TUAN 
Ceeesetece POIKT OM ORIGIN SIDE. 
VOT Ine, ans] 
RETYUAK 
Ceesccvese POrmT ON PLANE. 
162 CALL EAPRNT(6HTNPLAN) 
RETURN 
Eug 


CSS 


alton 


FuacTic¢sx IAS€s¢s) 


COCO 0006 60066 O8 OF OEDHEESSAOECOEEEO OO OEE FESS OHOEOHEEESEEPECOSEEEEOOOOCE 


Ceccsceoes FLUCti im INSEG CETEMSINES IF POINTER eVEeZI) IS INSITE 


Coesceesoce SF GHEtT MUR==Q Je 
Ceeeeseseee IeSeGs] FOX VES. 
Cesesetece 1255625 FOH °.7, 


COCR 000000800008 SOC OOS SEOEESEOEEPOOEDEO CEASE OH OES ESEETOOEECES EOC EESLEED 


CO-~. /65/ ASC0i290255) e&LR2v(569) ent RAY(563) 
I eCMO6TS (6058) «=n E (4902951 sEGASE (290) oESC4(200200) 
2 e&ARAy (SEO) eRED1200230) e TITLE (12) ovat (20) 

3 ov (2085 oxLati2¢0) eXHav(SEO) eVeav(S9e) 
@ e7Rav50e) 

Copecn sAS=8VS/Z BLK(200) eALPrA (40) oBE« (200) 

1 eRETA(4E) eDEV(265) e04(296) eEK (200) 

« oF LAN(200) eELameT (20) eat (20) eEXESSv (200) 
3 eGAK(200) eGares (46) eREF (20) eSCrtat (2c) 
& oS'Pa4V (262) eoSlU"SGa (286) oK(60) 0%" (260) 

5S evi4e) 0 VK&(209) 0Z(49) o2& (260) 

@ e7L44(20) 

COwwox 7 ] aRRAY/ IFL26(408) oe IPC INT (59) ° 
®I1S"axktSe) © 0 1S62(28) owF (206) oSFL (200) e 
(200) ot(286) oe (206) 0 JP (200) ° 
@F2V(200) 0 JS0R (5¢) eB WT (909) oMFL 466250) ry 
enP (790) eKTYP (296) eLarr (59) eLARGLE (26) e 

, @LREFL(200) eL TYPE (20) oat (266) oMFL (20) ® 
onwaT(2c) _ eUMRAY (208) 

COwenQe 7? 22S aL oAgM 248 eALRF AL 3 ° 
Bia eGBER eeRF e0IS eOMl& 002 eEPS e 
6a oGas eGar eGanrF ePER e 
eSurana eX] eXa0u eAC.D elkEr eX TRAM evi e 
evagu eVl.d eVREF eVTRay ry 4 | e 
e2uou PY Samery) oZReF oZTtang 

COC ys P As +7 FrAK olPce eIkOu ® 

- @]epaay gliarxT = 91k olGaV oJvax oJNOd eJOLD ¢ 
oxy: eth Ad eKnGu ohm ak ehhou eNSOR 


TacteGEe FoFMAR 


Ceoseoeeoes GET IRDEX FOR BCUADARY ARRAY FOR SEGHWENT J. 


Se=aP (J) 


Ceesssceos GET IRDEX K FOR ORE BCUNDLAY OF SEGYEKT J. 


106 Kei 2AS (KEN (me)? 
Ceeenetoce SET GECWETRY TYPLe 
% KTe<TyP (a) 
KS I CXcaFAD (bh) 
Cessectece granCn OV’ GEOWETAY TYFE. 


GO TOC1O6 02900395 0606 e550 ef C60 760 oe BGs 0 SOG) KT 


Ceoaevcereeece [Ff PCI*T Of DESIALO SILE CF 
100 1F (ISPLAN(K) ORS TEN) TOZ0162059) 


, Ceaseetsee IF POINT Ow DESIRED SICE OF 


200 IF CIWCYL (K) OKSIGA) 10201620161 
Ceesectooe IF POIKT CX DESIFED SIDE 
300 IF C1eCCm(K) ORS IGRIIO2Z01020161 
Ceececeoos IF POIRT OV DESIRED SIDE 
400 IF CIASHPH(K) ORS IGW) 10201620181 
Cesessesese IF POIST O% DESIXED SILE 
See IF c1wCa (K) ©<S1G4) 16201620161 
Cosceetece [F POINT O% DESIRED SIVE 
600 IF CIMPAR(K) KS TEX? 1062016020161 
Ceaceetsce IF POINT Ow DESIRED SIDE 
700 IF ( TeetvP (CK) OKSISR) 10201020101 
Ceecceetece IF POINT ON DESIRED SICE 
see TF CIMELL (K) OMS 164) 16201620161 
{ 


ng 
ae IF POIRT o« OES 2 SIDE CF 
Pits @ IF ey (€) OKSIGaD IG ol&2elul 
eestcce 1 AT REKT BOUAOANY 
101 Reaoy ny 


IF : 
F (SAD (41110651035 168 


% 
103 Imsecky | 1S INSIOE Seorear, 
é 6 10 Vie 
eeese@oee EGInt 1 
202 lasecae S CUTSIDE SEGrEat. 
1l@ Costiace 
| RETURN 
Eup 
Gq 
\, C56 


RRIRRA 


Prartese GO TO 161 
CYLINDER. GO TO 101 
COMICs GO TC 161 
SPHEQEs GC TO 101s 
COME ¢ GO TO lél. 
PARABCLOIDs SO TO 181. 
MYPENBOLOIOe GO 10 1¢1. 
ELLIPSOIDe GO TO lols 


MELIXe GO To 101. 


“ 


wr 


FUACTICN IASPH(A) 


ce sesesercececscececesessecscoccccesescocceeseeesenccccesecescecscseeeee 


Coccoveeeoe FUACTION 15SPm FINOS IF POIAT IS INSIDE SPHERE, 


Coccccccccccccccccccccccccecceececcnntccccccecccoceocccetetececccccccces 


Comeda 7697 A£6C0 (200200) sALQAy (500) 
1 0C7020S16050) §=9£ (400200) sEGaSE (203) 
2 eGeRAy (500) eRED (205200) oTITLE C22) 
3 eV¥NL 1200) eXLAM(200) © eXRay¢(5605 
4 7241 (500) 

Co“vCs saGaaySs ai«i200) eo AL P48 (469) 
1 eRETA(460) eDEvi209) 9DKx (200) 

2 s&L4™(200) sELamal (20) oEmal (20) 
3 064 (200) eGAv~vA (40) eREF (20) 
@ oSMHAY(Z00) «= e SLU*SOGRI 256) 0X42) 

S e¥<40) e¥K 6200) 0Z (40) 

6 e7L4m(25) 

Commons Tass ays TFLAG (40) 
154x159) 61808029) o JF (200) 

JL (209) aJ4(200) oJ (206) 
@JQav12c0) -@JSCRI5S) 9Ka~D (909) 
@KP (200) oXTYP (299) olan? (59) 
@LREFL (200) eLTYPE (20) 04aT(200) 
MAT (20) oNUMRAY (200) 

COmsCysPeaauss — aL o ALN eALR 
eBon eAzZR = GERF s01S oOMIn 
6a oGAn eGaRr o Sarr ePER 
*Suvaa ox! eXN0e eXCi0 oXREF 
evans eVCre o VREF eVTRAN eZ] 
oznoe eZCLO oZREF eZTRAN 

Covvcns{FLesus F max olwax olnOw 
elcsaay .fCPRAT =o IR ol ty o JHAX 
KYAT oX¥AX oe KNOw ot AX oNNDe 


InwtFGER Fo Fuan 


* o SFL (200) 


sBEHAY (500) 
vESCR (200200) 
svat (20) 

eo VRaY (560) 


oBEK (200) 
sEK(200) - 
sERERGY (200) 
oSuRMal (26) 
oXK (206) 

02K (260) 


sIPCINT (56) 


i 
® 
9J? (200), A 
oKFLAG(200) yg 
eLAAGLE (20) s 
oMFL (20) ‘ 


eALAF AE 


yy 4 EPS r 


i 
(eXTRANg¥I , 
8 


r 
oJkOu eJOLD 9 
oNSOR 


Ceeeeeoee CaLCLLATE OISTASCE To CENTER CF SPnzzE. 
SSs(XT“KK IK) 0OZe (VI “VK (K) ) 0020 (ZZ m7ZK (Kd) 002 
Cecoeceeeee IF GISTAY.CE GREATER TRAN RADIUSe POINT 15 OUTSIC=, 


IF (SSoOK (KI 002) LAT eTOL e162 
Ceeeeeteooe POIAT IS INSIDE} 
1601 Ingonsy 
RETUAN 
Ceeeeceooe POINT IS CUTISITE. 
102 INS©~s-3 
RETURN 
Eno 


SUAROUTINE INTCN (Ke DIT) 
COS SOOO 0806 C0 OOO O8OS OOOOH OOOES SOS OOS SSOOS ODS OSOSOOSSSOOOOOSOOHOESOOOOSOSS 
Ceescessoes JATCN FINDS INTERSECTION CF Ray nlTe A CONE. SEE eRITEUP 


Ceeecetees FCA EQUATIONS USED. 
(009000000 0000000002 000000000280 00000 Se ncncneessessensesseesesoecqsooree 


Cowvon 7697 46C0120+200) eaLeay (500) eBERAY (500) 
I eC2ORTS(6e50) 2E (400200) eEcaS= (200) £508 (260200) 
2 eGanay (S500) eREO (200200) eTITLECI2) evmal (20) 

3 eV9L (200) eXLAnk(200) eXRav(500) eYHAY (500) 
4 e7RAY1(S500) 

Cowto~ saRRAYS/ ALK(200) eo AL Pra (40) eBEK (200) 

I eAETA(460) eDEV(200) e0K (260) 0EK (200) 

2 e€LAm( 200) oELA™AT (201 eEmat (20) oENERSY (200) 

3 eG4K (200) eGArHaled) eREF (20) eSCaAPAat (20) 
& eSU”Ray (200) eSU“SOR (200) eX60) AK (200) 
Sev(4o) «ss e YK (200) 0Z(40) eZK(200) 
© eZLais(2e) / / 

CoewOus 1 aRRAYs i IFLAE (40) e IPOINT (50) e 
@1S*+<156) eISOR(20) oe JF (200) o JFL (200) ® 
A (208) eJ46200) 93h (200) eJP (200) ° 
@RMayvi2o0e) e 50) #KBND (960) oKFLAG(200) e 
KP (290) eXTYP (260) eLaeP (50) eLAAGLE (20) e 
@LREFL (200) eL TYPE (20! ePat (200) oMFL (20) ° 
nmat (20) eMUMRAY (260) : 

CO=~0%/PaRaMss aL ALN eALR oe ALRE one ° 
@BEN eBER eSERF 001s Pd 002 | o€PS ° 
*G6A . eGan oGaRr eGanFr ePER e 
eSue2a = og XY eAnOw eACLO eXZREF eXT&&aN eV! ° 
OYNC ds eVOLD eYREF eYTRAN eZ! Ca 
eZn0u eZOLD oZREF eZTRAN - 

Cosa 7 1PaReMs Fmak olwax eo INOw e 
elorvay ,ICORAT ea eyRay eax eJrCe eJOLD 6 
@KIKT eKPak oKieOw oMMAX eNnOw ehSOR 


INTEGER F oF “ax 
Ceceneeoce Carcu_aTE PARAMETERS AEEDED FOR SISTERCE EQUATSCK, 
OEL Rex TouK(K) . 
DEL YsVIovK(K) 
DEL 7=ZFoZA Ky 
Care, Zee2-0K (K) © (DEL X©P2CDEL VOC 2epE| 2002) 
B=) .0-DK (K) ) ©GACDELZ “DK (K) © (CAL °CELXeBE*DELY) 
AzG Ae 2-DK(K) 
$Q2R002-4°C 
IF (SO) 10101020102 
1e2 RonTs$¢RT(S0)} 
Ceeceeteoe CALCULATE SOTH ROOTS AND CHOOSE THE SMALLEST POSITIVE 
Ceeeeseeos ACOT THAT IS 1% THE DESINED WALF CF CONE 
On} T15¢-BeROOT) 7A 
O01) T2x (-B-ROOTI 7A 
IF (Net ET1) 10101032103 
103 EAL x(ZJeDHlT)] OSA ZKIK) DOEK(K) : | 
Empat ZePHl T2¢GA~ZK(K)) OER (K) f 
IF (OH1T2) 10501050106 
104 IF (FT) 1050105,106 
105 IF (EPLI10101010107 
106 On] TsDe1T2 
GO TO 108 
107 OnyTs0r]T) 
Go TO 108 
101 Ow{T=) £021 
106 Cali OEBUG(OnINTCm ) 
RETURN 2 


Eun 


C58 


SUARQCUTIAE INTCOW (Ke DHIT) 


Coccccscccccccccnccccncsc: cocnesssescecesoee exssscenesesccooooscsonsose 


Ceesccecee SURROUTIKE INTCON CALCULATES Tet 


OISTAMCE FROM POINT 


Ceeeccecece KieVieZI TO Trt POINT emEME Te M&Y IATERSECIS COAIC 
Cee © ROUNDA2Y a. SEE wat TEue FOr EAaCT EQUATIONS USED. 
(0000000000000 00 0000000000 0000 000000880000008000SE8EES 2800S E5008 00000588 


. Cnw2Q~ /69/ ABCO(2092A0) eALGay (500) sHERAT (500) 

1 0€9027S(6050) = s E (400200) o€9aSE (200) £504 (20.200) 
2 eGaRary (S506) eRED (266200) oTITLE Cle) eovmal (20) 

3 ev (700) oSL 441200) oxRay (500) oe VYRaY (500) 

& e72471500) 3 

Co~“ou sARRZYS/ = ALK(200) eALPA(4E) eBEx (200) 

I se@ETA(49} sDEV (200) 9Dn1200) o€K (200) 

2 oF Lam 200) sELAMAT (20) Ea (20) sEKERGY (200) 

3 264*xt2700) GAs (46) eREF (20) o SGAMRAT (20) 

& oS'™Aay (200) 9SU"SOR (260) oX(4e) oAK (260) 
= Phaey )) oY (200) oZ(40) 02% (200) 
o7L 4" (20) 

gone Re IFLAG (40) el POINT (50) ° 
e1svax se) 9150R (20) oJF (200) oJFL (200) ° 
5 (2008) oJ™ (200) oJu(208) oJP (200) ; 
®Mavi230) 9J50R (50) 9KEnD (990) oXFLAG (200) r) 
onP (700) exTYP (200) eta (50) eLAhGLE (20) ° 
@LREFL (200) eL TYPE (20) oMaTi2ee) o™FL (20) ° 
omaaTizer stemMAY (200) 

CO~~<0a /PARRMS/ aL oALy sALR oALRF oWE ° 
@BES | @@&a oGERF 9018 oDAin 002 s€PS ° 
@Ga 1 eGAn eGaR oGaRF sPEa ® v 
@SurRa oXI oXhOa 9xou0 oXREF oXTRAN ev! ° 
eC¥nnd eYOLO o YREF eVTRAN ry 4 | ° 
eZuew eZOLD 0 ZREF e2TRan . 

Courdas[PaRass FOAR olmak 9 160u ® 
elopaay eICoaNT =o AA elaat oJMAX 0 JNOM sJOLD 9 
exyaT oKMAK oKnOe ooAd obetetea oMSOR 


INTEGER FeFran 


Ceeeeeeoce CALCULATE COEFFICIENTS FOR \DISTANCE EQUATION. 


Oxek{[orK (a) 

OvsviovacK) 

IF 60350 (6A) ©°2-] 00) -1 Eb) qeeeneeeset 
101 COSTHsSCAT (1 .6-GAee2) 

ALLSAL/COSTM 

GEL =BE/COSTH 

CALL CHECK (4+ Iu TCON CALL otkl 0000) 

E€sNZeaLt oOYe-E=L 

FF eOKOALK (XK) oDY°BERK (K) 

GGsALK (K) ALL MER (0) OBEL 

BSeNhx0ee2e0yee2 

PP 21 o0-(EK(K} 902) © (66992) 

IF (EAS (PP) —1 E20) 10601549102 
102 beesFEWEK (K) 0820660 (FF eDK IK) ) 

O0sAS—E am (K) 929 (FF DAK (K) OO? 

MS$EGs-01002-PPe00 

IF (4SPC) 1942103103 


Cecceseese FIM BOTs ROOTS OF Tne EQUATION ahD CHOOSE uy SHORTEST 


Ceeesetese POSITIVE ROOT. 
103 ROOTeSCRT (HSPQ) 
Det gy TE (-teOOT) / (PPOCCST ot) 
007 T22 (ot oNOOT) / (PPOCCSTH) 
_ ONT Teas {NS (CHIT Ne 112) 
IF (00117) 10501050107 
105 Ooo} Taawaxy] (Oni T) ed T2) 


IF ("43 1920401640107 
100 Ony Tel E21 
107 COATIAAE 
CALL DERUGIEHINTCON) 
RE Titty 
Eno 


C59 


SUMROUTIBE INTCVL (Ke DHIT) 


Ceceee POOOSSCOSE SOS O8S8OSO8S0 006 00606002055 500 0060000006090 0000000000088 


Ceocenectece SLARQUTINE INTCYL CALCULATES Tok CISTANCE FROW PolaT 
Cosccoeoes ApeVIeZI IW Tee OIRECTICN BLeceGA TO TE POIN? oF ERE JHE 
Cesceeesees Ray IRTEASECTS CYLIADESCA, bOuNvaGY Ke SEE eRTifue FOR 
Ceaseuseee EXACT ECUATIOKS USEC. 


Céesseesccoesecesooaons Seeccesccesececccces SONSHSSEC HSS SOKSESCoSeessee®te 


Cor~0ON /59/ AB8CO(280266) = pa Rav«See@) eBERKT (500) 

1 eCOORDS(6050) 0& (400200) ‘eEGASE (206) ef (20.200) 
2 eGaRAy (S06) e8€0 (20.280) eTITLE«l2) e T¢2e0) 

3 ev0L(200) oKLAn(206) oKRay (S50) eVRAY (S00) 

4 eF72avi5¢ce) 

Coesty sARRAYS/ ALa(200) eaLpna(4e) ° (203) 

1 eAF Taree) eDEVI200) e0K (290) eEK i200) 

2 c€L a1 206) oELANAT (29) e€r al (26) eEMERGY (200) 
3 eG4x(200) eGarrun (48) eREF (20) ry AT (20) 

© eS'mRav(200) «= o SU*SON (208) oX(60) _ ekk(200) 

5 evtse) eo¥K (200) oZ(4e) 02K (206) 

& eZLAmI20) 

CO0e~)us 1] ARRAYS IFLAG (46) eIPOINT (50) e 
@1Srax(Se) e¥S0R 128) edF (200) eJFL( 260); e 
5. (200) oJ%1200) oJh (200) oP (230) | e 
@snav (200) eJSOR 158) eKHAD (900) oXFLAG (200) e 
@KP (260) eoXTYP (208! ela (S55) eL AAGLE (28) e 
@LREFL (200) eoL TYPE (26) omat (200) oMFL (20) | e 
enonaT (20) eMewraaYi20e) a. 
Coon sPagaNs/ aL oALee oALR eALF gE 
@B8Ew eBER eBERF e01S oD IN eDZ eEPS 
SGA Gan eGaR eGanFt ePER ° . 
eSuuaa =o KI oXhOw eACLO » @AREF oX TRAN eYI 
eYnge eYOD eYREF eVraan ry 4 | e 

eZac4 eZOLD o ZREF oZTeotanm 

Cow <n TP aRnsns Fax oTmax oInow e 

C1OPAAY e1OPRNT 869 If eTaay oJMAk e eJOLO 
@KIAY oXPAax okuOw ri ¢ eed e 


IMTEGER Fef "2k 
Cevevetcce CaLCULSLTE COEFFICIENTS FOR DISTAWCE EQUATION. 
KXs6L AL a (K) oBE°BEK(K) (GLOEAK (a) 
Ovssx] .g-x3Fea2 
1F (08S (OWXS) -} .E-6) LoTeT@1e1e2 
102 DEL xsxzougi<) 
OELYsYJoYa(«) 
OELZ2Z1-ZK(K) 
VYs al PEL KeuE “DEL YeGA°DELZ 
ZZsAL<(K) OPEL ABER (K) DEL Y eGAK (K) ODELZ 
BSI2NE| LOO 2eDELYS°® 2eDELZ¢02 
SOT#AM= (YYZZEXL) CO2— (CHKS) © (ASO°7ZZ6°2—DK (K) eee} 
IF (SQTER) 10101030163 
Ceeseteoce FAD BOTH ROOTS OF Tee EQUATION AAD CrOOSE Te SHCRTEST 
Ceceseocce POSITIVE M00T. 
163 G0o0T=Scat (STERN) 
ZxwyeZ7Zexx-Ty 
Ony 7) s(ZEMyeROOT) OMS 
Duy T2e(Zxev-ROOT) SONS 
Det] Tear ym] (O4IT 1 DN TT2) 
IF (OMIT) 104010460305 
104 Ony724/ax1 (DOKIT] eOnIT2) 
IF (0617) 10101010105 
362 Ong T71.Ee21 
105 COnTINuEe 
CALL DEBUG I6HINTCYL) 


RETURN 
Emp 


C60 


en 


SURQQUTIAE INTDATCIT SLT) 


\ 


Crecceresccoccccoces ose ces esos eSSeeesecoooesoocconsenesetecenseesesete 
Coceecesoe SLAZ0ITINE IkTUaT CALLS TrE DESIRED ROUTInE TO 
Cease © CALCULATE MAVELENGTW CEPZADERT Way INTENSITIES. 


COSSSSS Se Ssosseseogsacogsoseseseoeo epecsesasooesseeseooe SSeSeeceecaaseoese 


COs=S5 449/ 48C0(26s290) oALRay (S06) 
I 0€92°TS(6050) E640 e200) oEB8aS€ (200) 
2 eG42av (500) s8ED 1202209) eTIILECI2) 
3 sv% (200) oKL4n(206) oKRAY (506) 
@ sZ8@av6500) } 

TARAAYSs ALK (200) oAL Prta (46) 
1 s@ETA¢se) eDEv«2ce) 004¢200) 
2 s&Lanc2ee) oEL anal (26) o€Mal (20) 
3 064«(200) oGarma (48) oREF (20) 
@ eSUPRay (206) 0 SU"S02 (206) oX(40) 
S evtee. o¥K (200) oZ(ee)- 
—~\ © oF8a¢2@) 

Coweny sl aaRays IFLAG (0) 

@1S~ac1S@) 03S0a8 (20) oJF (208) 
OM 1726) 0381200) .osu(200) 
Rey (200) 0 J508 (50) oKBn0 (900) 
@n? (200) e<TY? (200) olan? (Sa) 
@LREFL (200) oL TYPE (20) oMaT (200) 

onmaY (28) ommmay (200) 

Compan 79 a6 3457 al oAL™ oALf 
oBE® sfEQ eveRF e01S eOMTn 
GA eG6ta 64a oGanFr ePER 
eSu~Ga = g KY pate eX0OL0 {| »XkREF 
eYndy eYCLO oTREF - oYTRan eZI 
oZupe eZ0LD oZuEF oZTRan 3 

Cow sus yranees Frax olmax oleow 
etc>2aav se IC7any elt oIgay oJMak 
OKYAT oKWLE oKbOw oMnan othOu 


IuvESEa £oF™Max 
OIvEASION TS1G6(200) 
Ceeesetece GO 10 NESIRED Call. 


GO 1061600290 030004002500 06000708 


YOO CALL INTENS (AL Ag F MaKe TSIGeLI) 
60 TO Iie 

200 CALL INTEM2(XLaweFmAKeTSTGoL I) 
60 TO lie 


- 900 CALL SATENZ(xLaneFMAXeTSIGeLT) 


Go TO 11e 

600 Cart INTENG (KL anteFmaxe TSIGOLI) 
60 To lle . 

Ste CACU INTENS (alan Fane TSTG LT) 
GO YO Ie 

eee C TATEMG (XL Ag FHAR GO TSTG LY) 
GO To ite = 

700 Cath IATENT(XLAMSFMAKS TSIGeLT) 
Go TO r1e 

MOC CALL INTE) (ELAM FHAKeTSTGoLT) 
Go To ire 

O09 CALL INTEMO(XLAMSFMARSTSIGSLTD 

110 OO 36x FatoFurx 


ESA ILT oF eTSTG(F) ! 
yer Coatinut 
REYURN 
Em 
Col 


000059003017 


o®turay (See) 
0£S0% (20.290) 
oV¥Mal (20) 
eo YRAY (500) 
oBEK(200) 
oEK (200) 
oEMERGY (299) 
oe SORMAT (26) 
oXK (200) 
02K (206) 
o IPCINT(58) ° 
oJFL (200) ° 
o JP (208) e { 
oKFLAG(200) ° ‘i 
eL ANGLE (26) e. \ 
oMFL (20) : 
0 ALRF RE ® ‘ 
002 e&PS r) 
° 
oXTRAn evi ° 
° 
° 
o INOW oJOLD 
oMSOR 
3 
z 
h i 


Suenout [ME INTELL (KeDH1T) 


Cocco neccsccccceasscsscesecees co scecsessssoescecoososesesessoneseccoseer 


eo IKTELL FINDS THE INTERSECTION DISTAKCE OF A RAY WITH 
© AW ELLICSOID. SEE eRITUP FOR EuURTIONS USED. 


COSESSS OS 00000 S805 OOS SOO SOOEOSOSSESSSSOSOS OSES SOOS ONS OSOOSE POS GOCCOO SS 


Coun 7697 a8C0 (200200) eALRSY 1509) et ntY (SCO) 
I eCOORDS (6050) 0£ 1900200) e£UsSE (200) eES34 (200200) 
2 eGaRAy(S06) oRED(262200) oTITLECI2) ovmMaT(20) 
3 ev0L (200) oXLan(2e00) oxRal 1506) evHaY (Seo) 
& e7fAvtS00) 
COmwon sARRAYSs ALK(200) eAL Pra (46) oBEK (260) 
1 oMFTa14@) eDEV (200) DK 1200) eEK(200) 
2 e€La™(200) eELAmat (20) eEMal (20) eERENSY (200) 
3 eGax(200) oGANMA (40) oREF (20) oSUimat i260) 
& oStmRay (260) eSUMSOR (200) oX(40) oX%& (200) 
P 5S ev(46) e¥K{200) 0Z(46) oZ% (200) 
/ 6 e7L 4120) 
| Comes 1 APRAY/ IFLAG (46) oe IPOIKT(SO) ° 
| @ISWAX (S60) eISOR(2e) oJF (206) 9SFLI200) e 
; PIL 17061 03(200) ‘o Ju (200) oJP (200) e 
f @JQav (200) 0 JSOR (50) X50 (900) oKFLAG (200) ° 
@XP (206) eKTVP (200) olan? (50) eLakGLE (20) ° 
, . @LREFL (200) eL TYPE (20) otal (200) oFL (20) e 
1 \ eomaTl (20) emuUmRay (2060) 
: Cowvotss/PARL4S/ aL oALn eALR ° ALR ry 
‘ CbEN GEN oBERF eD1S oDMI is odZ ofPS oe 
' °6a eGan e6aR oGanF oPER e :, 
| eSum2a 8g XI oAnOw eX0L0 oXREF oXTRAN evI oe 
i \ evn eYOLO eo YREF eYTRAN eZ! ° 
! : OZuC a eZOLO oZREF oZlTaan = 
e\e COws0Ns TP anaes Fmax o Imax olwOw r) 
i \ @®10>cay eIOPRuT 8 oIh eJRatT 0 IMEX o INCE eJOLO 
\ @KIAT “okhax oKNOd oMeax ohare eKSOR ww 


=X 


162 


Cesesececce CALCULATE COEFFICIENTS REEDED 1m EQUATIONS» 
Se beeen vee 1 DEC ROE PER ae ae 


INTEGER FoFmMAK 


FV arms cee ne § ana ney gieme ad cm NT h 


OEL VevIovK(n) 

OFL 72927 L-ZK (K) EK (K) 

OScsEK (Kx) e#2 

Ce] .@<(DELKOO2ZeNELY®#2) /DK (K) -DELZ60°2/8S0 
Bze ( (AL ONEL KOBE *DELY) /0K (K) eCGA®DELZB/ASO) 
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COwMOS /HELCC¥/ Ao ABAR +B 9B 0B 1P oBPoCAPoDEL T= oDELX eDELY0CELZ oCul Ty 

OPS 12 oREV—R~O eS THe THe TH] oSCHKS eRKS oCSTu gShTmoACSTNGRSATH BS Tote 

COST ANS THZ00THZ : 

Cait Tors 

COATIAUE 

Ceeser*ece CALCULATE PaRaMETEAS REEDED 

EReSTweat-CELZeRCSIN 

YRaSTMOPRE<CELY-ASNTN 

ZAsSTHeGA-CEL2—-CaPpetin ~ 

FusQzxnee2evRee2e7Ree2 « 

DOws (AL@ RR oEE CYA GAZA) ; / 

OE sF S702,C°ABAR . ( : 

IF (ABS (OE) 02. T2280] °ABAR) RETURK : 

O€sCE/cow F : 

O68 (1 .-CcoSTane©2) /DDm ‘ » > : 

OELSTH20€/(1,°SCAT (1 e°CE*06)) ps. © . \ 
Ceorsereee CALCWATE DISTANCE TO INTERSECTION. 
7 "STusS Tne ST" , : 
Coeseessce tooo stm “ANGLE TH, ; ide: “ ‘ 

Cait THCFS ‘ , ; ; 


6oro1 |. * 
Ex te 
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SUSROUTIAE REF OAT 

i Concsctcossorccoccssescsssescoosceteseocccocoesocosesecose GOS SOD COCEEte 
Ceeccstese AEFLAT READS REFLECTICH Oa Ta, ‘ 

C@seeeesces 8000800 008600886880 008CESSOCORREOSS -« SCetseceesetoecoovecesesese 


Common 7697 SBCO(200200) = sat mavisea) eBERAY (509) 
I eCOCGRDS (6050) 0 1460266) e€BaSE (256) 0ES07 (200200) 
2 eGaMay (590) oRED1(20 02061 eTIVLE C12) ovmaTi2ze) 
3 evOL (200) oXLAnt290) oXRay (506) oYRAY (Se@) 
@ e%tav¢See) 
COm"ON sARRAYS/ ALK(206) oAL Ona (49) eBEK (290) 
I ePETA 146) eOEv (206) o0K (2¢0) EK (200) 
2 c€L4"( 200) oEFL "aT (26) eEMat 20) oERERGY (298) 
3 0G4K(200) oGAn~s (46) oRES (20) oS02"at (26) 
© eSumaay (206) oSU*SOR (266) oX (66) oXK (206) 
5S ov(66) e¥K(268) 0Z (48) oZ&(260) 
@ eZLawi2@) 
COeruQets 1 ARRAY / IFL66(e0) |» 3POaNT (50) ° 

| @lSeax ise) ©1S50R(20) oF (200) \eoSFL (200) ° 
OSL (209) oF11280) oJut200) jo SP (200) ° 
@JRav i200) oJSOR (S56) eX870 (996) oKFLAG(208) ° 

{ oxp (200) eaTyP(ze0) eLanh tse) oLAAGLE(200 = sg 

: *LOEFLIZ00) «= oL TYPE (26) > omaT(200) oMFL (20) ° 

t onmali2@) oRURAY (2068) 
COmnQasPaRans/ aL , 0fUn eALn oALRF oE ° 
e0En 0SEa oBtaF e01S vOntn eDZ * oEPS oo 
eca Gans eGaa geane PER ° ° 

eSupaa =o K1 oinOw 20.0 oXREF okTRAn evi ° 
e¥noa eYOLD oVREF oYTaag ry 4 | r) 
oZupd eZ0L0 oZREF oZtaan 
Conupus 1 FP cRARs Fmax olmax - elnow e ; 

‘ ; ®lomasy ,107RKT 869 §R olcay oJmax oJnOe eJ0LD 6 
eX int oKMAX othOu oman ohadw eluSoa 


INTESER F oF max 
Ol~cKS10m LTAPLE (10) TEMP C30) ; 
OATAILTASLES6"TABLE e6hCONST OGOIFFUSSGHT LACT oGswaVCaL, 
Che AVIAP SC 4164 » - ' 

! Ceeecccose PAINT PAGE HEADING. 
‘ Print le 

10 FOR4AT(149053X0 ISHMEFLECTICN OaTAs/11m REFLECTIONS GR» 
in PONT YPE/ 11 et MNQER/) 
{ : Coseccesce AAD OME GEFLECTION Oata Caro, 

109 READ 11 ot ol Aare. TEMP 


i 11 FOQ4AT (16010086 010F6.4) 
Cocevooces REFLECTION MUWBER OF ZERO INDICATES ERD OF REFLECTION pata, 
1F (1911001100193 : 
| . 108 List-l6 . 
' Ceecsctece STOIE REFLECTION Data, 


{ BO 103 w=i,le 
: 103 MEDILL oh) =TERO 14) .) 
°° OETERMImE AEFLECTION TYPE, : 
BO 106 17210168 ‘ 
CT vee tL) LT 
IF (L&APELTASLE (LT)? 10601050106 \ 
106 CorTINEe 
Coeccococe REFLECTICN TYPZ WOT Im TABLE, 
Call ERPant (g-acrpat) * 
Stae 
mE 105 IF (LT-5) 10601070106 
306 IF CLT —6) 22207920131 
UT PRINT 13eLeLAAnMEcTERS . 
13 FOMAATCIM 011004X046010F 10.6) 


GO TO 169 
107 PRINT 160l ol APE 
16 FORTAT Ce 0110. 4Ke86) 
GO 10 169. : . . 
Ceocesseose CrarraTEe REFLECTION FAACTIONS aS 2 FUACTICN OF WavELEnctTN, 
2127 Cece MEF IWR) ; ; : 
PRINT 150L cL MA<Es (REDILL OF) oF 2] of wag? 
! wn BS FOR~ATIIN 110 c4eAGs U1 e203 1 CFIO,4)) 
6 10 169 a - 
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Cececcsoes REFIAP CALLS THE CESIGED GCUTIAL FOR GEKERATIAG 


Cogesoccce AFFLECTION FRACTICAS £5 A FUNCTION CF wivELERGTIN, 
Cocoponcocesconcocoococovesovecvesseccsosnccoooscosessscssesccccoconence 


Comvsa /69/ ABCOC260266) eal Rcv(See) oBERAT(596) 
1 eCOCRES16050) 0E 6460200) oEba*z (290) eESC#62002¢2) 
2 shARAy (566) eRED(28 0200) oTITL E612) eweal coc) 

3 ov (200) oR 46206) oXRAY (S68) -aYPAT(Se@s).. 
& o®8314500) . 

‘ Cowens sJARGAYS/ ALK(206) oALPrA(6@) oBEa (200) 

1 e@€7T4¢4@) efEv¢200) oDe (269) oEM (200) 

2 eFL4™6200) eEL ara i2e) oEMal 26) oEREGy (200) 
3 2G64K(206) oGS0A (46) oREF (20) oSGle £1 (29; 
@ eSUPRayv (206) oStS2-2(266) oX(46) oXK (206) 

S evcse) o¥K (206) oztee) oZ& (206) 

6 e%tanc2e) 

COs=0u,/,sRRAys IFLAGt C8) el POinT (56) ° 
@1Ssax¢e) elsonc2er oF (296) oSFt (266) ° 
JL (206) ott206) oJh (266) oP (266) e 
Laavizees oyscaSe) eKGnc (990) oXFLAG (200) ° 

(206) oKTYP (208) elem (S50) eL BAGLE (26) Py 
‘OUREFLI(Z00) «= oL TYPE C29) omal 6206) oMFL ize) ° 
@eheaT (25) oh IPALY (260) : 

Common sPahe4Ss aL oAth o42 ofl SF 3 e 
*8€a e€Ea oBERF e01S eDerya 002 of PS ° 
Ga e€as Gea oGanF ePER ® 
Suna ox! oXnCw eXCLc oRREF oRiaan eV e 
eTNQe eYCLo eo YREF oVTRAn ry $j e 
olan eZCr.0 eZRtF oZTlin = =- 

CO~906 71 Pans Frex olwax elkte e 
@10+a2y ICHRNT 86g hk oIgay o SUE owacu eJLD eo 


e@xiar oth sr ohne ohmax ohkiie e&SOR 
INTFGEQ F oFvax 
DIMEASICA TST6(2e@) bn 
oa 106188 02600306600 05800660070900C0 0906! oil. 
100 COLL REFLIL CAL AMormaxelSIG e0u*) 


60 To 116 
200 CALL REFLIZCKLA™eFraxclSTs eDun) 
60 TO ile = 
300 Cate SEFLIZ(IL Ar FRA ST STC DUM) 
GO TO ile 
400 CALL TEFLIG (XLANSF MAK STSTG oDumty 
* 6€© TO ile 
SOO CALL REFLIS(XLAWeF wARLTSTS eDun) 
60 TO ile 
600 CALL REFLIGCAL AM FRAC TSTS eDun) 
60 TO ike 
yee Cait REFLITCALAM Frakes ISTE eDun) 
60, TO 116 ° 
B00 CALL REFLIOIKLamreFHansTSTG — oOun) 
6O TO ile 


908 CALL REFLIOCAL Ae FRAXeTSTG eOuny 
116 00 1601 Fol Fran 

ES SLL oF DOTSTGCE) 
$01 CoaTIWe 

RETURY 
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oe 
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Ceecoccccccccecoosss ecceeeseeoooece ceacescedecocesoqsoeroososess eesseeto 


Ceseeercse QAEFLCT CALCLUATES THE DIAECTIONAL COSIAES OF Trt 
Ceeccesces AEFLECTED RAVe 
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Cecesesece CALC 


Conn. 1594 a8CO(200230) eataay (Seo) eBEAAT (500) 
I efCars (6256) o£ (6022068) e€GaSE (200) e£SOR (700200) 
2 s@saay (560) onL0 (200200) eTITLE (329 owmal (20) 

3 ev (200) og ant200) ox2av(50e) eveaY (500) 
@ 07287(580) ; 
Comman saRRAYS/ ALK(2C0) oA pra (60) 28EK (230) 

1 emeTatse) eDEV (200) epa(200) EK (290) 

2 c€L 491206) e€Larat(2o) eEmal (20) eERERGY (206) 
3 9cex(200) 2GAewh (48) oREF (20) eSORMAT (20) 
@ osumaay (253) oSU"S58 (200) oX(e0) 23K (200) 

5 e¥(4e) ovK (Zee) 02160) 20241200) 

© oF 29120) 

CowedusTaRRavs TFL A648), o SPOINT (58) PY 
o1Srr2(50) e15ca(2e) oJ (200) oJFL (200) ° 
oN (700) ott 200) evn (200) oP (200) e 
oxlav (200) 0 JS0258) | 258™D (900) o8FL26(200) e 
oap (29%) eatyr(zee) oL as? (50) el A GLE (20) e 
@LReFL (200) ei tyPze (23) - o™aT( 200) oh. 620) e 
onmar(ze) et tRAV (200) 

COP~OR /PARAKS/ au oAL™ ad oALRF PL P 
8s e2ER e6E FF. oDIS Py) Cl 20Z 2e€PS le 
oGa oGas 2622 oGawr oPER 
eSueaa oX3 oXndw eXCLo eRREF oXTAAN ev! ) 
eorugd evVQD eVREF oVTRan oZ3 ® 

Tene eZtL0 oZeeF eztaaw - 

/ Comwoustesasns fmak oSeak o Inu ry 
/ efoecay eIC7=8T e1a oxymaY oak 2 JNO eT 6 
euyat eax ekadeu orem ax ote eMSOR 


INTEGER F oF PAR 
Ceeeeetoce caccuat 


E COEFFICIENT NEECEO. 


STwP=-2.8¢ (AL CAL boBEC3EmeGSOG4R) 


BL abt cB THO Clin 


BENsBE oBTOOCFES 


GameGschT Poca 
RECTIO“AL CCSIRES. 
CALL CHECK (OREFLCT eal A oBER 0SAK) 
CALL DEBUSIGMREFLCT) | 


Ceessetees Cr€c< oL 


GE TURN 
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SuaacuT ihe REFREC 
Coscecececoncces se negocccccecccoscencccsccoococonorosorcooooooosc ee ee’ 


Ceecceetsee REFAA caccuLate THe DIRECTIO«aL COSISES OF IrK€ 


Ceeescesce RE FRACTEOCTRANSWITTED) MATe 
Covrcccascencsnnscessscacccscocooosscoerscesres ee ececconcesssesece 


Cowon /69/ aBCO(200266) eaLnay(50e) eGERAT (500) 
1 eCC0<nS (6058) oF (000206) e£0aS£ 6200) 2£SC4(205200) 
2 efanay (506) eRED (260200) eTITLECI2) eveaTt (20) 

% ewct 1200) eRLAM(2C5) eXRay (500) eYRAT (500) 

@ e7Rav¢50e) 

Coer~Om sARRAYS/ aLx(260) AL Ora (e6) eBEK (200) 

1 onmeTacse) eoeuczee) eBK (200) eEK(209) 

2 e€LaNc200; eEL Amat (26) ofmaTc2e) eEKERGY (200) 

3 eeaxcze0) eGarma (40) eREF (20) e@ScaMat (20) 

& eSepmnay (206) 2 SUMS.Q4 (206) X40) oXK(208) 

S e¥(46) e¥K(200) oztee) oZK(200) 

© eanc2c) : 

CownonsTARRAY/ =. 1FLAG (40) eIPOINT(S9) ‘ 
e1swan tse) o190a120) eJF (206) eJFL (200) e 
1200) ot 1200) oJ (2004 om (200) - e 
egmav (206) 0 gS0R (50) 0X80 (900) eKFLAG (200) ; 
exp (206) exTvP (200) oLacw? (Se) aL anGLe (20) ° 
@LREFL (200) eLTYPECZe) oMat(200) eMei (20) . 
enwaT (20) enamemnay (206) ' 

Cows0n 7? 28 ASS aL eALM oAtR eo ALRE ev€ °. 
osEn eRtR edtRF edIS eOmter 002 . ths e 
oss Gan Gan oGaaF oPER e 
@eSuvaa eX! okra excLd oXREF oATAan evi e 
eoYe.td eYOuo eVREF e¥Taan ry a ° 
elusw eZ0L0 ofktF ezTaam - 

Céeugus 1 PARaMs Fuak o Tmax o INOW ° : 
e1Ne2ay efOPKHT IA e1nav 9 JMAR euhOe eOL0 8 
@nqaT @kPaxk eKnOw olmak ghee euSOR 

IATEGCER F oF MAL 

REaL LCOTN 

1=1.0u \ 

Caecssecee SET ImDICES FCR MATERIALS On EACH SIDE OF BOUNDARY. 

MhewaT¢JOLC) 

e2e~w7aT ( prow) 

cCeesceseed FIND AATIO OF REFRACTIVE INDICES. 
ATepsAer cml) /REF (M2) 


eoccee CALCULATE COEFFICIENTS REEDED: , 
102 LecTassle AL meSECBERCGACGAN 
Ceesccccce IF COSIAE OF ANGLE OF INCIOENCE LESS THAN TeE~6e 60 TO Ite 
IF (ABS (LOOTHI <1 o£ -6) 10101020303 
@ecee SET DIRECTIONAL COSINES FOR: TRANSPITTEO RAYe 
fel AL QF eAL ¢ : 
O€SFaaE 
Gast =sGa : 
PEaz1~-6 \ 
* Go 70 1¢7 
163 soe reas sec (ieS-aTemoez) /((aTHPee2) © (LOOTNOSZ)) 
"YR (TERM) 10402650305 
204 PEaSE 6 
Ceceesseces IF OEFLECTION pUMBER GREATER Tuan 16¢ GO TO 109¢ 
Ie CLREFL (ENDED A 1O) LOT OT OID. 
T12 ALSE ORL . 
' BEQFeE ; 2 
Gaat=Ga 
Go vO 109 
Coeesesece FRESNEL REFLECTIONS - 
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105 BlaPag THP OL DOTAS (~) .9eSCKT (GTENMS ) 
Ceeovoecce CALCULATE COSINES OF REFRACTED Rave 
ALAS SATHP OAL CBI MP OALN 
BEEF eATuPege oATHP OREN ‘? 
gBAF SA THD OG ACH TMPOGAN 
CALl CHECK (OHAEFRAC. ALRF o DERE sGARF) 
106 1F (LEEFL (40m) -1091080308,109 
Cooceecces FIKD PERCENT OF RAY TO BE TRAASHITTED. 
106 CALL “PERCNT 
GO TO 107 
Coecceeoce GET REFLECTION TYPZILTle 
109 LeLAEFe (xnOw) 
LTsL TYPE CL) - 
_ULsL-le 
Coecceteoce AAANCH ON REFLECTION TYPE, 
- 60 1061000200 0 300046000599.600) LT. 
Ceoccceeooce TaARLE LOCKUP REFLECTION. 
2100 Catt TaSReF tLe» 
Go TO 107 
Ceecee®sce CCOAKSTANT REFLECTION. 
2ee Cart CCAAEF (LL) 
60 10 107 
Ceoceeeecce OIFFUSE REFLECTION, 
300 CaLL OFFAEFILL) 
60 10 107 
Cececsecce FUACTION REFLECTIGNI PERCENTa} .00£9(1 00°CCS(TM)DD 
400 Cart FUAZEF ILL) 
GO 70 1067 
See Cart wsaveEF ices 
400 OO 113 Fu) oFmax 
EUL~OXoF PSREDILL oF) OE LT oF} . 
ECL oF BEC LOFICECIMAXSF) 
213 Costiace 
GO TO 310° 
107 DO 112. F=leF max 
E014426F) 5 (1 PER) OE CT oF) 
ECT oF DePEROE (TSF) . 
212 Coatiace 
Ceeeeeecce STORE COSIMES OF TRANSMITTED Rav, 
120 AL O-~4A(}) caL AF 
BETACT) eBERF 
GAr+h 13) 2GARF 
CALL DEBUG (6HREFRAC) 
RETR 
Eup 
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CCOSSS ESF OSCOSECEOENSLOOOSESSOCONSERESCE 1EE00088 SSSSCOCSSOESSESCE 

COs=t0m 769/ 48C0(280288) oAd RAY(SOC) ebERAY(ISE6) 

1 sCO02NS 16.50) eof (680200) e€BASE(200) | o£ SON (2ee200) 
2 eGamay (500) eoREO(ZE e200) eTIVLECI2) otmali26) 

3 svOL (200) oXL AH (200) oxRay¢S¢s) ° eVkAT (See) 
@ o788Y(500) 

Com"Ou S4RQAYS/ AL K&(206) oh Prd (469) oBEn (290) 

1 e@€TAc4e@) oDEVI200) o0< (200) ofS (200) 

2 of LA5(200) eELAnati2e) o€at (26) eERERGY 1200) 
3 06641200) oGares (40) oREF 620) sSCRMAT(2¢3 
4 eSrmavi2ee) «6s SU*SGaI 290) oXlse) XK (200) 

" 9¥(40) \e¥u (200) oZ(4e) oZK (200% 

6 s7Lant2e) \ 

CO="Ons TAARAY/ FFLAS (46) olPCINT«S8») =, 

@1S<2x150) © »02502(20) oJ (263) oFLIZEG) je 
ul (700) oPt(200) om (208) 9571200) _ 
@mMay (700) 0 JS (50) 088A 19900) eMFL46(2¢9) ° 
On? ¢2700) ealyP (Zee) ols (Se) ol e867 € (20) ° 
@ULREFL (200) oL TYPE (26) omal (269) : oMFL (20) ry 
omnatiZe) ohwetavi2ee) 

COON PaRanss au ome oALA oAL SF RE ° 
°CEn e8ER estar 01S eDu}a ebZ o) 9FPS ° 
eGa . eG e640 ous oPER ° 
OSum2a okl ozatu oKQ.D oXAGF eXTaan v¥i ° 
CVs eYCLD oVaar e¥Taag eZ] r) 
e2aoa oft otkee oLTtay 

Cor~ons 1F £527 Fear = clMAX - »ginOu ° » 
*lopmay .iteael 8 of@ olmay osmak o JuOu oJ0LD » 
OuTNT = _ g KMAK Ot o fe otSon 


IMTEGES F oFman : 


COP~ON 74 CON Ae BARB eB1 0€ 1P cP eCAPs DEL Tete DEL 2 oDELYOCEL 7 oDU1T» 
_ SPSTT BEV eRe DoSTH9 THe Ter] oSCKAS -RESp CST ne SATNOKCST Ho RSMTN BST oto 


“CosTame InZoDtnz 
1Cse@ 
1F ( (BABAR) © (31-4822) ) 100. 1000202 
200 Cart FIepe 
GO TO ile 
101 IF (18-2848) SBP eDEL TH) 11001040166 
166 952" ag = 
Toon, = Tes 
pag 
TogeTos 
Taga Tit} 
B4e9) , 
Sercsl1° 
12 CORT INUE 
, OTe ( tte 14) 7 (BP 84) OOP 
THastn 7 ry 
. 8453 ° 
Lr 
Totes Titus) Tog t ; 
. JE CCT THD CT THG)) 13018018 
16 ConTinge 
Tits . 5° ( TH] oTH6) 
13 CONT Tae 
Cau ScaLc 
@n{2-4842) /(50-4BAR) 
IF ¢R) 10602660198 . 
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sa 


wee 


Feed oo — —e 


19 COATISLE 
WF (2eLT. O05) GO TO 14 


IF (aCereP) Lae2Ge20 
2a COAKTINUE 
LF (25S (Tal oTr6) oSTec01) 6S TS 44 


332 CONT INCE 
Tra 145 
8225 
bP 2-5? 

RETURN 

14 COST INE 
IF (20ea,P.LT.0) GO TO 15 
14] oT 
Bisa . 
8}9=92 
GO TO 16 

15 COAT INUE 
TaqaTat 

16 COAT INE 
IcetCeo} 


IFCICoL TSO) GS TO 12 
Patel 17e Top hebPoTrtLed 10S 1P 0Th40869d4PoTHSeBSeB SP 0 ING 
iT FoavaT(o~ REXS 17efl 12E10033) 
60 TO 132 ? 
106 Cait FINOR 
Toy eto 
84<8 
BlpszReP 
Tra THS 
Beas 
BP arse 
Carl FINOR 
le COATINE 
RETURN 
Eup 
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Coes. “69/ aBpcCO(200205) eataeyi500) eBeaavi¢scec) 
1 eCOGRTS 16050) ef (400206) efcaSt (260) e£SO%1262205) 
2 eGanay (Seo) eREO (260200) eT) WE C12) eveal (20) 
3 ov (206) eXLam(2G0) exRav (S00) eYraY(csoo) 
@ 9 ZAAV(59068) ; 
Cows0n JARRAYS/ ALK (260) ekLOratso) eBEK (260) 
5 gRETA4O) ebEvi2es) eon (265) oEK C256) 
2 ettanc200) oF aral(20) e€Mar (20) eE&CRGY (200) 
3 9eaKkt200) 2Garea 40) eREF (20) esoamaT (20) 
@ .sueemay (200) eSUMSSAIZO0) oXte0) oe XK (2065) 
S ev(s6) e¥K (205) oZtae) oZK (260) 
© e7Lawi2o) 

Cow40~s 1] ARRAYS TFLAG (40) elPCiey (507 . 
@1S5van(So) 01568120) oUF (200) e JFL (206) . 
oF. 6200) om (200) 0 Ju (200) esp (2c) e 
@jeay (200) eJSCR (SO) eKBrd (900) eKFLAG(250) e 
exe (705) extyr 1206) et an? (50) eLAnGLF. (20) e 
@LRFEFL (200) eL TYPE (20) oMalT (200) eMFL (26) e 
oh (20) eheMRAY (200) 

COw“ONsPARAYS/ at eALN oAtf eALRF eRe e 
@oSEa eBER eBERF eOrs eDmy® 002 of PS e 
@6a eGAaN 96aR eGanF eFEa , 

@SuwRa eX) eXnOa eic.0 e AREF oX TRAN evi e 
eYntd eYOLD e Yer oe YTREN ZI t 
ezacd Py 4G) 9 ZREF eZTaan 

COvOns 1 PAK ees Fuad olmak . glhoe . 
elerrey 10°RNT ola eyRay oJMAX esne Pe, 
@K1nY ekMad eKNOe eo RM AK omNOu ehSOR 

INTEGER © oF “AK 

Cowsonsezsi 2 FRAC elOEare 

100 CO 20) FaleF Man 
IF ce (1oF) 20) o£ 6) SESASELF)) 10201020102 
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elcoazy 41 °° aNT ela olay osmak , oJ0u eJOLD 
exist oC Ax 2fROu orm Ak oWOw ensca 


IATSSES F,Fean 
CGe ons Sul COMsRNOP AR (12) 
Coeveteeoe Silsi 1s ul Cala 
AGC O<)2aQrogaty) 
Ve (4) 520°222(2) 
26 (43 53N7724(3) 
Rle2n7F 0214) 
Y1s8.5>42(5) 
Z1sAN3=42(6) 
ZF wAxDraa(y) 
SF x7F <7 (4) 
Mhn2.9° (21-20 (<)) OSFo(Z1-7Z« (a) ) 002 
BHz( T1191) ) CO2e LY LeYK(K)) OO2 
CCz4.0°SFO(Z21-2a1K)) 
Coceecetece SFE SF CATA DESCAIBES ad ELLIPSOIC. 
IF (40-92! 191019101603 
VOL LF (83200) 102010920163 
Coecceeteee Sala COES KOT DESCRIBE ELLIPSOID. 
303 PANT 120K 
12 FOa~AaTC 141 o21rDaTa FOR BOULDARY MOeelbe 
O27 DOES WOT FIT an ELLIPSC1O.) 
tine 
162 CCaAYInce 
Ceoeeeeeooe CaLCULATE COCZFFICIERTS REECED,. 
WStz (K1-1K (4) ) COZe (VI VK (K) 1,002 
ZSTaZlezZa (a) 
RzSIAT ((Z7ST<SF) 0020050) ©SIEN(1.00SF) 
Ceeeseesoo CarcUATE 6. 
Ext) S(ZST9SF R17 (4 .0eSOs(ZSTeSF)) 


Ceeesetece CALCULATE acee2, 
Om (4) SEK (IK) ©92- (EM (K) @SF) 002 
Ceaeeetooe PAINT 1kPUT DATAe 
Past LL oR oLREFLIK) 0 (5X>P A(R) oheye7) 
Tl FOcwat¢in TS elGe2K oGrELL IPS 026k 0 7F 8.3) 
RE TURE 
Eup 
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SU22OUTINE SETet! (a) 


Céeoceseteces SPHSSHSHSVSHSHSSOSHSSSHSHSSESSHSSHSSESSSEHOSEHGEHSHSSHESENHOCEEEOED BLE 
Céerectcacs SESCOHOEASOODOOSOE SHES OOOCOP SHES EHOCOCOO SOE OO SE SCC EDEOCOHOE ROO 

COeetn 69/ £369(250200) escaey (S05) wsb5h41555) 

Y eCOl3T$i6e58) = sE Lee 2ti) e€sece (205) or ScsteC ess 

2 e642hv(5C635 eRE JI 200255 eTITLE ti 2° evPst (27) 

3 ev (209) eRL S255) eXaev (Stay oVae¥ (505) 

& e7=av(509) 

COm~d% sBARAYS/ = ALKI200) ef 2-4 (69) eBE< (260) 

1 eAETA(40) eDEvV(200) o0<(296) eEK (258) 

2 0€L 41200) e€L AMAT (29) oF" 2T(20) oLENERGY (290) 

3 e€2K (268) oGarrk (60) eREF (20) oSea"4T (26) 

& oSurzay (206) oSuvStat2se) ox(40) eR C250) 

S svcad) e¥K(7 28) 0Z(eS) eZK (209) 

© ot 120) 

COvr"-/1TaRReys IFLAG (40) oe TPCONT(SS) ° 

@1Swer 559) 91So0a¢29) 9JF (295) eJFL (208; * 
ON, 1750) 0 4(209) eh (200) odP (266) ‘* 
eJaav(2e@) 0 JSOR(56) oMan0 (G00) oKFLAG(250) ® 
@KP (700) oKTYP (206) ols~P (Se) eLANGLE (20) ° 
LAF FL (200) eL TYPE (20) om2T (206) oMFL (26) ® 
oneat (25) eMUMARY (220) 

Cowon sPe,ae~ss aL eALN e4LR o SLAF eM ° 
o8Ea etE2 oSERF oOIs oDMIh 002 oEPS Py 
oGs oGin 9642 oGaeF ePEa * 
oSuw2e og KY othe oXCid oXHEF slTR&N eS ® 
eYnr a eYC.0 oVnkF sYTRin Py 4 | * 
orlarda eZC_D oZHEF oZTRAN 

Cowwons {P25 16/ Foan olmak olhOu ° 
OlOprAyY gil Pant otk oTRAY =. 9 IMAX aoe et0Ll of 
@KAT oh SX eKnOe emaAx oWhCa ehSor 


le 


INTEGER Far MRE 
COerOws RSC OMs SNP LRT 2) 

Ka (4) B=%OP LACT) 

YR (<3 2£GN57£2(2) 

ZK (a) BPMDPAS (3) 

ALE (KK) z7ARTPLER (6) 

BE w(K) s§ADE62 (5) 

Ow (K) 22e-TP AR 6) 

PAY. 6 19K el REFL (IK) o (BNCP RAK) gho1 96) 
FOSMZT¢ = elSolGe2Xe6M HELI Xe24X06F 6.3) 
RETURN 

Eno 
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Sus2cuTIZ€ SE 


Ceeeeecsasececessoe 


Ceeecstocccecacsess 


THYP (4) 
SOTSOSSHESOSSCHSEOOODE OS OSDOSHHSOHHSSSOSCE PeseecesSesare 


etcese SeTRYP SETS UP Tek annaYS FCH A MYPENSOLOIDe 


PPPPTTT TTL ar 


Covers 4697 28C09120+209) eRLR4Y(500) eHERAY(590) 
1 eCPL2I7S (6050) «=o E1400200) eE£8aSE (200) eESOR (20.200) 
2 sG8eayvescgs eRED(20020C) eTITLECI2) evwal(20) 

3 evL (760) eXLAM(200) okHav (506) eVRAY (500) 
@ 979270502) 
Corw-3\ saaRaYS/ ax (200) oALPHA(40) oBEK(200) 
1 e8F Tase) eCEvi200) 00K (200) 0EX (200) 
2 oF LA4(200) eELAMAT(20) eERaT (20) efRERGY (200) 
3 eract2co) eGAKMA(40) oREF (20) oSGUMAT (20) 
& eSewaev (205) oSu@SCAR(290) ox(e0) eXK (200) 
5S evis0> eK (200) 07140) ZK (200) 
@ e244 (26) 

COwetnsTaRRays IFLAG(40) eo LOCINT(S50) ° 
@1Svan(Sa) e1SCR(20) eJF (200) oJFL(200) ° 
@9L (799) oJ (200) oJh (200) r 03° (200) ° 
eJaayv(200) eyscaiso) oKB%5 1900) oKFLAG (200) e 
x? (200) ekTYP (200) el aX? (S50) eLaXGLE (20) ° 
@LAFEFLI209) eK TYPE (20) eal (200) oI. (26) ° 
enuaT(76) eKUMRAY (200) 

Co“voxsPaRaMsSs AL eALN oALA oALOF oBE ° 
edEa eRER eSERF eO01S oDM Is 00Z ofPS ° 
eGa , eGar oGanx eGarF oPER ° 
eSumaa =e XI eXhOw eXCLO oXREF © oXTRAN e¥I e 
eTNne eYCLO eYREF eoYTRAN o2i °. 
eZnca eZOLC eo ZREF oZTRAN F 

Coe-cns [PARews Fuax olmax olNOd ° 
elcsaay eIChRNT oleh olRay o-MMaX oJNOW eJOLD o 
exIAT eKMSK eKNOw oMnAX onhOg etSOR 


INTEGER F eFrAN 

Cov“sonsaNDCOM/54N9P 42 (12) 
Ceeceetece STORE INPUT DATAe 

RK (K) ADP AR (1) 

Ya (aK) SENN AR (2) 

ZK (<) 28%7P AR (3) 

Gan (KP =F. SPAR ITPOZKUK) 


X)aAe.0F425(4) 
V1 sANDF 2215) 
Z12980P AR (6) 


Ang  OPGAK(KD O(ZLOZKERKD) 
Ba Clexc(K)) O20 (Y1L2VK(K)) O82 

Ceeeceeece SEE IF DATA DESCRIBES A HYPERBOLOIO. 
IF (AD 1OT 21020102 


102 PAIAT 100K 


10 FOOMAT()HI eZ: 


“HDATA FOR BOUNDARY "Ovel be > 


©32~5CES NOT CESCRIBE A HYPERBOLOIDe) 


stoe 


Ceecceecee CALCULATE COEFFICIENTS NEEDED. 


1Ol DELX2x1<xX(K) 


DELYZYV}LoYKIK) opie: 


WSCADELK®*Z2eDELY°°2 same 


ZSTaZ 197K IK) 


Aesgat((ZSTGAK(K) ) @O2* wSQ) CSIGN(1000ZST) 
Cesececeeteee CALCULATE AND STORE 6. 

EK (K) 2 (ZSTOGAK(K) eRD/ (WSO7 (ZSTOGAK(K) ) 9400) 
Ceeeeeeece CALCULATE AND STOKE A0®2. 

OK () BIEK(K) CGAK(K) ) POZE (A) PO? 
Ceeceeeeee PAINT INPUT DATAs & 


PAINT TleKel 
11 FORMATCIN ol 


RETURN 
Eno 


REFL(K) 0 (BAOPARIN) cN2,07) 
51 16e2K eGret (PERG 244s TF 853) 
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SUAROUTINE SETPAR(K) 


CEOS 08 SOE o OOOH OSEETHOOOOOS SOLES SESH OL OS OES OO OO SESS HOS ESES OSE DES EDEN ESOS 


ce ee SFTPZA SETS UP ThE £ha4£YS FCR &2 PARARCLOID. 
COS SS SSS eCESHHHEHHSHHFOHHSSEHESSHSSEHHS SHH SGAHHHHSHSHHSSHSLOS ESC EE EOOSESS 

Ccoe~08 7697 ABCO(200200) efL RSV (550) eEENAY (505) 
1 eCAOENS (6059) 0 E 1400200) eEvast (200) e£S0%(205200) 
2 eG2RayviSoo) eRED(26 0206) _ oTTILe 12) eval (20) 
3 ev (200) oXL44(200) oXRayvi(S6o) eY&AY(S500) 
& e7RAY(506) 

COMMON /ARRAYS/ ALK{(200) eo ALPHA(40) eBEK (200) 
1 e®ETA(40) eDEVI200) eOK (200) eEK (200) 
< oF L4™¢290) eEL Amat (20) oF al(20) eEXERGY (200) 
3 068% (200) oGAMMA (40) eREF (20) oSCRMAT (20) 
@ eSuvRayi2co) eSUMSOR(2003 = on A040) oXK (200) 
5S e¥(40) 0 YK(200) 0Z(40) 0ZK (200) 
6 e7Lan(20) 

COmM*ON/ 1 ARRAY/ IFL4G6 (40) eo TETINT(SO) ° 
@lSwax(Sg) #1S08(20) o JF (200) eJFL (200) ° 
©RN (700) 0341200) oJ (200) oJP (260) ° 
@JRavi2e0) eJS0R (59) oKB*) (900) oKFLAG (2003 e 
KP (700) eKTYP (200) eLam? (Se) eLaNXGLE (20) ‘6 
@LREFL (700) oL TYPE (20) eomaT(200) oMFL (20) ° 
eoneaT (29) eMUHRAY (200) 

CO™“ONsPARLMS/ © al wALN e4ta »SLRF eRe 5 
BEX ¢ e3ER eBERF e01S eC¥Ih ie 4 eEPS ° 
6a oGan Gan oGakF ePER e 
Suma eX! oxk0e eXOLD oXREF eATRAN ev e 
eYKCs eYOLO eo YREF eo YTRAN Zl e 
eZane eZCLO eZREF oZTRAN 

CO~40N7 1PARAKs Fuax olmax «Inow e 
*10P2ay = FTOPANT =o JR oUkay o JMAX sJA0W eJOLD 
eKINT oKPAX oKhOd oMmMax ohkOu ehSOR 


IKTEGER FoFHAK 
CO#40xsRNDCOM/BNDPAR (12) 
Ceeeeeseoe STORE IKPUT OATAL 
RK (K) eAMHPAR CL) 
Ya (%) SANDPAR (2) 
Za (0) wRt4)P AR (3) 
Dx (X) TRANDPAR (4) ZK CK) 
Ceeseeeece PRINT INPUT DATAe 
PRINT 11eKeLREFLIK) 0 (BMDPAR(N) 9MB1 04) gOK(K) 
11 FORMAT (IM 0350169 3X0 SHPARAD 0 242 04F 823 024X oF 603) 
RETURN 
En 
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SyRIDuTPre SETPLA {(K} 


Ceeeeereo see gavbevencecssacansseusecsceesacsdtetecsarstorsNeseecereesss 


cee e SeTFeLa SETS UP a FLane BCU Cade . 
en rrsurcesiectuanguesuesscnoncenswencnonsosecsetreete 27 222 22 ShS00"" 

CowrvGn 7597 BACOLZ0¢e 206) eat raviSo0) edERaY ("Ney 
L ef PCRS 1605C? o£ (4002C0) e€saSE (295) sESOR(20.c00) 

2 ocarev(506) eREOIZGe290) oTITLE (32) ovMaT¢20) 

3 ew (700) otLan(200) ekatay (S00) eo VRaY (500) 

@ 9758V(550) 

Coesse sacRarS/ aL<{200) eAL Pra (40) e8Ex (200) 

1 esr TA(69) eDEvV(2c0) e0a 1260) oEK (200) 

2 e&caar2ce) eErawat¢20) efrati20) eEnERGY (200) 
3 ef ant200) eGev“a led) oREF (20) 9 SORMAT (20) 
& oSmaavi200) oSUNSOR (200) ox(40) oXK (2068) 

S ev(se) o VK (200) 0Z(4@) oZK(200) 

6 ezLsv(20) 

Cow “Ors TRRRAY/ TFL a6 (68) 0 TPOINT(SO) ° 
o1S~ax(Se) e1S0R(26) oUF (2706) 0 JEL (200) ° 
@#59.1200) 0 41290) ewe i200) 0 JP (209) ° 
egnav( 220) e JSR 150) oKB~D (900) oKFLAG(200) ° 
exe (760) exTyP (260) eLan? (56) eLanGLe (20) ° 
oLacrL (209) eL TYPE (20) omat (200) oMFL (20) ° 
enue T( 2c) etureray (200) 

CcwsonsPanays/ AL oAL™ ALR o ALAF o8E e 
enka e8ER eBERF eO1S eDaIn 002 eEPS ° 
eGs eGeN 96aR oGanF oPER r) 
eSuv2a eX? oxnOd eACLO eXREF oXTRAN evi ° 
eCVure eVCtD e YREF oYTRAn ry 2 ° 
Cy ar e210 eZaHEF oZTRAN 

CowatastPscis/ FMAX eIMAK 6. 9 SNOd ° 
elorasy 9 1Or= ST 018 eo TRAY o JAK o JNOW eJOLO 9 
enyat ouwak okxOe ommAk oteOu ehSOR 


IRTEGES Fe Fearn 
COw@Cn sPAOCO¥/ SDP ARTZ) 
Ceecectece STGRE goundaRy DATAs 
ALa (K) zALDPARCL) 
BEa(K)s3N0PAR(2) 
Gan (CM) s=6HP RR 13) 
Catt CHECK (OnSETPLAG ALK (KD BER (K) GAR (KD) 
Ka (<) sPADPAT (4S) 
Ye (<) sAaxpPar(S) 
Za (<) seAnOPAR (6) 
Ceeseseree CaLcucaTe OISTaMcE TO BOUADARY, 
Ox tc) SALE (KD OX (K) OBER (KD OYA CK) OGAK (CD 02K (8) 
Caceeseece IF OISTANMCE LESS TrAn Oa EQUAL TO ZEROe PRINT AND STOP. 
IF CDK (KP DTOLO LOL IO2 
wel PRINT Lh eKeDK(K) 
i FORMAT (1HO 9 T2“30UNDARY nO, cGeTSt HAS A WEGATIVE® 
20» OR ZERO DISTANCE OF p IPEL24) 
S$Tco 
102 COKTINME 
Ceesseesee PAINT BOUNDARY OATA FCH PLAAKEo 
Print LecKLREFL (KD oALK (KD oBEK(K) oGAK (KD 0X (KD 9YRIKD 9 
©2K (<) eOK IK) / 
Te FoauaT( in 9150 1603KeSMPLANE 96F 8030 IZ2Xe2F B03) 
RE TURN ? 
Exp 
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weeane 


SUQROUTINE SETSPH IK) 
Cecccesccesencccescescccesecenes coscceeoes SPeCCODOOD ELC EEODCL CDE DCCCCREDS 


Ceeceetsoe SETSPH SETS UP a SPrERICAL SOUNGAaY. 
Cocccerrecoccccccccecesesconcesecesececccencecesee cee sontcetcecssceces te 


COe“0nN /69/ 4BC0(209200) eALRayisSoo) eBERAY (S00; 

1 eCOORDS (6050) 0€ (400200) eEUASE (200) 0€S50% (20.200) 
2 e€4RAY(500)  oRED(20,200) eTITLE (12) oVPAT (20) 

3 ev0L(2e0) eXL4n(200}) oKRay (S09) eYRAY (S00) 
4 e7R4Y(500)} 

COm~Om saRRAYSs aLK(206) AL Pra (40) BEX (200) 

1 eRETA(4Q) eDEV(200) oDK (206) 0EK (260) 
2 cFLam1200) sELAmaT(2@) e£MaT (20) ENERGY (200) 

3 064K (200) eo Gla (40) eREF (20) oSGamay (2@) 
@ oS'*Rav(200) oe SU"SOR (200) oX(69) ©XK (200) 

S ev(4e) o¥K (206) 02 (40) 02K (200) 
@ o7Lan(2e@) 

Cowxcns1 aRRAYs 1FLAG (40) eIPOINT (50) e 
*1S¥ax(S0) eISOR(2@) 0 JF (200) eo JFL (200) e 
©1200) 0 JM(200) oJN(200) oJP (200) ° 
*JRAY(200) e JSOR(S0} #KBN0 (900) oKFLAG(200) ° 
OXKP (706) eXTYP (200) oLaMP (S59) eL ANGLE (20) e 
@LAEFL (200) eoL TYPE (20) eMaT (209) o9FL (20) ° 
Snuati2e} oMUMRAY (200) 

CO~vOnsPaRAXS/s aL oALN " 9ALR eo ALRF oRE ° 
BEA ofER eBERF e01S oDuin oDZ eEPS e 
e6A eGan Gan oGakFr ePER ° 
eSuera, eX! . oANGw eX0LO eo AREF oXTRAn eVI ° 
©YNQw eYOLo eYREF eYTRAN ZI ° 
e2nou eZOLO oZReF oZTHan 

COunQns [Pan ans Fomax olmax - eI] now ° a 
*lopray .J0fanT oIR okRay o Sax oJNOM eJOLO 6 
OKINT oKMax oKneOw oMmAX o MeO eNSOR 

"8 INTEGER F oFmax 
CO“40N sANDCOM/BNDPAR ¢12) 
Ceeseeccce STORE INPUT DATAs 

AK (K) sRNDPAR (1) 

Yq (©) sBNDPAR (2) 

ZK (%) SANDPAR (3) 

Ox (<) eQMDPAR (4) 


Ceeseseces PRINT INPUT DTA, 
PRINT 11 oKoLREFL IK) oXKCK) oVK(K) oZKIK) SDK (KD 
11 FOQMAT(IM 0150 1602X eGHSPMEDE 0 24K 9 3F Bi 30 32K 0F 8,3) 
RETURN 
Eup 


a 
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SUE2“MUTINE SE TUPP 


Ceeeesese SSS SCHSSSHSOCHSECECEE SSS SOS OOSOLES CEOS OOEOEHEEOCE SOS ECOCEOOCOEOD 


© SURELETISE SETUPP REACS Tré INPUT CATs AmO HAROLES ALL 


eBERAY (S98) 
eESOR (205290) 
eveali2e) 
eVRAY (See) 


oBEx(2¢0) 
EK (200) 
vERERGY (200) 
9 SGRMaT (28) 
oXK (206) 

oZK (200) 


oe IPCINT(58) ° 
e JFL (200) ° 
oJ (200) ° 
eXFLAG (200) ° 
eLAAGLE (26) ) 
oMFL (20) ° 


o ALRF eh ° 
00Z oEPS ® 


° 
oATRAN evil ® 


o Jnow eJOLD » 


Cee 
Cee oe SETUP AND PACALE INLTIALIZAT190%, 
COOOOO8S 0080006 SOOO ES + OHSS OOOS FOOSE OOESS OOS ESOSEO SESH OOOO CECE OEESOROEE 
COw+tn 769/ €2051200290) eALRezy (509) 
| 1 eC 20916050) = E 1400200) #E€8aSE (200) 
2 eGakavi569: eREO1 200200) eT1TLe £12) 
3 ov 1700) oXL4%4(206) eXRav (500) 
@ o7a8av1500) 
CO@“ON saRAavLs = ALKI200) s 1L P48 (48) 
| 1 eat Tareg) eDEvi20G6) 9M" (200) 
2 oF LaMf206) EL Amat (20) oEmatt28), 
3 eax 1200) eGarwA (40) eREFI20) 
& eSu~aayv(206) eo SU"S$5R (200) oX (40) 
5S ov(se) e¥K(206) 0Z (60) 
@ 071 4%(20) 
COe~64, 1aRRAY/ 1FLAG (45) 
©15"4x(69) 01508120) o JF (200) 
ON 1206) 041200) 0 Jh (200) 
@JRaV 6200) 0JS0215C) 0KBA0 (506) 
| @KP (760) eK TYP (200) ela? (SQ) 
@LREFL (200) oL TYPE (20) omal (206) 
enum aT(26) ener: 71200) 
COe0N /PARANS/ ar oALN eALR 
@BEN e6Ea eBERF oD1S oDeryn 
esa o6a% e62a eGarF ePEa 
eSueaa eX eXhOd eXCLD eo RREF 
ervnoe e¥C.9 oe YREF oVTtan eZl 
elnee eZ0LD . ep 2REF eZYAaw. 
Coprons1Pantes Fuax olmwax olhow 
e1locasy .10°RkNT 018 eIRay oJMAX 
@K1KT oKMax oKknOw oMmak oMnOe 
INTEGER Ferman 
Comme, tases FRAC e LOESUG 
Coceceetece LAITIALIZE PARAMETERS, 
IRaveo * 
| ION} 
€Pss1.€-3 
Sw2ase.8 
Imatse 
PayxT 12 
12 FoRvaT(1H1e214GEMERAL InPUT FOR ZaP7/) 


Ceeeceteee AFAD INPUT DalThe 
REar LOsTITLE 
le FOa~aT(1286) 
REAM LleFFaCoDZewayniNne 
31 FOR 44T(6E12.4) 


wAVPAK 


READ 16_F MAXs LOPANT s JOEBUGSIA 


15 Foaq~aT 63160116) 


Ceecceccce IF MULTIPLE WAVELEMGIMS® CALCULATE AMD STORE THEM, 


1F (Feaz]e}) 10901100111 
{ 109 Feaxs] 
GO 10 116 


111 GF s(wayMakouaVel hy) /(FMAXO]) 


WAyF rd aye 
OO 112 FsleF max 
A_as(F) awave 
wAyF suave eOF 
112 COwTInue 
GO To 113 
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Pare ern 


=e 


11e KLaw(}) swayed 


113 Coa; 

If «F 
165 Frac 
107 Coat 


IME 
RAC) 1050105e107 
Bee 


Int 


PRINT 130eTI TLE 
13 FOR4AT (165128677) 

PAINT 14 cFRACeDZ o WAVMIN, BAUMAZ oF MAX, LOFRANTS (OEHUGEIR 
16 FORMATCIH e1BPROULETTE FRACTION=eE12.4/ 


Ceeeeeesoe 


100 Care REFOAT 

eeeceseece READ MATERIAL OATA AND SET UP MATERIAL ARFAYS. 
CALL waTGAaT 

Ceeseeeocees FFaD BOURCARY DATA ANC SET UP BOUADARY LAGAYS. 
Catt S8apo0aT 

Ceeeeesces READ SEGMENT DATA AND SET UP SEG ENT ARNAYS. 
Cart SFGeoaT 

Coceceesose arap STRUCTURE OATA AND SET UP SEGHENT STRUCTURE. 
Catt STRUCT 

Ceeseeeeoes READ LAMP OATA AND SET UP LAMP ARRAYS. 
Cal LawnaT 

Cesceeseeee MF4D SOURCE OATAe 
Care Scanalt 

Ceeceeteoee RELI) TEST POINT OaTA AKO TEST SETUP. 
Cait vesy : 

Ceeesstece GENERATE PRINTER PLOTS OF GEOMETRY. 
Catt GEOPLT 
CALL DEBUG (6NSETUPP) 
RETURA 
Esp 


OISTAMCE BETWEEN CONTIAUITIVE BOUNCARIES® oE1200/ 
MIR LMU WAVELENGTHZ E12.4/ 

MARLON WAVELEAGTNI o£ 1204/ 

NUYBER OF BAVELENGINS=.16/ 

KAY PRINT OPT10%Z016/20" DEBUG PRIAT OPTICNSe16/ 
STAaATIKG RANDOM INTEGEN20116) : 

EAD REFLECTION OATA 
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ao 


ney are Gamma) 


ose 


—_— 


SURGOUTINE SOF Tr 
CO ¢ 00000000000 000000 00 000000800 O08808 5006000000 OOO0COOHOOCOS OR OCRNOEBOOO 
Ceeeeceeoe SOE Tn CALCULATES ST ateD OSIn~ 


@ooe PERE CF MELIX IMTERSTCTICR CALCULATION. 
CO OO0 O60 00 00 0088020 GOODES OEE OS OOOOODODEOHA HOOD OHEOEDEHOS COHEOOHOOHCOEOOD 
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